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Introduction 


Throughout this module you will use computer software called Maxima to 
investigate mathematical concepts and perform calculations that may be 
too laborious or difficult to do by hand. Maxima is a computer algebra 
system (often abbreviated to CAS); that is, a computer program that 
can perform algebraic manipulations, such as expanding brackets, 
factorising expressions and solving equations, as well as performing 
numerical calculations and plotting graphs. It is typical of the type of 
computer software used by mathematicians today. 


This Guide explains how to use Maxima for some of the mathematics in 
MST124. It is designed to be read in conjunction with the main study 
texts and is not meant to be read all at once because its later sections use 
mathematical concepts that you may not have yet met. As you work 
through the main study texts in Books A—D, you will come across sections 
or activities marked with the icon shown here. Some of these ask you to 
study part of this Guide. Other activities ask you to solve a problem using 
Maxima commands that you have previously learned. While the 
mathematical answers to such activities are included within the 
appropriate book, the computer methods needed to solve them are given in 
the section ‘Computer methods for CAS activities in Books A-D’ towards 
the end of this Guide. 


The final section of this Guide contains a summary of all the Maxima 
commands introduced. This should be useful for reference, both as you 
work through the module and afterwards. 


Note that this Guide does not attempt to describe all the features of 
Maxima. If you need more information, please consult the extensive 
documentation provided on the Maxima website at 
maxima.sourceforge.net or use Maxima’s help system, which is 
described in Subsection 2.8 of this Guide. 


The ‘Maxima accessibility guide’ (see page 94) is primarily aimed at those 
who may have difficulty using Maxima because of a disability, and provides 
some advice and suggestions. 


William Schelter (1947-2001) 


Introduction 


A word of warning 


Using a computer algebra system is not a substitute for learning the 
mathematical techniques taught in this module! You need to become 
proficient at using the methods taught so that you can apply them in the 
examination and in your future studies, and also to help you develop your 
mathematical understanding and intuition. 


Maxima is not intelligent — it is just a tool, similar to your calculator. You 
need to use your own mathematical skills to guide it and interpret its 
results critically. 


Conventions 


Here are the various computing terms used in this Guide, and how they 


should be interpreted on Microsoft Windows and Apple Mac computers. “Rather than learning how to solve that, 
shouldn't | be learning how to operate 


software that can solve that?” 


Computer terminology 


Term Microsoft Windows Apple Mac 
Click Click the left mouse/touchpad button Click the mouse/touchpad button 
Right click | Click the right mouse/touchpad button Hold the key while clicking 
the mouse/touchpad button 
Select Either click to select an option, Either click to select an option, 
or select text by dragging the cursor over or select text by dragging the cursor over 
the text while pressing the left mouse button | the text while pressing the mouse button 
Ctr1-q Hold down while pressing (Q) Hold down while pressing 
Etc. Etc. Etc. — but there are exceptions to this 
pattern that we'll mention as they arise. 


On some computer keyboards the key might be labelled : 


or simply 


In this Guide, these keys are referred to as (Enter) and 
respectively. 


, and the key might be labelled 


Sometimes = (Return ) are two separate keys. For such keyboards 
references in this Guide to the[ Enter }key should be interpreted as 
references to the key. 


1 Installing Maxima 


Maxima can be used on many different types of computer, and there are 
several different intermediary systems, known as interfaces or 
front-ends, that provide a means of communication between you, the 
user, and the Maxima system itself. This is illustrated in Figure 1. The 
interfaces enable you to enter commands that are then processed by 
Maxima, and they display the results obtained. 


Interfaces 


Figure 1 Different interfaces for Maxima 


This Guide concentrates on using Maxima on a Microsoft Windows 
personal computer through the wxMaxima interface, which is shown in 
Figure 2. This is the interface that you are recommended to use. (The 
letters ‘wx’ in the name of this interface refer to the software used to 
create it.) 


(¥i2) expand ( (2*x+1) * (5-3*x%2)) 7 
(tol) -6 9-3 x7+10 x45 


Figure 2) The wxMaxima interface 


If you need to use an interface other than wxMaxima, for example, if you 
are not using a Windows computer and you are unable to install 
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Installing Maxima 


wxMaxima, then you should try to achieve all the outcomes shown in this 


Guide using your chosen interface. 


All the Maxima commands mentioned in this Guide can be used with any 
interface to Maxima, unless indicated otherwise. 


If you use a screenreader the basic command-line interface shown in 


‘igure 3 may be the easiest interface for you to use. 


There are more details on choosing an interface, configuring the 
wxMa 
interface in the ‘Maxima accessibility guide’ (see page 94). If you think 


ma interface to meet your needs and using the command-line 


you need this advice, you should read that section now 


Maxima 5.30.0 http: //maxima.sourceforge.net 
lusing Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) 
Distributed under the GNU Public License. See the file COPYING 
Dedicated to the memory of William Schelter 

Ihe function bug_report() provides bug reporting information 
(#42) 


Figure 3. The command-line interface to Maxima 


Before you can use Maxima, you need to install it on your computer. You 


will do this in the next activity. 


Computer activity 1 /nstalling Maxima 


(a) Go to the OU Maxima website at learn1.open.ac.uk/site/maxima 
and follow the instructions there to download and install Maxima on 
your computer. (You can also find a link to the OU Maxima website 
on the MST 


24 website.) 


(b) Create a folder on your computer in which to store your Maxima work 
for this module. You might like to put this on your computer desktop, 


or in your usual documents area. 


If you are studying within a prison or other closed institution, you can 
install Maxima from the Offline resources disc. 
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2 Getting started with Maxima and 
wxMaxima 


In this section you will learn about the basics of the wxMaxima interface 
and how to use it to perform simple calculations with Maxima. 


If you encounter unexpected problems when working through this section, 
remember to check the Frequently asked questions (FAQ) section on the 
OU Maxima website: learni.open.ac.uk/site/maxima. 


If you are using the command-line interface you should not study this 
section, but study the alternative material in the ‘Maxima accessibility 
guide’ (see page 94) instead. If you are using another interface, you should 
make sure that you can achieve similar outcomes with that interface. 


2.1 Using wxMaxima 


Computer activity 2 Starting wxMaxima 


Start wxMaxima on your computer, and keep it open to work with as you 
study this section. 


If you cannot remember how to do this on your computer, see the 
OU Maxima website: learni.open.ac.uk/site/maxima. 


The wxMaxima interface to Maxima, illustrated in Figure 4, consists 
mainly of a large white working area (called the worksheet) where you 
can type Maxima commands. These commands either instruct Maxima to 
perform calculations or change system settings. 


2. Getting started with Maxima and wxMaxima 


Figure 4 The initial wxMaxima window 


When the software first starts, a secondary window showing the “Tip of the 
Day’ is also displayed. You can close it by clicking the ‘Close’ button. If 
you want to prevent this window from appearing the next time you start 
wxMaxima, click the box in its bottom left-hand corner, so that it is no 
longer ticked, before clicking ‘Close’. (You can always reinstate these 
messages at any time by selecting Show tips from the Help menu.) 


Along the top of the main window (or along the top of the screen in the 
case of an Apple Mac) is a list of menu names, including the File menu 
which allows you to save or print your work and re-open previously saved 
sessions (there is more on this later). You can also close the software from 
the File menu by selecting Exit. (On an Apple Mac you can do this by 
selecting Quit wxMaxima from the wxMaxima menu.) 


Below the row of menu names is a bar called the toolbar, containing a set 
of icons that provide shortcuts to commonly-used functions. 


At the bottom right-hand corner of the wxMaxima window is a message 
that indicates the current status of the Maxima system. When Maxima is 
waiting for instructions from you this reads ‘Ready for user input’. Other 
possible messages include ‘Maxima is calculating’, or ‘Reading Maxima 
output’. 3 


The following activity shows you how to use wxMaxima to perform a 
simple calculation. Follow the steps using wxMaxima as you read the 
activity. 


10 


2 Getting started with Maxima and wxMaxima 


Follow these steps to use Maxima to calculate 2+ 3 x 4 


1, Click anywhere on the large white area of the wxMaxima window, to 
ensure that the text you type next is directed to the correct place. 


2. Type the following exactly as it appears here: 
2+3*4/5; 


including the semicolon (;) at the end. 


Notice that as you type your command, it is shown in 
wxMaxima preceded by a red arrow: ~~>. This is the input prompt 
and shows where commands can be typed. 


To the left of the input prompt is the cell marker 


[ 


which will be described in more detail later. 


3. Press Shift-Enter, that is, hold down the key while pressing 
. (Remember: if, as on some Macs, your keyboard has 


separate and (Return ) keys, then you should hold down the 
([shift_) key while pressing (Return ) .) 


If you make a mistake, try clicking elsewhere on the workspace, below 
your error, and try again. 


Pressing Shift-Enter instructs wxMaxima to send the expression to 
Maxima to calculate. This is known as evaluating the expression. 


The input prompt (-->) is now replaced by (%i1). This is how 
Maxima labels lines. The i stands for ‘input’: this is input line 
number 1. The percentage symbol (%) is used by Maxima to indicate 
objects built into the system. You will learn more about this later. 


The status in the bottom right-hand corner of the window first 
changes to ‘Maxima is calculating’, then ‘Reading Maxima output’, 
before the result of the calculation is displayed, as shown below. 
This may take a little time for the first calculation in any session. 
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The result is displayed under the input line, preceded by the label 
(401). The o stands for ‘output’, so this label identifies the line as 
output line 1. It gives the output corresponding to input line 1. 
Throughout each Maxima session the line numbering starts at 1 and 
then increases by one for each new input (and output). 


(Ril) 24+3"4/57 


22 
sel} — 
(sol) 


In wxMaxima, matching input and output lines form a cell. Each cell is 
identified with the marker shown in Figure 5. 


Figure 5 The wxMaxima cell marker 


The top part of the cell marker (between the upper two horizontal marks) . 
corresponds to the input line, and the bottom part corresponds to the 

output line. The cell marker turns red while an input line is being edited. 

and is shown surrounded by a box while the cell is being evaluated, as 

illustrated in Figure 6. 


Figure 6 The cell marker during evaluation 
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1. Click on the triangle at the top of the cell marker. 


This collapses the cell to a single line. This can be useful if a 
calculation or result is very long and so makes the worksheet difficult 
to read. 


The triangle at the top of the cell marker becomes filled. 


2. Click the filled triangle at the top of the cell marker. 


This expands the cell again. 


3. Click part of the cell marker other than the top triangle to highlight 
it, then press the keyboard key. 


This deletes the cell from the worksheet. 


In Computer activity 3, you used Shift-Enter rather than to 
instruct Maxima to evaluate expressions, and the input prompt (-->) 
appeared only once you started typing. Both these behaviours can be 
disconcerting, but can be changed, as described in Computer activity 5. It 
is recommended that you configure wxMaxima to work in this modified 
way, and in the remainder of this Guide it is assumed that you have done 


so 
1. From the Edit menu, select Configure (or on an Apple Mac 
computer, from the wxMaxima menu, select Preferences). 


Alternatively, click the following icon on the toolbar. 


a 


| This opens the ‘wxMaxima configuration’ window. 


2. Tick the ‘Enter evaluates cells’ box in the Configuration window by 
clicking it. 


a3 
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This changes the behaviour so that pressing alone evaluates 


commands. 


3. Tick the ‘Open a cell when Maxima expects input’ box in the 
Configuration window by clicking it. 


This ensures that the input prompt is shown whenever 
Maxima is waiting for input. 


4. Click ‘OK’ (or close the window if using an Apple Mac). 


i oS == ~ - o* 


Troubleshooting Maxima 


If you are using Maxima and it does not seem to be responding, try the 
following suggestions. 


If you encounter other problems, check the Frequently asked questions on 
the OU Maxima website: learni.open.ac.uk/site/maxima. 


Troubleshooting Maxima 


Cell Marker | Status message Comments 


Boxed: C Maxima is calculating Maxima is evaluating a command that is 
taking a long time to complete. 


Boxed: (| Ready for user input Maxima is waiting for you to type 
something before continuing. 


In either case you might wish to abort the 
command in one of the following ways. 


¢ Click @ on the toolbar. 


e Select Interrupt from the Maxima menu. 
e Type Ctr1-G (on an Apple Mac hold 

down the key while pressing yy 
Maxima might take some time to respond to 
this request. 


Maxima process terminated | Maxima (but not the wxMaxima interface) 
has stopped working. 

Select Restart Maxima from the Maxima 
menu. 
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Closing wxMaxima 


When you have finished using Maxima, you can close Maxima and the 
wxMaxima interface by using one of the following methods. 


e Select Exit from the File menu (or Quit wxMaxima from the 
wxMaxima menu on an Apple Mac computer). 


e Press Ctrl-Q. (Here, Q stands for quit.) 


e Ona Windows computer, click the usual small cross button at the top 
right-hand corner of the wxMaxima window. 


You may be asked if you want to save your changes to the worksheet 
before closing. At this stage, click ‘No’. You will learn how to save your 
work in Subsection 2.7. 


2.2 Basic calculations using Maxima 


In Computer activity 3 you were asked to enter a numerical expression into 
Maxima using the symbols + for addition, * for multiplication and / for 
division. Maxima evaluated it using the rules for the correct order of 
mathematical operations; that is, the BIDMAS rules. The way in which a 
calculation or a command has to be entered in Maxima is known as its 
syntax. 


You should end all Maxima commands with either a semicolon (;) (as you 
saw in Computer activity 3) or a dollar sign ($). Ending a command with a 
dollar sign tells Maxima to perform the instruction, but not to display the 
answer. When using wxMaxima, if you forget to finish your command with 
a semicolon or dollar sign, then the system will add a semicolon for you. 


The Maxima syntax for basic mathematical operations is summarised as 
‘ollows. 


Basic mathematical operations 


Operation Example 
Addition + 4+3; 
Subtraction - 4-3; 
Multiplication * 4x3; 
Division / 4/3; 
Brackets Cand) | 2%(3+4); 
Powers, for example 2° * or *# | 243; 
2**3; 
Square root, for example V5 | sqrt (™) | sqrt (5); 
Magnitude (absolute value) | abs(™) | abs(-3); 


In this Guide, we will use * for powers. Note that many Maxima 
commands take the form of a command name, such as sqrt or abs, 
followed by a pair of round brackets containing one or more objects, such 
as numbers, that the command operates on. Such an object is called an 
argument of the command. In the summary tables in this Guide, 
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arguments of commands are often represented by jf, as illustrated in the 
previous table. 


Warning: Do not omit multiplication signs! 


A common error when first using Maxima is to omit multiplication 
signs. For example, you might type 


e = =6sqrt(2) rather than 2*sqrt (2), or 

e =2(3+4) rather than 2*(3+4). 

Doing this will result in an error. The displayed error message might 
include one of the following statements: 

e parser: incomplete number; missing exponent? 

© = incorrect syntax: Syntax error 


incorrect syntax: Too many )’s 
incorrect syntax: Premature termination of input at ;. 


e ... is not an infix operator 


Tf you receive one of these error messages, first check for missing 
multiplication signs! 


Warning: Take care with question marks (?) 


Do not type question marks when asking Maxima to evaluate an 
expression. The symbol ? is a special symbol in Maxima which 
provides access to underlying systems. 


The one exception to this is when you are using the commands for 
obtaining help from the system, which you will see in Subsection 2.8. 


Maxima generally ignores all spaces that you type while entering a 
calculation or command, so you can use spaces to make a command easier 
to read. You can also enter line breaks, using Shift-Enter, which are 
similarly ignored by Maxima. 


Whenever you enter an opening bracket, the wxMaxima interface 
automatically adds a closing one for you. This behaviour can be turned off 
by unticking Match parenthesis in text controls in the Configuration 
window. (See Computer activity 5 for how to open the Configuration 
window.) 


You can calculate roots other.than square roots, such as cube roots, by 
using the index laws. 
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(a) Remembering that %/27 = 273, calculate {27 by entering the 
following line into Maxima. 


27°(1/3) ; 


Remember to press (Enter) to evaluate the command. 
(b) What happens if you enter 
OT 1/33 


instead? Can you explain this? 


Maxima usually tries to calculate quantities exactly rather than 
approximating them by decimal numbers. It does this by manipulating the 
symbols in the expression according to the rules of mathematics. This is 
sometimes known as symbolic computation. You can see an example of 
this in the next activity, which also shows you two different ways to 
instruct Maxima to output a decimal answer instead. 


In Maxima, calculate 1309/3451 in three different ways, by entering each 
of the following expressions in turn. 


(a) 1309/3451;  (b) 1309.0/3451; — (c) float(1309/3451); 


—— 
Maxima returns an exact answer for part (a). 


Since part (b) involves a decimal number, 1309.0, Maxima returns a 
decimal approximation to the answer. 


In part (c), the command float instructs Maxima to return a 
decimal answer. 


Converting to decimal numbers 


Operation Command | Example 


Convert to a decimal number | float ({§) | float (sqrt (2)); 


The name of the command float arises from the fact that the method 
used by computers to store decimal numbers internally using the binary 
digits 0 and 1 is called a floating-point representation. The float 
command instructs Maxima to convert a number stored symbolically as a 
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mathematical expression to one stored using a floating-point 
representation. 


Maxima usually displays decimal numbers to 16 significant figures (though 
it suppresses trailing zeros at the end of the decimal part of a number; for 
example 0.125 is displayed simply as 0.125 rather than as 
0.1250000000000000). This is also the precision to which Maxima 
performs decimal calculations. You will see in Subsection 2.5 of this Guide 
how to change the number of significant figures displayed. 


For each of the following expressions, use Maxima to simplify it and then 
find a decimal approximation for it. 
8 


2/3 
(a) /i47 —(b) (3) (c) 2500 


The expression in part (c) doesn’t involve decimal numbers and so the 
result is evaluated exactly. Notice, however, that the middle 91 digits 
of the answer are not displayed because the system does not expect 
you to read them all! 


The decimal form of the result, 3.2733906078961419 10*®°, is given 
in wxMaxima’s way of displaying scientific notation. It means 
3.2733906078961419 x 101°. 


Maxima knows the values of standard mathematical constants such as the 
irrational numbers e and 7. It treats them symbolically unless you 
explicitly ask for them in floating point form. 


The constants e and 7 


As with the input and output line labels, the “%’ symbol indicates the 
name of a quantity built into Maxima. 


When displaying mathematics in an output line, wxMaxima displays %pi 
as T. 5 
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Use Maxima to calculate the value of the expression 7 + : +e, both 
exactly and as a decimal number. 


Notice that in the exact answer, %pi is displayed as 7, but %e is 
unchanged. 


2.3 Reusing and editing commands 


In Computer activity 9, you may have entered the expression twice, once 
to evaluate it exactly and once for a decimal answer. This is, however, not 
necessary. When you have entered and evaluated an expression, you can 
edit it and then evaluate it again, as shown in the following activity. 


(a) Calculate 2V3 + 527 by entering the following line into Maxima. 
2x*sqrt (3)+5*sqrt (27) ; 


(b) Edit your entry in part (a) to calculate 2/3 — 5/27 as follows. 


1. Click on the expression you entered in part (a), or move to it 
using the up and down keyboard arrow keys. 


As you do this the cell marker will turn red and you will see a small 
flashing vertical line (the editing cursor) appear in the expression. 


(%i1) 2* sqrt (3) +45*sqrt (27); 
($01) 1713 


2. Edit the expression by typing new text and using the 
(backspace) or (Delete ) key to make deletions. 


3. Press : 


Pressing re-evaluates the cell. When a cell is re-evaluated in 
this way, the input and output line numbers are updated as if a new 
cell had been added. 
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Should you click in the space between two cells, a horizontal cursor will be 
shown, as in Figure 7. This indicates the position at which a new cell will 
be created, if you begin typing here. This is also shown as you move 
between cells using the up and down keyboard arrow keys. 


($i1) Spi + %pi/3 + %e*2; 


4 
(%01) tte? 


(%i2) float (%pi + tpi/3 + %e*2); 
($02) 11.57784630371704 


Figure 7 The wxMaxima horizontal cursor 


In wxMaxima, you can also copy and paste expressions, or parts of 
expressions, as you might copy and paste when using other software; that 
is, using Copy and Paste from the Edit menu, or the Ctrl-C and Ctr1-V 
keyboard shortcuts. 


Maxima also permits the result of one calculation to be used in another 
calculation. The symbol % represents the result of the last cell evaluated. 
Due to the updating of cell numbers noted above, this might not be the 
last cell in the worksheet, but it will be the cell with the highest line 
numbers. So, for example, if you enter sqrt (%) ; Maxima calculates the 
square root of the last evaluated output. You can also use the output (or 
input) of other cells, by typing the cell line label in a calculation. For 
example, if you enter 405%3; Maxima calculates the cube of the expression 
in output line 5. 


(a) Enter 6+2; into a cell at the bottom of your worksheet to evaluate it. 
(b) On the following line type 3+%; and enter it. 


Maxima returns 11, the sum of 3 and the previous answer. 


(c) Edit the line where you typed 6+2;, change it to 6*2; and re-evaluate 


the cell by pressing B 


Notice that the cell number and output of the line you edited is 
updated, but the result of 3+%; does not change. This is because 
when the 3+%; cell was evaluated, the symbol % represented the result 
of the old version of the edited command. 
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(d) From the Cell menu, select the Evaluate All Visible Cells 
option, or press Ctr1-R. 


This re-evaluates all visible cells in the worksheet, in the order in 
which they appear, updating all the line numbers as it does. 


Notice that the result of 3+%; has now changed. At the time it was 
re-evaluated the symbol % represented the value of the last evaluation, 
which was the revised version of the previous line. 


(e) In a new cell at the bottom of your worksheet enter Zoff+5; 
with §§ replaced by the output line number corresponding to the result 
of 6*2 entered earlier. 


This command adds 5 to the result of the line referenced. 


2.4 Assigning values to variables 


Maxima allows you to assign a value to a variable, and then use the 
variable in further calculations. 


For example, you can assign the value 23 to the variable a, and then 
evaluate an expression in a, such as a? — 3a +2. You can also assign 
expressions, both numerical and (as you will see later) algebraic, and 
indeed any other Maxima object, to a variable, and use that variable in 
subsequent commands. 


To assign a value to a variable you use a colon (:). For example, the 
command 


a:23; 


assigns the value 23 to the variable a. The colon can be read as ‘is 
assigned the value’. 
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Variable names can be any combination of letters and numbers that begins 
with a letter. For example, Maxima will accept any of the following 
variable names. 


a A solution solution2 x2b 


Variable names are case-sensitive; for example, the variable x is different 
from the variable X. If a variable name is the name of a Greek letter, then 
wxMaxima will display it in output lines as the proper Greek character; for 
example, the variable alpha will be displayed as a. 


Once a variable has been assigned a value, Maxima will remember this 
value for the rest of your session. This may sometimes surprise you. It is 
easy to forget that earlier in the day you assigned a value to, say, x; you 
may then be confused later on when a calculation involving x gives an 
unexpected result. You can, however, tell Maxima to ‘forget’ about a 
variable to which you have previously assigned a value. You can also ask it 
to list all the variables that you have assigned values to. These commands 
are demonstrated in the following activity. 


(a) Assign the value of V8 to the variable a by entering 
a: sqrt(8); 


The value of a is displayed as output, in an equivalent form. 


(b) Enter a; 


This displays the current value of a. 


(c) Assign the value of a2 to b. 


Don’t forget to include a multiplication sign when typing aV2. 
The value of b is simplified when displayed. 


(d) Edit the line where you assigned V8 to a, so that the value V7 
becomes assigned to a. 
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(e) Enter b; to display the value of b. 


Notice that the value of b has not changed. The variable b was 
defined when a had its original value. 


(f) Enter values; to list the names of all the variables that are currently 
assigned values. 


The variable names are given within square brackets, separated by 
commas. This is how Maxima displays lists. 


(g) Enter kill(a); to remove the variable a from the system. 


You should obtain the output done. 


(h) Enter a; to display the value of a. 


The variable name is output, as a no longer has a value. 


(i) Enter b; to display the value of b. 


b is still defined. | 


eae ee eee 
Notice that in Computer activity 12(d), when the value of a was changed, 
the value of b was not affected, since it was defined when a had its original 
value. If you want to update the value of b using a different value of a, 
then you should enter a new value of a then re-evaluate the line defining b 
(by selecting the line and pressing (Enter_)), or re-evaluate the whole 
worksheet by selecting the option Evaluate All Visible Cells from the 
Cell menu. 


Working with variables 


Operation Command | Example 
Assign a value to a variable : a:23; 
Display the value of a variable | Wariable) a; 

List all user assigned variables | values values; 
Remove an assigned variable | kill ((WaeHabI@)) | kill(a); 
Remove all assigned variables | kil1(al1) kill(all); 


Note: here the placeholder (¥a#@abI) represents any variable name. 
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You can practise using variables in the following activity. 


(a) Define the variable a to have the value 42. 
(b) Define the variable b to be equal to a® — 3a + 4. 


(c) Define the variable c to be the decimal approximation of the square 
root of b. 


2.5 System variables 


There are some variables built into Maxima whose values affect the 
behaviour of the system. These are called system variables. You met 
one system variable in Computer activity 12; the variable values holds a 
list of the names of all the variables you have defined. 


Another system variable, fpprintprec, specifies the number of significant 
figures of decimal numbers that are displayed. The name fpprintprec is 
an abbreviation of ‘floating point print precision’. 


To change the system behaviour so that, for example, only 4 significant 
figures are displayed, use the command fpprintprec:4; to assign the 
value 4 to the variable fpprintprec. 


You can set the variable fpprintprec to any value between 2 and 16. 
Also, setting it to 0 restores the default behaviour of displaying 16 
significant figures. The value of fpprintprec sets the requested number of 
significant figures, but Maxima does not always exactly meet the request. 


Display the decimal approximation of 7 to 8 significant figures by doing 
the following. 


(a) Set the system variable fpprintprec to be 8. 
(b) Display the decimal approximation of z. 


You can reset the values of fpprintprec and all other system variables to 
their original values by using the reset() command. Try entering 
reset(); followed by float (%pi) ; and check that the value of 7 is 
displayed to 16 significant figures. 
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Resetting system variables 


| Operation Command | Example | 


| Reset all system variables | reset () | reset(); 


Note that the list of variables displayed when you enter values; does not 
include any system variables. If one of your variables does not appear in 
the list displayed by values;, then you have probably used a variable 
name that is also the name of a system variable and hence changed the 
value of that system variable. In extreme circumstances this may change 
the behaviour of Maxima, which can be restored by resetting all system 
variables as described above. 


2.6 Annotating your wxMaxima worksheet 


It is important to be able to clearly explain and present your mathematical 
work, both so that someone else can understand it, and so you can 
understand what you have done if you return to your work after some 
time. To help you do this, wxMaxima allows you to enter text comments 
within your worksheet, as illustrated in Figure 8. 


(%il) 4:37 
(tol) 3 


(#12) ¢: float(2*tpi*r): 
(#02) 18.84955592153876 


(814) a: float (Spite*2); 
($04) 28.27433388230814 


Figure 8 Adding text to a wxMaxima worksheet 
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To enter text into your worksheet, first you need to create a suitable cell in 
which to type. To do this first you position the cursor, then you select the 
appropriate one of the following commands from the Cell menu. 


e Insert Text cell 

e Insert Title cell 

e Insert Section cell 

e Insert Subsection cell 

The new cell will be created under the cell in which the editing cursor 
appears, or between two existing cells at the position indicated by the 
horizontal cursor, if shown. Remember, the horizontal cursor appears 
whenever you click between two cells, or move within the worksheet using 
the up and down keyboard arrow keys. 

As you type your text, to start a new line within the cell press 

When you have finished typing the text, click elsewhere on the works Bee 
or use the up or down keyboard arrow keys to move out of the text cell. If 
your text cell is at the bottom of the worksheet, and you click or move 
below it, an input prompt may not be shown until you start typing again. 


Section and Subsection cells are automatically numbered. 


(a) Create the worksheet shown in Figure 8. Use a Title cell for the title 
‘Circles’ and Text cells for the remaining text. Remember to click 
elsewhere on the worksheet or move position with the up or down 
keyboard arrow keys when you have finished entering each piece of 
text. 


After evaluating a Maxima command, you will need to correctly 
position the horizontal cursor before inserting the next text cell. 


(b) Use your worksheet to find the circumference and area of a circle of 
radius 5. (Remember to re-evaluate the worksheet, for example, by 
using Ctr1-R, after changing the value of r.) 


Each Title cell, Section cell and Subsection cell has a small square in its 
top left corner. Clicking this square hides (or reveals) the contents of the 
worksheet from that point onward, until the next cell of the same type. 
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2.7 Saving and printing your work 


After working with Maxima you will probably want to save your 
calculations so that you can reuse or read them at a later date. You can do 
this in wxMaxima using the Save or Save As options from the File menu. 


There are three different formats in which you can save your work. These 

are as follows. 

e wxMaxima document, with the .wxm file type. This saves all your 
input but not any calculated output. You have to recalculate the 
output when you re-open the file. 


e wxMaxima xml document, with the .wxmx file type. This saves all 
the input and the output. 


e Maxima batch file, with the .mac file type. This is for more 
experienced users of Maxima. 


You can open a saved wxMaxima worksheet by using the Open option from 
the File menu, or by double-clicking on the saved file on your computer. 


Save your current worksheet then reload it, by following these steps. 
(a) Select Save As from the File menu. 
(b) In the window that appears: 


(i) Choose the folder in which to save your work. You might like to 
save it in the folder that you created in Computer activity 1. 
(ii) Choose a name for the file in which the worksheet will be saved. 
(iii) Make sure the file type is set to wrMazrima xml document 
(*.wamx). 
(iv) Click ‘Save’. 
(c) Close wxMaxima, by selecting Exit from the File menu, or typing 


Ctr1-Q, or clicking the small ‘x’ at the top right of the window. (On a 
Mac select Quit wxMaxima from the wxMaxima menu or hold down the 


key while pressing (Q).) 
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(d) Open your worksheet in wxMaxima again, either by 
(i) double clicking on the file you just saved, or 
(ii) following the steps below: 
e Start wxMaxima again. 


e Select Open from the File menu. Maxima will ask if you 
want to save your just opened, blank worksheet. Click ‘No’. 


e Find the file containing your work and click ‘Open’. 
CSS 2 eee 


You can print a wxMaxima worksheet by selecting Print from the File 
menu. 


Alternatively, you can obtain an image of the entire wxMaxima window by 
taking a ‘screenshot’. To do this on a Microsoft Windows computer, first 
click on the wxMaxima window, then press Alt-PrintScreen (that is, 
hold down the keyboard key while pressing the PrintScreen key, 
which may be labelled “ or something similar). This stores 
the image, which you can then paste into a suitable application. 


If you are using Windows Vista, Windows 7 or Windows 8, you might like 
to use the ‘Snipping Tool’, available by typing ‘snip’ into the Start menu. 


To take a screenshot on an Apple Mac computer, press Cmd-Ctr1-Shift-4 
(that is, hold down the {command} eon and (@) keyboard keys while pressing 
the key). Next press the space bar, move aK camera pointer over the 

‘axima window and click. This stores the image which you can then 
paste into a suitable application. You may also like to use the Grab 
application from the Utilities folder. 


If you want to copy a single output line from your worksheet into, for 
example, a word-processed document, then you can select that part of the 
worksheet, and choose Copy As Image from the menu obtained by right 
clicking on what you have selected. You can then paste the resulting image 
into a suitable application. 
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2.8 Getting help 


There are several ways in which you can obtain help with using Maxima, 
beyond this Guide. 


In wxMaxima, you can access the complete Maxima manual by selecting 
Maxima Help from the Help menu. This opens the window shown in 
Figure 9. You can use the Contents, Index and Search tabs to look for the 
information that you want. 


| Maxima 5.30.0 Manual 
Maxima is a computer algebra system, implemented in Lisp_ 


‘Maxima is derived from the Macsyma system, developed at 
MIT in the years 1968 through 1982 as part of Project MAC 
MIT tumed over a copy of the Macsyma source code to the 
‘Department of Energy in 1982; that version is now known as 
DOE Macsyma A copy of DOE Macsyma was maintained by 
Professor William F. Schelter of the University of Texas from 
1982 until his death in 2001. In 1998, Schelter obtained 

|| permission from the Department of Energy to release the DOE 
‘Macsyma source code under the GNU Public License, and in 
2000 he initiated the Maxima project at SourceForge to 


Maxima infrastructure 


Sy eey ‘Sample Maxima sessions. 


fee 


(=e 


Figure 9 The Maxima help system 


Maxitna also includes commands that you can use to obtain help. These 
are demonstrated in the following activity and summarised in the table 
that follows. 


(a) Enter ? float; to show the help for the float command. 


Note the space between the ? and float, which is needed. This is one 
of the few times in Maxima when spaces matter! 


The help information for the command may include more detail than 
you need! 


After giving the help information, Maxima displays the output true. 
This means that the command you entered was successful. 
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(b) Suppose that you cannot remember the Maxima command for square 
root, but you do remember that it was something like sqr. 


Type ?? sqr; to list all the help information titles that contain the 
letter sequence ‘sqr’. 


A list of possible titles is displayed. 


Title number 1 is what you were looking for, so enter 1; to display the 
relevant information. 


If more than one title is displayed in response to a 7? command, then 
Maxima will not continue until you enter one of the following: 


e the appropriate number, to display your chosen information 


e several numbers, separated by spaces, to see help information 
on several topics 


e all to display the information on all topics listed, or 
e none to see no information. 


A box will be shown around the cell marker while Maxima is waiting 
for your response. 


Help commands 


Operation | Command Example 
Get help on a command ? (command ? float; 
or, describe ((€ommand ) describe (float) ; 


Find help information whose | ?? ?? flo; 
title contains the given text | or, describe(J, inexact) | describe(flo, inexact); 


Note 1: the space after the ? is needed. 


Note 2: the second pair of commands displays a list of the titles of all 
help information pages whose title contains the given text. To obtain 
help on a particular topic listed, type the number corresponding to the 
required title in the list, followed by a semi-colon. 


There is further help with Maxima in the Frequently asked questions 
section of the OU Maxima website, at learn1.open.ac.uk/site/maxima. 


3 Algebra 


In Section 2 of this Guide you used Maxima to perform numerical 
calculations. Much of the power of Maxima, however, lies in its ability to 
perform algebraic operations. In this section, you will learn how to do this. 


3 Algebra 
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3.1 Algebraic manipulations 


You can enter algebraic expressions into Maxima in a similar way to 
numerical expressions, using the symbols +, -, *, / and *, and the 
command sqrt (J) for square roots. You can also assign algebraic 
expressions to variables, using the assign (:) command, in the same way-as 
numerical expressions. 


In the next activity, you will meet some commands for manipulating 
algebraic expressions. 


(a) Enter 
expand( (2+x)*(3*x+4)%3 ); 
to expand the expression (2 + «)(3a + 4)°. 


The command expand multiplies out the expression that is its 
argument. 


Spaces have been included simply to help with the reading of the 
command. 


(b) Enter 
factor( 2*x”*3-7*x"2-10*x+24 ); 


to factorise the expression 22° — 7a? — 10x + 24. 


The command factor factorises, where possible, the expression that 


is its argument. 


x41 
a+1- 


(c) Define p to be the expression 


Remember to use brackets as necessary when typing this expression. 


(d) Enter 
fullratsimp(p) ; 


to simplify the expression assigned to the variable p. 


3 Algebra 
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3 Algebra 


The command fullratsimp simplifies, where possible, the expression 
that is its argument. It does this by using techniques for multiplying 
out brackets, collecting terms over a common denominator and 
cancelling out common factors. Simplifications requiring other 
techniques may not be handled by fullratsimp in the way you would 
expect. 


The name fullratsimp is short for ‘full rational simplify’. 


(e) Enter 
subst (3,x,p); 


to find the value of the expression p (defined above) when x = 3. 


The command name subst is short for ‘substitute’. 


The command subst (3,x,p); substitutes the value 3 for the variable 
x in the expression p. 


The commands introduced in the previous activity are summarised below. 


Algebraic manipulation commands 


SC i 


Expand brackets | expand (J) expand( (x+1)%2 ); 
Factorise factor(™) factor( 2*x+4*x"2 ); 
Simplify fullratsimp(() fullratsimp( (2*x+4*x%2)/x ); 
Substitute subst (Value), WarHabIe), | subst (4, x, x*2+1); 
‘expression ) which substitutes 4 for x in x42+1 


It is hard for a computer to simplify algebraic expressions, because often 
there is no single ‘simplest form’ for a particular expression, and the 
computer has no way of knowing which form you prefer. You can use the 
commands in the table above to simplify an expression in the way you 
want. 


You can access many of these commands from wxMaxima menus. For 
example, to factorise an expression you can first enter the expression in the 
worksheet, select it by highlighting the text using the mouse or clicking on 
the corresponding output line number, then select Factor Expression 
from the Simplify menu (or from the menu obtained by right clicking on 
a highlighted output line number). 


Other options on these menus, such as Substitute... open a window for 
you to enter the command options. Options that work in this way are 
identified by ... at the end of the option name. 
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3.2 Equations and their solutions 


Maxima can manipulate and solve equations. To type an equation in 
Maxima, you use an equals sign in the usual way. For example, you can 
type 3#x+2=4, 


You can assign equations to variables in a similar way to expressions. 
For example, by entering 


q: 3*x+2=4; 
you assign the equation 3x + 2 = 4 to the variable q. 


You can extract the left- and right-hand sides of an equation by using the 
lhs and rhs commands, as demonstrated in the following activity. 


(a) Assign the equation x? + 2x = 4 to the variable r. 
(b) Extract the left-hand side of the equation by entering 
lhs(r); 


(c) Extract the right-hand side of the equation by entering 
rhs(r); 


Extracting the left- and right-hand sides of an equation 


Left-hand side of an equation | lhs ((@qwati6n)) | lhs(4*x+1= 
Right-hand side of an equation | rhs ((@quatien) | rhs (4*x+1=2*x 


Note: (@q@@ation) indicates that the argument of the command must be 
an equation, or a variable whose value is an equation. 


Maxima can find the exact solution(s) of equations, as demonstrated in the 
following activity. 
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a) Solve the equation x? + 2x = 1, by entering 
ju 


solve( x°2+2*x=1 ); 


Notice that Maxima displays the two solutions of the quadratic 
equation within square brackets and separated with a comma: 


[x=-V2-1, x=/2-1] 


This is a Maxima list. 


(b) Use Maxima to solve 4a? — 127 +9 = 0. 


Notice that although there is only one solution of this quadratic 
equation, Maxima still gives it as a list, containing just one element. 


c) Use Maxima to solve 2x + x? —5a+2=0. 


———— 
This is a cubic equation, since the highest power of x is 3. Each 
cubic equation has at most three solutions. 


(d) Try to use Maxima to solve x° + 2x —4 = 0. 


Here, Maxima returns (within a list) the equation, slightly 
rearranged, but unsolved. 
The solve command always tries to find an exact solution of the 


equation using algebraic manipulations. If it is unable to do this, then 
it returns the original equation, as here. 


In Section 6 of this Guide you will see how to find approximate 
solutions of such equations. 


(e) Use Maxima to solve x? — 4a +5 =0. 


————<—— 
This equation has no real solutions, since the discriminant, 
(-4)? -4x 1x 5 =—4, is negative. 


However, Maxima finds solutions that are complex numbers. These 
involve the imaginary number, usually denoted by i, whose square 
is —1. Maxima denotes this number by %i. 


Complex numbers were mentioned in Unit 1 and you will learn more 
about them in Unit 12. 


3 Algebra 
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The solve command is summarised below. 


Solving an equation 


| Operation Command Example 


| Solve an equation, exactly | solve ((@quatien)) | solve(2*x*2-1=0) ; 


You saw in Computer activity 20 that the output of the solve command is 
a list. In Maxima, a list is a collection of individual elements separated 
by commas and enclosed within square brackets. You can learn more 
about using lists in the following activity. 


(a) Assign the list containing the first five prime numbers to the variable A 
by entering 


A; 1253,5,7523); 
(b) Display the fourth element of the list by entering 
A(4]; 


The position number of an element in a list is known as its index. 
You can obtain individual elements of a list by typing the index of the 


element in square brackets after the list, or after the name of a 
variable to which the list is assigned. 


(c) (i) Assign the solutions of the equation x? — x — 1 = 0 to the 
variable solns. 


(ii) Display decimal approximations to the answers by entering the 
command 


float (solns) ; 
(iii) Display the exact value of the first solution by entering the 
command 


rhs(solns[1]); 


The command solns[1] refers to the first element of the list solns 
Tee ont 
which is Soe: 


This is an equation. The value needed is the right-hand side of this, 
which is extracted using the rhs command. 


3. Algebra 
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3 Algebra 


The operations introduced in the above activity are summarised below. 


Lists 
Operation Command 
Form a list CQ. @. --- 1] A: [a,b,c]; 
Extract an element of a list | (ist) [Gindex] 
Note: ... indicates the list of arguments could continue. 


As well as solving equations, you can also manipulate equations using 
Maxima. The next activity shows how you can do this. 


(a) Assign the equation 3y + 4 = 3x + 92 to the variable eq. 
(b) Enter the following command. 


eq-4; 


This subtracts 4 from both sides of the equation eq. 


(c) Now divide both sides of the rearranged equation by 3, by entering 
h/3; 


Remember that % represents the result of the last evaluated cell. 
The variable y is now alone on the left-hand side of the equation. 


We have rearranged the equation to make y the subject. 


You can also use the solve command to change the subject of an 
equation: you solve the equation for the variable that you want to be the 
subject, as shown in the next activity. 
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2a+b+3c 
a—b 
(b) Rearrange the equation to make a the subject by entering 


(a) Assign the equation ¢ = to the variable eqn. 


solve(eqn, a); 


Here, the second argument of the solve command is the variable that 


you want to solve for; that is, the variable that you want to be the 
subject of the equation. 


This use of the solve command is summarised below. 
Changing the subject of an equation 


Operation Command Example 


Solve an equation for a variable | solve (@@uatien, WaFGabIe)) solve(2*a*b-3*b=0, a); 


Another use of the solve command is to solve two or more simultaneous 
equations. To use the command in this way, you must enter the equations 
as a list, and enter the variables to be solved for as another list, as 
indicated below. 


Solving simultaneous equations 


Operation Command Example 


3 Algebra 


Solve simultaneous equations 


solve (QistJof equations), | solve([2*x+y=4, 3+#x-2*y=-1], 
list of variables ) (x,y]); 


| 


Use Maxima to solve the simultaneous equations 
3a? + 2y = 20 
4x — 3y = —4 
as described below. 


(a) Assign the first equation to the variable eqi and the second equation 
to the variable eq2. 


Maxima might reorder the terms of the equations when displaying the 
output. 
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(b) Solve the pair of equations (eqi and eq2) for the two unknowns 
(x and y) by entering 


solve( [eqi, eq2], [x,y] ); 


Maxima’s output is [[x=2, y=4] , [ x= 26 a y= $8 Jd 


This tells you that there are two solutions to the equations: one 
solution is = 2, y =4 and the other is x = —26/9, y = —68/27. 
Notice that Maxima’s output here is a list of lists. 


3.3 Plotting graphs 


You can plot graphs in wxMaxima by using the wxplot2d command. This 
command inserts your graph into your worksheet, as demonstrated in the 
activity below. Maxima plots the graph by calculating a large number of 
points on the curve and joining them up. The ‘2d’ in the name of this 
command comes from the fact that its graphs are two-dimensional plots. 
The ‘wx’ is because it is a wxMaxima command. 


If you are using another interface to Maxima, you should use the command 
plot2d instead. You use it in the same way as wxplot2d, but it will plot 
the graph in a different computer window, which you then need to close 
before you can create another plot. 


Plot the graph of the equation y = 4x? — 8x + 2 for values of x such that 
—1 <2 < 2, as described below. 
(a) Enter the command 

wxplot2d( 4*x%2-8*x+2, [x,-1,2]); 
=D———$—<—— 
The first argument of the command is the expression to be plotted, 
and the second argument is a list that specifies the variable in the 
expression together with the least and greatest values that you want 
it to take in the graph. 
wxMaxima displays the graph next to a line label beginning 7%t, such 
as (%t1). This is how wxMaxima shows intermediate output, that is, 
output other than the final result of the command. 
The wxplot2d command has no final result, just intermediate output, 
which is the graph. 


3 Algebra 
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(b) Change the vertical range of the graph to be —5 < y < 5 by including 
an additional argument [y,-5,5]. 


That is, change the command to 
wxplot2d( 4*x%2-8*x+2, [x,-1,2], [y,-5,5]); 


Notice the warning: 
plot2d: some values were clipped. 


This means that some of the points on the graph with z-coordinates 
between —1 and 2 could not be displayed due to the range of values of 
y specified. 


(c) Change the colour of the curve to green, by including the additional 
argument [color, green] after the other arguments of the wxplot2d 
command, but before the closing round bracket. 


This argument is a list whose first element is a keyword indicating 
the property to be set, in this case color (with the American 


spelling!). The second element is the value of the property that you 
want to choose. 


You can plot multiple graphs in the same diagram by including a list of 
expressions as the first argument of the command, as shown in the 
following activity. 


Suppose that the displacements s; and s (in km) of two cars along a 
straight road from a given starting point are given by s; = 15t® + 10¢ and 
8 = 60t respectively, where ¢ is the time (in hours) and 0 <t < 2. 


(a) Plot graphs of the displacement of the cars for 0 < t < 2 by entering 
wxplot2d( [15+#t*3+10+*t, 60*t], [t,0,2] ); 


Notice that the first argument of the command wxplot2d is a list 
containing the two expressions to be plotted. In this case, the 
expressions involve the variable t, so the second argument specifies 
the range of values of t to be plotted. 


The curves are automatically plotted in different colours, and a legend 
showing which expression corresponds to which curve is included at 
the top right-hand corner of the graph. 


3 Algebra 
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3 Algebra 


(b) Change the colours of the curves by adding the following additional 
argument to the command used above: [color, green, magenta]. 


That is, change the command to 


wxplot2d( [15+*t*3+10+*t, 60*t], [t,0,2], [color, green, magenta] ); 


Maxima uses the first colour listed for the graph of the first 
expression listed and so on. If you do not include enough colours, 
Maxima starts using the listed colours again. 


(c) Now change the labelling within the legend to be the variable names 
given to the expressions in the question, rather than the expressions 
themselves. Do this by including another argument 


[legend, "si" ,"s2"] 


in a similar way to the colour argument; that is, within the round 
brackets of the wxplot2d command and separated from the other 
arguments by a comma. 


In general, to change the legend text, you include an extra argument 


of wxplot2d: a list that contains the keyword legend, followed by the 
new pieces of text, each enclosed in double quote marks. 


(d) Finally, change the labelling of the horizontal and vertical axes, by 
adding two more arguments to the command: 


[xlabel, “time, t (h)"] 
and 


[ylabel, "displacement, s (km)"]. 


In general, to change the z-axis label, you include a further argument 
of wxplot2d: a list that contains the keyword xlabel, followed by the 


new label text, enclosed in double quote marks. You can change the 
y-axis label in a similar way, using the keyword ylabel. 


The different ways in which you can use the wxplot2d command are 
summarised as follows, along with some of the optional arguments. If you 
are using a different interface to Maxima, you can use the plot2d 
command in the same ways. 
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4 Functions 


Plotting graphs 


Operation Command Example 


Plot a graph wxplot2d (jexpression , wxplot2d(x*2, [x,0,1]); 
lHorizontalrange,...) | which plots the graph of x? 


for x between 0 and 1. 


wxplot2d (x2, [x,0,1],[y,0,2]); 
which plots the graph of x? 

for « between 0 and 1, and with 
vertical axis values between 0 and 2. 


Plot several graphs | wxplot2d (QHSt\OE expressions), | wxplot2d([x*2, 2*x], [x,0,1]); 
(Horizontal ‘range’, ...) which plots graphs of x? and 2x 


for « between 0 and 1. 


Graph plotting options 


Option Argument 
Vertical range ly, a. 


Example 


ly, 0, 5] 


Curve colour [color, 0 [color, red] 
Legend text [legend,"/",... [legend, "Car A"] 

(a label should be listed for each curve plotted) 
Turn legend off [legend, false] [legend, false] 


Horizontal axis label | [xlabel, ")j"] 
Vertical axis label [ylabel, "§§"] 


[xlabel, "time, t (h)"] 
[ylabel, "displacement, s (km)"] 


Note: for a full list of the possible options, view the Maxima help for 
plot2d, using the ? plot2d; command. 


If you began studying Sections 2 and 3 of this Guide from Activity 31 of 
Unit 2, this now completes your study of the unit. 


4 Functions 


You can define a function in Maxima using the := operator, as shown in 
the following activity. 
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(a) Define the function f in Maxima, with the rule f(x) = «° — 3x? + 4, 
by entering the following command. (Use the characters : followed 
by = to type the := operator.) 


f(x) :=x%5-3#x%2+4; 


(b) Calculate the value of f(2). by entering £(2) ; 
(c) Now define the function g with rule 


and find the value of (3). 


Remember to use brackets appropriately in the rule for g. 


(d) Now try to use Maxima to find g(2). 


This gives an error. 


This is because when Maxima tries to evaluate g(2), it finds that this 
involves dividing by zero, so it returns an error message. The function 
g is undefined at x2 = 2. 


(e) Enter h(1); 


Maxima returns the expression h(1), unevaluated. 


Since the function h has not been defined, Maxima cannot calculate a 
value. It knows only that the value is h(1). 


(f) Enter functions; 


This displays the value of the system variable functions, which is a 
list of all the functions that you have defined in this Maxima session. 


Notice the difference between defining a function, which has an argument 
and which can be evaluated at different values of the argument, using :=, 
and assigning an expression to a variable (which has no argument) using :. 
As with variables, once you have defined a function in Maxima it is 
remembered by the system for the rest of your session, unless you remove 
it using the kill command that you met in Subsection 2.4 of this Guide. 


4 Functions 
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Working with functions 


Operation 


Define a function 
Evaluate a function at a value 
List all user defined functions 
Remove a defined function 


= £(x) :=2*x+3; 
function (9) = £(1); 


functions functions; 


kill ((ietiom) | kill(£); 


You can plot the graph of a function using the wxplot2d command 
introduced in Subsection 3.3 of this Guide, as demonstrated in the next 
activity. 


(a) Plot the function f defined in Computer activity 27 over the range 
—2 <2 <2, by entering 


wxplot2d (f(x), [x,-2,2]); 


Note that the first argument must be f(r) rather than just f. 


(b) Plot the function g defined in Computer activity 27 over the range 
=5.5 mS 5. 


Notice the vertical scale of the graph, and the large spike near x = 2. 


The function g has the rule g(x) = — 3? 8° it is not defined at 


ax = 2. Maxima has evaluated it at values of x close to 2, obtained 
large values and attempted to join up the corresponding points. 


To obtain a better graph, you should restrict the range of the vertical 
axis. You are asked to do this next. 


(c) Plot g again, for values of « in the range —5 < a <5, but restricting 
the range of the vertical axis to between —5 and 5. 


Remember, to restrict the vertical range to —5 < y < 5, include the 
argument [y,-5,5] in the wxplot2d command. 


If you began studying this section from Unit 3, Activity 11, now return to 
the unit to continue your study of functions. 
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5 Exponential functions and logarithms 


5 Exponential functions and logarithms 


You saw in Subsection 2.2 of this Guide that the mathematical constant e 
is represented in Maxima by %e. The % symbol indicates that this is the 
name of a special quantity built into Maxima. 


Computer activity 29 Working with exponentials and logarithms 
(a) Find e° using Maxima, by entering 


ne; 


Since Maxima works symbolically, it does not return a decimal 
approximation. 


(b) Find a decimal approximation to e° by entering 
float (4e*5) ; 
(c) An alternative way to find powers of e is to use the exp function. 
Find e° by entering 
exp(5); 


Again, this is expressed exactly, as a power of e. 


(d) The Maxima command for a natural logarithm is log(™). 
Use this to find a decimal approximation for In 8. 


(e) Use the solve command introduced in Subsection 3.2 of this Guide to 
solve the equation 


os, 
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5 Exponential functions and logarithms 


You may be surprised to see two solutions to this equation. 


One solution, Jog(5) is the sort of solution that you would expect. 
2log(3) 
log(—V5) 


The other solution, Toa(3) 


, involves the logarithm of a negative 
number. 

You have seen that only positive numbers have logarithms, but that 
applies when you are working only with real numbers. If you are 


working with complex numbers, then negative numbers do have 
logarithms. 


log(—v5) 


log(3) 
clearly by obtaining a decimal approximation for it. It contains the 
imaginary number i, which Maxima represents as i. 


The solution is a complex number, as you can see more 


There is only one real solution to the equation. 


The Maxima commands for the exponential function and logarithms are 
summarised below. Maxima has no command 1n(/§); remember that the 
command for a natural logarithm is log (§). 


Exponentials and logarithms 


Operation Example 

Exponential, for example e? | %e he*3; 
or exp(§) | exp(3); 

Natural logarithm, In log() log (8) ; 


When you are working with logarithms in Maxima, the algebraic 
simplification commands that you met in Section 3.1 of this Guide can be 
useful. You can also use the commands shown below. 


Simplifying expressions involving exponentials and logarithms 


Operation Command Example 

Simplify something radcan())) radcan( log(x”2) ); 
involving exponentials and 
logarithms 

Combine logarithms logcontract (jf) | logcontract (log(a)+1log(b)) ; 


The name of the radcan command arises from the fact that it converts 
expressions involving radicals, that is, roots, powers and logarithms into a 
canonical (or standard) form. 


These commands are demonstrated in the following activity. 
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5 Exponential functions and logarithms 


Simplify the expressions 
In(w? + 22? + 2) — In(a) 
In(x + 1) 
by following these steps. 


and In(a) + 21In(b) — In(c). 


(a) Define u to be the expression 
In(w? + 2a? + x) — In(x) 
In(a + 1) 


Remember that the Maxima command for a natural logarithm is log. 


(b) Simplify u using the following command. 


radcan(u) ; 


Can you simplify this expression and obtain the same answer by 


hand? 


(c) Define v to be the expression In(a) + 2In(b) — In(e). 
(d) Simplify v by using the following command. 
logcontract (v) ; 


If you began studying this section from Unit 3, Activity 48, now return to 
the unit to continue your study of functions. 
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6 Trigonometry 


The Maxima commands for the trigonometric functions are given below. 
Like many computer systems, Maxima assumes that all angles are 
measured in radians. Notice that the names of the inverse trigonometric’ 
functions sin~!, cos~! and tan~! in Maxima are asin, acos and atan 
respectively. This notation is used by many computer systems; it is short 


for arcsin, arccos and arctan, the alternative names often used for sin~', 
cos”! and tan7!. 
Trigonometric functions 

Function | Syntax | Example 

sin sin(™) | sin(1); 

cos cos(§) | cos(3*%pi/2) ; 
tan tan(§ | tan(%pi/4); 
cosec csc() | csc(2); 

sec sec(§) | sec(%pi/3) ; 

cot cot (ff) | cot (%pi/4) ; 

sin“! asin( | asin(sqrt (3) /2); 
cos”! acos({) | acos(1/sqrt(2)); 
tan-! atan(§) | atan(1/2); 
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Calculate the following using Maxima. 


1 «=(%) 


Remember that to enter the constant 7 in Maxima, you type %pi- 


Maxima returns an exact answer, in the form of a surd. Remember 
that Maxima always uses exact values where possible, rather than 
decimal approximations. 


(b) sin 45° 


| Remember to convert the angle to radians, by multiplying 45° by a 


(c) sec (3.4) 


Here, since the angle is entered as a decimal number, Maxima returns 


a decimal approximation to the answer. 


d) sin~! (2) 
2 
Maxima returns an exact answer, as a multiple of 7. 


(e) cosec (=) 


There is no simple exact answer here, so Maxima leaves this 
expression unevaluated. 


To find a decimal approximation, you can use the command 


float( csc(%pi/7) ); 


(f) tan (-3) 


Again, Maxima leaves this expression unevaluated, but it simplifies it 
slightly, using the fact that tan(—0@) = — tan é. 
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(g) cos“ (2.0) 


Here, since the argument of the function is a decimal number, 
Maxima returns a decimal approximation to the answer. 


You might be surprised that it returns an answer at all, since 2.0 is 


not in the domain of cos~!. 


Notice that the answer is a multiple of 4i. In Maxima, %i represents 
the imaginary number i whose square is —1, which you will learn 
about in Unit 12. Although there is no real number y such that 
cosy = 2, there is an imaginary number with this property, and 
Maxima returns this value. 


If you are working only with real numbers, you should ignore any 
Maxima answers that contain the constant %i. 


(h) cos! (2) 


This expression is equivalent to the expression in part (g), but does 
not contain a decimal number, so Maxima tries to find an exact 
answer. Since there is no simple exact answer, Maxima leaves the 
expression unevaluated. 


This hides the fact that the answer is not a real number. You need to 
take care when using inverse trigonometric functions. 


You can solve equations involving trigonometric functions in Maxima, 
using the solve command that you met in Subsection 3.2 of this Guide. 
However, Maxima will return only one solution to such an equation, even 
though many may exist. It warns you of this, as shown in the following 
activity. 


6 Trigonometry 
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(a) Solve the equation cos(x) — 0.6 = 0 by entering 
solve (cos(x)-0.6=0) ; 


Since the solve command always tries to find exact solutions, 
Maxima first converts the decimal number —0.6 to the exact 
fraction —3. It tells you that it has done this. 


It also displays the following warning to tell you that there might be 
other solutions. 


solve: using arc-trig functions to get a solution. 
Some solutions will be lost. 


Then it returns a single solution. 


(b) Find a decimal approximation to the solution obtained in part (a), by 
entering 


float (%) 


Remember, % represents the answer of the last evaluated calculation. 


(c) Solve the equation sin(3x +4) = 3. 
(d) Solve the equation cos(2x + 4) — sin(a/3) = 0. 


Maxima returns the original equation, slightly rearranged but 
unsolved. This is because the solve command could not find an exact 
solution to this equation. 


In circumstances such as Computer activity 32(d) when the solve 
command fails to find the exact solution of an equation, you can use 
another Maxima command, find_root, to try to find an approximate 
solution. 


Whereas the solve command manipulates the equation mathematically to 
find a solution, the find_root command performs a series of numerical 
calculations to find a decimal approximation to a solution. This is known 
as solving the equation numerically, or finding a numerical 
approximation to a solution. 
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You know from Unit 2 that the solutions of an equation of the form 
f(x) =0 are the z-coordinates of the points at which the graph of 
y = f(x) crosses the x-axis, as illustrated in Figure 10. 


y 


Figure 10 The graph of y = f(x) with the solutions of f(x) = 0 marked 


The find_root command uses the following property of functions. If f is a 
function and [a,b] is an interval included in its domain such that one of 
f(a) and f(b) is negative and the other is positive, then there is at least 
one solution of the equation f(a) = 0 in the interval (a,b). Figure 11 
illustrates this property for two different functions. This property holds as 
long as the graph of f is continuous on the interval {a, 6]; informally, this 
means that you can draw the graph of f for input values in the interval 

{a, b] without taking your pen off the paper. 


yw y 
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Figure 11 There is at least one solution of f(x) = 0 in the interval (a, 6) 
if f(a) and f(b) have different signs 


The find_root command looks for a single solution of an equation of the 
form f(z) = 0 that lies within an interval that you specify. The command 
requires that the values of f at each endpoint of the interval have different 
signs. If this is not the case, then Maxima gives an error message. 
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6 Trigonometry 


Solving an equation numerically 


Operation Command Example 

Solve an equation numerically | find_root (@xpression, | find_root (2*x*2-1,x,0,1); 
variable, 
interval start value , 
interval end value ) 


Note 1: the command finds a value of the variable at which the 
expression is zero. 


Note 2: the values that the expression takes at the endpoints of the 
interval must have different signs. 


(a) Find a numerical approximation to a solution of the equation 


cos(2a + 4) — sin (5) 30; 


as follows. 
(i) First, define the function f in Maxima, with the rule 


f(x) = cos(2a + 4) — sin (5) ? 


Remember, to define a function use :=. 


(ii) Next, plot the graph of f(x) over the interval —2 < x < 2, say. 
This will help you to find an interval that contains a solution. 


Remember, you can plot graphs using the wxplot2d command. 


From the graph, you can see that f(0) < 0 and f(1) > 0, so there is a 
solution in the interval (0,1). 


If there were no solution in the interval over which you plotted the 
function, then you could plot the function over a larger interval, and 
so on, until you find a suitable interval, if one exists. 
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(iii) Find a numerical approximation to the solution of the equation in 
the interval (0,1) by entering 


find_root (f(x), x, 0, 1); 


Remember to type the underscore symbol (_) in the middle of the 
command name. 


Notice that the first argument of the find_root command is an 

expression (or a function defined by an expression) not an equation. 
Maxima finds an input value at which the expression takes the 
value 0. 


(b) Using similar methods, find a solution of the equation sin?(« + 2) = 3x. 


Hint: Consider the function g(a) = sin?(x + 2) — 3x and find a value 
of x for which g(x) = 0. 


Remember that sin?(« +2) means (sin(a + 2))*. You can enter this 
expression in Maxima by typing either (sin(x+2))*2 or sin(x+2)*2. 


Maxima can use trigonometric identities to rewrite expressions involving 
trigonometric functions in different forms. 


These behave like the algebraic simplification commands you met in 
Section 3 of this Guide. In particular, you need to guide Maxima in how to 
conduct such rearrangements. Commands for doing this are demonstrated 
in Computer activities 34 and 35. 


(a) Enter 
trigreduce( sin(x)*2 ); 


The trigreduce command attempts to express powers of sines and 
cosines of x in terms of sines and cosines of multiples of x. It reduces 
the powers of the functions. 
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6 Trigonometry 


(b) Enter 
trigexpand( % ); 


Remember that % represents the result of the previous calculation. 


The trigexpand command attempts to express trigonometric 
functions of sums, differences and multiples of angles in terms of 
trigonometric functions of the individual angles. 


(c) Enter 
trigsimp( % ); 


The trigsimp command attempts to simplify expressions involving 
trigonometric functions. 


This returns you to the expression that you started with! 


G2 eiy SSS aes — ee ee 
Simplify ae ) by entering 


trigrat( sin(4+*x)/sin(x) ); 


The trigrat command attempts to simplify algebraic fractions 
involving trigonometric functions. 


Trigonometric manipulation commands 


Operation Command 


Expand trigonometric functions trigexpand(j) | trigexpand(sin(A+B)) ; 
of sums, differences and 
multiples of angles 

Express powers of sines and 
cosines in terms of sines and 
cosines of multiple angles 

Simplify trigonometric expressions 

Simplify algebraic fractions 
containing trigonometric 
functions 


trigreduce(§§) | trigreduce(sin(x)%2) ; 


trigsimp(sin(x)*2+cos(x)’2) ; 
trigrat ((sin(2x))/cos(x)); 


trigsimp(™§) 
trigrat(™) 


7 Plotting circles 


You can use these commands, and others that you learned previously, in 
the final activity of this section. 


(a) Find an expression for cos(5a) in terms of products of sin(x) and 
cos(xr). 

(b) Plot the graphs of cos(5z) and the expression you found in (a), to 
check that they look the same. 


If you began studying this section from Activity 37 of Unit 4, this now 
completes your study of the Unit. 


7 Plotting circles 


In this section you will learn how to use Maxima to plot circles. 


In Subsection 3.3 of this Guide, you learned how to use the wxplot2d 
command to plot curves. You can use this command to plot a curve only if 
the equation of the curve is in explicit form; that is, only if it is written in 
the form y = f(x), where f is a function. 


To plot a curve represented by an equation in implicit form, such as 
x? + y? = 1, you can use the wximplicit_plot command. This is 
demonstrated in Computer activity 37. 


The wximplicit_plot command is not part of the basic Maxima collection 
of commands, so you have to add it by loading an additional Maxima 
package, the implicit_plot package, which was installed on your 
computer when you installed Maxima. The Maxima command to load an 
additional package is load, which is also demonstrated in Computer 
activity 37. Each package needs to be loaded only once per Maxima 
session. 


If you are using a Maxima interface other than wxMaxima, then you 
should use the implicit_plot command instead of the wximplicit_plot 
command. You load and use it in exactly the same way as 
wximplicit_plot. 
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Plot the circle represented by the equation x? + y? = 1, as follows. 
(a) Load the implicit_plot package by entering 
load(implicit_plot) ; 


This package provides the wximplicit_plot command. 


(b) Plot the circle by entering 
wximplicit_plot (x‘2+y*2=1, [x,-2,2], [y,-2,2]); 


The first argument of the command is the equation to be plotted. 
The second argument specifies the variable to be plotted on the 
horizontal axis, and the range of values of that variable to be plotted. 
The third argument specifies the variable (and its range) to be 
plotted on the vertical axis. You must include this third argument. 


The circle with equation x? + y? = 1 has centre (0,0) and radius 1. 
So the circle lies within —2 < x < 2 and —2 < y < 2 and hence these 
seem suitable ranges to use. 


The curve that Maxima plots, however, does not look very circular! 
This is discussed below. 


Note that the first argument of the wximplicit_plot command is an 
equation, whereas the first argument of the wxplot2d command is an 
expression. 


The curve that you plotted in Computer activity 37 does not look circular 
because it has different scales on the axes. Whatever ranges of «- and 
y-values you specify for a plot, by default Maxima makes the y-axis appear 
3 as long as the a-axis. 


To specify equal scales, you can add the following additional argument to 
the wximplicit_plot command: 


{gnuplot_preamble, "set size ratio -1"] 


Gnuplot is the underlying software that Maxima uses to plot graphs. In 
general, if you set the size ratio to be —r, where r is a positive number, 
then the length that represents 1 unit on the y-axis is r times as long as 
the length that represents 1 unit on the z-axis. On the other hand, if you 
set the size ratio to be r, where r is a positive number, then the y-axis 
appears 7 times as long as the z-axis. 
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7 Plotting circles 


Plot the circle x? + y* = 1, using equal scales on the x- and y-axes by 
entering 


wximplicit_plot (x*2+y*2=1, [x,-2,2], [y,-2,2], [gnuplot_preamble,"set size ratio -1"]); 


Including the gnuplot_preamble argument in the wximplicit_plot 
command can be tiresome if you are plotting lots of curves and need 
equally scaled x- and y-axes for each. An alternative is to change the 
default behaviour of Maxima so that it uses axes with equal scales for the 
rest of your session. You will see how to do this in the following activity. 


Tf you have configured Maxima to change the font and line styles used 
when plotting graphs, as described in the Accessibility section of the OU 
Maxima website, then the command given in the following activity will 
revert these back to their original styles. To retain your graph style 
settings, use 


set_plot_option([gnu_preamble, OUEqualScales]) ; 
instead. 


Plot the circle (a — 1)? + y? = 3, using equal scales on the xr- and y-axes as 
follows. 


(a) Tell Maxima to always use equal scales on the x- and y-axes by 
entering 


set_plot_option([gnuplot_preamble, "set size ratio -1"])$ 


This sets [gnuplot_preamble, "set size ratio -1"] to bea 
default option for all plots during the rest of your Maxima session. 


The line ends with $ to prevent the display of the output of the 
command. When present this output consists of a long list of all the 
current default settings, which can safely be ignored. 


(b) Plot the circle by entering 
wximplicit_plot((x-1)%2+y*2=3, [x,-1,3], [y,-2,2]); 
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In the remainder of this section, it is assumed you have changed the 
default plotting behaviour as described in Computer activity 39. 


To plot more than one curve you can list their equations within square 
brackets as the first argument of wximplicit_plot, as you did in 
Subsection 3.3 of this Guide when using wxplot2d. The equations do not 
necessarily have to be in implicit form. This is demonstrated in the 
following activity. 


Plot the line y = 3a —5 and the circle (x — 2)? + (y — 3)? = 2 as follows. 


(a) For simplicity, assign the equation of the line to the variable line, and 
the equation of the circle to the variable circle. 


Remember, to assign an expression to a variable, use the colon (:) 
operator. So to assign the equation of the line to the variable line, 
use 


line: y=3*x-5; 


(b) Plot the line and circle on axes with equal scales by entering 
wximplicit_plot([line, circle], [x,0,5], [y,0,5]); 


The first argument is a list of the equations to be plotted. Since you 
have assigned the equations to variables, you can simply list the 
names of the variables. 


Since the circle has centre (2,3) and radius V2 ~ 1.41, the ranges 
0<a<5and0<y<5 are large enough for the whole circle to be 
plotted. 


The commands introduced in this section are summarised as follows. 


Loading additional packages 


Operation Command Example 


Load a package | load (package/mame)) | load(implicit_plot); 
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7 Plotting circles 


Plotting curves represented by equations in implicit form 


Operation Command Example 

Plot a curve represented | wximplicit_plot ((@qiation, wximplicit_plot(x*2+y’2=1, 

by an equation in implicit | @razentalyrangel, (x,-1,1], [y,-1,1]); 

form Werticalizange,...) 

Plot several curves wximplicit_plot( wximplicit_plot( 

represented by equations list of equations , [x42+2#y%2=4, y=x’3], 

in implicit form horizontal range , (x,-1,1], [y,-1,1]); 
Nerticanirange,...) 


Note 1: the first argument is an equation or a list of equations. 
Alternatively it can be an expression, or a list of expressions, in which 
case the equation formed by setting the expression equal to zero is 
plotted. 


Note 2: equations in explicit form can also be plotted using this 
command. 


Graph plotting options 


Option 


Change the relative lengths/scales 
of graph axes. 

(A positive number specifies the 
ratio of axis lengths, a negative 
number specifies the ratio of axis 
scales.) 


Example 


Argument 


[gnuplot_preamble, 

"set size ratio -1"] 
(which sets the axes to have 
equal scales) 


{gnuplot_preamble, 


"set size ratio Humber") 


Setting default plotting options 


Operation Command Example 


Set the default plotting options | set_plot_option(@pti6H) | set_plot_option( 
[gnuplot_preamble, 
"set size ratio -1"]); 


Circles are not the only type of curve that can be represented by equations 
in implicit form. In the next activity you will plot some other equations of 
this type and see the curves that they represent. 


59 


Plot the curves represented by the following equations, using equal axis 
scales for each curve. Make sure you use suitable ranges of x and y. 


Remember, you can use the editing facilities of wxMaxima to avoid 


having to type the whole wximplicit_plot command each time. 


(a) a@+4y2=1 (b) 2?@—4y?=1 8 (c) y2=a2?-2r+1 
(d) yY=22-22+5 8 (e) yYo=a?-c®  (f) (x? +y?- 1) = 27? 


If you began studying this section from Activity 16 of Unit 5, now return to 
the unit to continue your study. 


8 Differentiation 


The Maxima command for differentiating an expression is diff, which is 
demonstrated in the following activity. This command finds derivatives by 
using the derivatives of standard functions together with the sum rule, 
constant multiple rule, product rule, quotient rule and chain rule. 


Sometimes Maxima may give the result in a different, but equivalent, form 
to the expression that you might obtain when differentiating by hand. If 
this happens, then you might like to check that the two answers are 
equivalent. You can do this by finding the difference between the two 
answers and checking that it simplifies to zero. A rougher check is to use 
Maxima to plot graphs of the two answers and check they appear to be the 
same. 


8 Differentiation 


(a) Find the derivative of f(x) = a4 + 5x? + 15, by entering 
diff (x*4+5*x%2+15, x); 


The first argument of the diff command is the expression to be 
differentiated, and the second argument is the variable to differentiate 
with respect to. 


You must include the second argument. If you do not, then the 
answer returned by Maxima will be the derivative multiplied by 
del(x). This notation is outside the scope of MST124. 


As you would expect, the first argument can also be a variable that 
you have previously defined to be an expression, or a function you 
have previously defined. 


2 
(b) Define the variable p to be equal to aaa and then find its 
derivative with respect to t by entering 
diff(p,t); 


(Remember, to assign a value to a variable you use a colon (:).) 


Using diff in this way can help you check the expression to be 


differentiated has been entered correctly before you differentiate it. 


(c) Define the function g(u) = (u? + 3) In(u?), and find its derivative using 
diff (g(u) ,u); 


(Remember, to define a function you use :=, and the Maxima 
command for natural logarithms is log.) 


Notice that the first argument given here is g(u), not just g. 


(d) Find the second derivative of the expression p defined above, by using 
diff(p,t,2); 


The third argument of diff specifies the number of times that the 
expression is to be differentiated. 


(e) Find the third derivative of g(u) defined above. 
GRRE ee ee Ser 
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The diff command is summarised in the following box. You can practise 
using it in the next activity. 


Differentiation 


Operation Command | Example 
ditt (\§, HaetaBIE) 


ditt (lj, HaRaBIE, 
positive integer) 


Differentiate diff (sin(x*2) ,x); 


Differentiate 
multiple times 


diff (log(x) ,x,3); 


Use Maxima to differentiate the following, simplifying the answer when 
possible. 
Remember, you can enter e” in Maxima either as Ze“x or exp (x). 


(a) 22° + (b) In(e*+2) (ce) a 


You can use the diff command to assign the derivative of a function to 
another named function. This is useful when you want to evaluate the 
derivative at a particular input value, for example. 


The next activity shows you how to do this — unfortunately it’s not 
completely straightforward. In the activity, the original function is denoted 
by f and its derivative is denoted by df, rather than by £’. This is 
because you cannot use f’ as a function name in Maxima. 


sin(x) 
+1 
(b) Try to define df (x) to be the derivative of f(a), by entering 


df (x) :=diff (f(x) ,x); 


(a) Define the function f(2) = 
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(c) Try to find the value of the derivative of f at x = 3 by evaluating 
df (x) at x = 3; that is, try entering 


df (3); 


This gives an error. 
The cause of this error lies in the previous command: 


df (x) :=diff (£(x) ,x); 


The problem is that when you enter this command, Maxima does not 
perform the differentiation. Instead, it just notes that df (x) means 
diff (f(x) ,x). Then, when you enter df (3), Maxima understands 
this to mean diff (f (3) ,3), which is meaningless, as you cannot 
differentiate with respect to a particular number. 


To force Maxima to perform the differentiation, you can use the 
syntax '' (9) (which contains two single quotation marks, not a 
double quotation mark, and a pair of round brackets). This is 
demonstrated below. You need to use this syntax whenever you want 
to assign the result of a Maxima command to a function. 


(d) Define df (x) to be the result of differentiating f with respect to x, by 
entering 


df (x) :='' (diff (f(x) ,x)); 


Remember that '' is two separate quotation marks. Use the upright 
quotation mark on your keyboard ('), not the sloping one (~). Don’t 
forget to include the outer round brackets too. 


Notice that this time the output line shows that df (x) is defined to 
be the result of differentiating f(x). 


(e) Find the value of the derivative of f at 2 =3 by entering 
df (3) ; 


If you want a decimal approximation to the result, use the float 


command. 


(f) Use a similar approach to find the derivative of sin(x)e°° at x = 1, 
giving your answer as a decimal number to 3 significant figures. 


The command used to force the evaluation of the derivative in Activity 44 
is summarised as follows. 
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9 Integration 


Forcing the evaluation of a command 


Operation Command | Example 


g(x) :='' (diff (x*4,x)); 


Force the evaluation of a | '' () 
command 


An alternative approach to Computer activity 44 parts (d) and (e) is to 
assign the derivative of f to a variable rather than use the derivative to 
define a function. You can do this by entering, for example 


df :diff (f(x) ,x); 


This calculates an expression for the derivative of f at x and assigns it to 
the variable df (the ''(§§) command is not required). 


You can then evaluate the derivative at 3 by substituting the value 3 for x, 
which you can do by entering 


subst (3,x,df); 


If you began studying this section from Activity 22 of Unit 7, now return to 
the unit to continue your study. 


9 Integration 


The Maxima command to integrate an expression is integrate. This is 
demonstrated in the following activity. 


As with differentiation, the results that you obtain from integrating using 
Maxima might be in a different algebraic form to the expression that you 
might obtain when integrating by hand. To show that the two answers are 
equivalent integrals, you could find the difference between the two answers 
and check that it simplifies to a constant. A rougher check is to use 
Maxima to plot graphs of the two answers (taking the arbitrary constant 
to be zero, for example) and check they seem to be vertical translations of 
each other. 
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(a) Evaluate / 2x? + €%* da by entering 


integrate (2*x’2+exp(3+x), x); 


The first argument of integrate is the expression to integrate, and 
the second is the variable to integrate with respect to. 


The output returned by Maxima does not include an arbitrary 


constant. That is, Maxima gives an antiderivative of the expression, 
not its indefinite integral. When you write down an indefinite integral 
that you have found using Maxima, you need to remember to add an 
arbitrary constant yourself. 


3 
(b) Evaluate the definite integral | = de by entering 
2 


integrate(1/x, x, 2, 3); 


The third and fourth arguments of integrate are the lower and 
upper limits of integration respectively. 


Notice that, as usual, Maxima returns an exact result. Use the float 
command to find a decimal approximation if required. 


(c) Define the function p(x) = nem and define q(x) to be an 
antiderivative of p. 


If you try to define q(x) using 


q(x) :=integrate (p(x) ,x); 
then Maxima does not actually do the integration, but just notes that 
q(x) means [r(0) dr. 
As with differentiation, to force Maxima to perform the integration 
you need to use the '' (§§) command and enter 


q(x) :='' (integrate (p(x) ,x)); 
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Integration 
Operation Command Example 
Integrate integrate ([), Wariable)) | integrate(x*2,x); 


(find an antiderivative) 


Evaluate a definite integral | integrate (),Wariable), | integrate(sin(x),x,0,1); 


You can practice using the integrate command in the following activity. 


Use Maxima to find the following. 
(a) [e-a0e- 1) dx (b) Iz dt (c) if ef sint dt 
> a aaa 


It is not always possible to express the result of integrating a function in 
terms of the mathematics functions you have previously met. The 
following activity illustrates some of the things that you may encounter 
when integrating functions using Maxima. 


Try to calculate the following using Maxima. 


(a) femar ) fewer ©) [Se 


Maxima was unable to find an antiderivative for part (a), and so 
returned the integral uncalculated. In parts (b) and (c) answers were 
returned containing the function erf, which is known as the Error 
function, and the function gamma_incomplete, the incomplete 


gamma function. These are examples of what are known as special 
functions in mathematics. They cannot be simply expressed in 
terms of functions you have already met, but many can be defined in 
terms of antiderivatives of simpler functions. 


Sometimes you may want to integrate a function that contains a variable 
other than the one you are integrating with respect to. For example, the 


integral | 2" da contains the variable n as well as the variable z. When 


you use Maxima to find such an integral, you sometimes need to provide 
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information about the other variable, as demonstrated in the following 
activity. 


(a) Try to use Maxima to find the integral / x" dx, as follows. 


(i) First, ensure that n has no predefined value, by entering 
kill(); 


(ii) Then, enter 


integrate(x%n, x); 


Maxima asks the question: Is n+1 zero or nonzero? 


This is because the answer to / a" dx depends on the value of n. 


Ifn = —1 (son+1=0), the integral is [5a which is equal to 


log |x| +c, whereas if n is any other value, then the integral is equal 
gntl 


n+1 
knowing more about n. 


to 


+c. So Maxima cannot calculate the integral without 


Enter the answer n; 
(which is short for nonzero) 


You could also answer the question with z; 


to mean zero. 


Maxima displays the expected result. 


(b) An alternative approach is to give Maxima information about n before 
you use the integrate command. You can do this as follows. 


(i) Tell Maxima that n # —1 by entering 


assume (notequal(n,-1)); 


The assume command tells Maxima information about variables. For 
example, assume(a>0) ; tells Maxima that a is positive. 


Maxima has no simple symbol for ‘not equal to’, so you have to use 
the notequal command for this, as above. 
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(ii) Now integrate x” using 


integrate(x“n, x); 


No question is asked this time, since Maxima knows all it needs to 
know about n. 


(iii) Enter 
facts(); 


to see all the additional facts about variables known to Maxima. 


This lists all such facts. The only one shown should be the one you 


gave earlier. 


(iv) Tell Maxima to forget about the assumption on n by entering 


forget (notequal(n,-1)); 


(v) Type facts(); again to see all the facts now known. 


No facts are displayed, since Maxima has now forgotten the 
assumption about n. 


The commands introduced in the previous activity are summarised below, 
together with the equal command, which is used in a similar fashion to 
notequal. 


Facts about variables 


assume (jf) 


assume (a>0) ; 


Make an assumption about a 
variable 


State two things are equal equal()), assume (equal (n,3)); 
State two things are not equal | notequal (/),)) assume (notequal(n,-1)); 
Forget a property forget (property) forget (a>0); 

List all known facts facts() facts(); 


List all known facts about facts(a); 


a particular variable 


facts (Wardable)) 
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Calculate the integral 


1 
if ErunAt where a > 0. 


You saw in Computer activity 47 that Maxima cannot always find an 
antiderivative for a function. Even if Maxima cannot find an antiderivative 
of a function f, it can still find an approximate value for a definite integral 


b 
of the form i f(x) de. 
a 


In Unit 8, you saw how the value of the definite integral can be 
approximated by breaking the area between the graph of the function and 
the z-axis into a number of rectangles. Maxima uses a more sophisticated 
version of this basic method to calculate approximate values of definite 
integrals. 


You will learn the Maxima command for this in the following activity. 


1 
(a) Try to calculate | In(cos(x?)) dx by entering 
0 


integrate(log(cos(x%2)), x, 0, 1); 


The integral is returned unevaluated, since Maxima cannot 


calculate it. 
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10 Matrices 


(b) Find an approximate value of the integral by entering 
quad_qags(log(cos(x*2)), x, 0, 1); 


The quad_qags command calculates an approximate value for a 
definite integral. 


The ‘quad’ part of the name of this command arises from the fact the 
process of finding the value of a definite integral (or any area) is 
sometimes called quadrature. The letters q, a, g and s in the second 
half of the name specify the particular method Maxima uses to find 
the approximate value of the definite integral. This method gives 
good results for a wide variety of integrals. 


Notice that the output of the quad_qags command is a list. The first 
element of the list is the approximate value of the definite integral, the 
second element is an estimate of the accuracy of the approximation, 
the third element is the number of values at which the function to be 
integrated was evaluated during the calculation, and the final element 
is an error code. An error code of 0 means that no errors occurred. 


1 
(c) Find an approximate value for the definite integral ' edt. 
0 


Finding approximate values for definite integrals 


Operation 
Find an approximate value | quad_gqags (jj, Wariable} , quad_qags(exp(-x*2) ,x,0,1); 
of a definite integral ‘lower limit , upper limit ) 


If you began studying this section from Activity 47 of Unit 8, now return to 
the unit to continue your study. 


10 Matrices 


In this section, you will learn how to use Maxima for matrices. 


10.1 Matrix algebra 


There are several ways to input matrices in Maxima. Some of these are 
demonstrated in the following activity. 
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(a) Assign the matrix 


13 5 2 a 
iE O a to the variable A by entering 


A:-matrixt [1,335], [2;4,76) J 


The matrix command allows a matrix to be defined. You give each 
row of the matrix as a list (that is, within square brackets with each 
element separated by a comma). 

Notice that wxMaxima displays matrices using square brackets rather 
than the round brackets used in the unit: 


13 5 
2.6)" 
If you are using a different interface for Maxima, or have chosen 
different display settings, matrices may be displayed as follows. 
[tee ey ator 


C ] 
2 446 J 


(b) If you are using wxMaxima, you can also input a matrix using an 
on-screen form. Use this form to assign the matrix 
1234 
4 1 2 3) to the variable B as follows. 
2413 


(i) From the Algebra menu select Enter Matrix... 


Make sure you select Enter Matrix..., not Generate Matrix... 


A window opens requesting properties of the matrix, as shown below. 
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(ii) The matrix given above has size 3 x 4, so enter 3 as the number 
of rows and 4 as the number of columns, replacing the default 
values shown. 


It is also possible to give information on the type of matrix. You have 
not learned about these types in MST124, so leave this set to 
general. 


Finally, enter B as the variable name the matrix will be assigned 
to, and click ‘OK’. 


This opens a second window requesting the elements of the matrix. 


(iii) Complete entry of the matrix by filling in the elements in this 
second window. When you have finished, click ‘OK’. 


This has the effect of including a matrix input line in your worksheet, 
followed by the appropriate output. 


Entering matrices 


Operation Command Example 


Specify a matrix | matrix (Ow, #6W,...) | A:matrix([1,2],[3,4]); 


Note: matrices can also be entered using the Enter Matrix... option 
of the Algebra menu. 


Once a matrix has been input, you can obtain its size and individual 
elements by using the commands demonstrated in the following activity. 
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(a) Find the size of the matrix A defined in Computer activity 51 by 
entering 


matrix_size(A); 


The output is a list containing the number of rows and the number of 


columns of the matrix. 


(b) Extract the element in the third row and second column of B by 
entering 


B[3,2]; 


Alternatively, B[3] [2] gives the same result. 
Notice the use of square brackets to reference elements of a matrix is 


the same syntax as used in Section 3.2 of this Guide to obtain 
individual elements of a list. 


(c) Extract the third row of B by entering 
B[3] ; 
(d) Change the element in the third row and second column of B to be 0 
by entering 
B[3,2]:0; 


then display the revised matrix by entering 


B; 


You can use commands of this type to correct any errors that you 
might make when entering a matrix. 


Using matrices 


Operation 


Size of a matrix matrix_size (matrix) 
Row of a matrix matrix [row] A(2]; 
Element of a matrix | Matrix [iow , eoluma) A[3,4]; 
or matrix [row] (eotuma) | A(3] [4]; 


The syntax to use for matrix operations in Maxima is given in the 
following box. You can practice using them in Computer activity 53. 


matrix_size(A); 
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Matrix operations 


Operation Syntax | Example 


Matrix addition A+B; 
where A and B are matrices of the same size 
A-B; 

where A and B are matrices of the same size 
2*A; 

-B* 

where A and B are matrices of appropriate sizes 
AMS; 

where A is a square matrix 


Matrix subtraction 


Scalar multiplication 
Matrix multiplication 


Matrix powers 


(a) Define the following matrices in Maxima 
2 0 0 3 13 
PH "| OSU Qo) ae ris 
-4 6 1 3 


b) Find P+Q. 


The matrices P and Q are the same size, so this calculation is 
possible. 


L 
(c) Find P—R. 


r 


The matrices P and R are not the same size, so this calculation is not 
possible. Maxima gives an error message: 


fullmap: arguments must have same formal structure. 
-- an error. To debug this try: debugmode(true) ; 
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(d) Find 2P. 


This is a scalar multiple of a matrix, so use the * operator. 


(e) Find PR. 


This is matrix multiplication, so use the . operator. 


The sizes of these matrices means the calculation is possible. 


(f) Find PQ. 


Remember to use the . operator. 


This calculation is not possible due to the sizes of the matrices. 
Maxima gives an error message: 


MULTIPLYMATRICES: attempt to multiply nonconformable 
matrices. 
-- an error. To debug this try: debugmode(true) ; 


(g) Find PS. 


The matrices P and S have sizes 3 x 2 and 2 x 1, respectively, so this 
is possible and the result is a 3 x 1 matrix, that is, a column vector. 
= 


(h) Find PT. 


The matrices P and T have sizes 3 x 2 and 1 x 2, respectively, so this 
product should not be possible. However, Maxima does return an 
answer! 


Maxima does not distinguish between single row matrices (which are 
sometimes called row vectors) and single column matrices (column 
vectors). It considers them interchangeable as circumstances require. 
You should take care to multiply only appropriately sized matrices. 


(i) Find R?. 


Since R is a matrix use the matrix power operator “”. 
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(j) Enter 
RSQ 


The result is not the same as in part (i). 


Here the * symbol was used and Maxima squared each individual 
element of R. 


(k) Enter 
PQ; 


The result is not PQ, since this product does not exist! 


Here the * symbol was used and Maxima multiplied corresponding 
elements of P and Q together. 


If you began studying this subsection from Activity 13 of Unit 9, now 
return to the unit to continue your study. 


10.2 Determinants and inverses 


The Maxima commands for finding the determinant and inverse of a 
matrix are demonstrated in the following activity. 


Consider the matrices 


120 x ee aa 
A=|-4 4 2 and B=1{5 6 2]. 
— a ak 221 


(a) First, enter the matrices A and B into Maxima. 
(b) Find the determinant of A by entering 


determinant (A) ; 


The determinant command calculates the determinant of a square 
matrix. 
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10 Matrices 


(c) Find the inverse of A, and assign the result to the variable C by 
entering 


C:invert (A) ; 


The invert command calculates the inverse of a square matrix. 


An alternative way to specify the inverse of the matrix A in Maxima 
is A” (-1), which accords with the notation usually used to denote a 
matrix inverse. 


(d) Check that C really is the inverse of A by calculating the products AC 
and CA. 


Remember to use the symbol . for matrix multiplication. 


Both products should be the 3 x 3 identity matrix. 


(e) Try to find the inverse of B. 


This gives an error: 


expt: undefined: 0 to a negative exponent. =a 
error. To debug this try: debugmode(true) ; 


This matrix does not have an inverse. You can check this by 
calculating its determinant. 


Another command that is sometimes useful when you are working with 
matrices is ident, which specifies an identity matrix. For example, 


ident (2) gives 3 Hp the 2 x 2 identity matrix. 


Matrix operations 


Operation Example 


Determinant (of a square matrix) | determinant (atrax) | determinant (A); 
Matrix inverse (of a square matrix) | invert ((Hatrix) invert (A) ; 

or matrix’ (-1) | a*(-1); 
Identity matrix ident ((Bize) ident (2) ; 
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For each of the following matrices, check whether it is invertible. If it is, 
find the inverse. 


2-11 pepe A se 

(a) {2 01 (b) 2 
A Oo 20 1-1 
-1 2 1 -1 


If you began studying this section from Activity 24 of Unit 9, now return to 
the unit to continue your study. 


11 Sequences and series 


In this section you will learn how to use Maxima to work with sequences 
and series. 


11.1 Plotting graphs of sequences 


In Subsection 3.3 of this Guide, you learned how to use the wxplot2d 
command to plot the graphs of algebraic expressions. Here, you will see 
how to use the same command to plot graphs of sequences. First, in 
Computer activities 56 and 57, you will learn how to define and work with 
sequences in Maxima. 


Consider the sequence (a,,) specified by the following closed form and 
range of values of n: 


a,=058° (n— 15,2532.) 


(a) Define this sequence in Maxima, by entering 


a[n] :=0.8%n; 


This command defines the general term a, of the sequence using the 
same define operator (:=) as is used when defining functions. 


Notice that you type the subscript n of the general term within square 
brackets, but the Maxima output uses a properly formatted subscript. 


11 Sequences and series 


78 


(b) Find the term ajo of the sequence by entering 
a[100]; 


Notice that when you defined the sequence in Computer activity 56 you’ 
did not need to enter the range of values of n. In fact, Maxima places no 
restriction on the values of n that you can use when calculating terms in a 
sequence defined by a closed form. For instance, for the sequence (ap) 
defined in Computer activity 56 it is possible to calculate values for 
a[3.5] and a[-4] using Maxima. 


However, these numbers are not terms of the sequence (a@,), which are 
defined only for integer values of n greater than or equal to 1. 


Consider the sequence (b,) specified by the recurrence relation 
bb =1, bh =2bp-it+1 (n=2,3,4,...). 


(a) Define this sequence in Maxima by entering the two commands below. 
(ot fel oe 
b[n) :=2*b[n-1] +1; 


Notice that you assign the value of the first term, b[1], using the : 
operator that is used to assign values to variables, and you specify the 
recurrence relation using the := operator that is used to define 
functions. 


Unlike for a sequence defined using a closed form, Maxima can 
calculate a term of a sequence defined using a recurrence relation only 
when the term number is an integer greater than or equal to the 
starting value, which in this case is 1. 


(b) Find the 5th and 100th terms of this sequence. 
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11 Sequences and series 


(c) Change the value of b; to be 10, by entering 
b[1] :10; 


and try to find the revised value of bs. 


Notice that the value of b; displayed is unchanged. This is because 
when Maxima uses sequences, it remembers all the values of the 
sequence that it has previously calculated, and bases subsequent 
calculations on them. This improves the efficiency of the calculations. 
However, it also means that sequences cannot be easily changed. 


To change the definition of a sequence, you first need to remove the 
sequence from Maxima’s memory using the kill command, and then 
redefine it. 


(d) Revise the sequence (b,,) correctly, by first entering 
kill(b); 
and then entering the two commands 
b[1]:10; 
bin) :=2*b [n-1] +1; 
to specify the revised sequence. 
Now calculate the value of bs. 


The commands for defining and working with sequences are summarised as 


follows. 
Sequences 
Operation Command | Example 
Define a sequence ‘sequence [n] :=expression in n a(n] :=n*2; 
using a closed form 
Define a sequence Sequence (FiFSt|LSrRIRUAUEE] - GIEHaNyWaNue | > (1) :1, 
using a recurrence system | Sequence [n] :~expression/in)/sequence(n=1]}) | b[n] :=2*b[n-1]; 
Calculate a term in ‘sequence [term number ] b[50] ; 
a sequence 


Note: to change the definition of a sequence, you first have to delete it 
using kill. 


In the following activity, you will learn how to plot a graph of a sequence. 
Remember, a graph of a sequence (a,,) is a set of points (n,a,,) plotted for 
a suitable range of values of n. 
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Plot the graph of the sequence (c,) specified by 
=n? (n=1,2,3) 
as follows. 
(a) First, define the sequence in Maxima by entering 
c[n] :=n*2; 
(b) Create a list of the coordinates of the points to be plotted; that is, the 


points (1,c1), (2,c2) and (3,c3) and assign it to a variable, say, pts 
(for ‘points’), by entering 


pts:( (1,c(1]], (2,c(2]), (3,c(3]] 1; 


Notice that you type the coordinates of each point as a list in square 
brackets. The variable pts is a list of lists. 


Do not use the variable name points for your list of points — this is a 
word used by Maxima to set plot styles, as you will see next. 


(c) Plot the graph of the sequence using the command 


Here, instead of the first argument of the wxplot2d command being 
an expression, or a list of expressions to plot, it is a list containing 
two elements. The first element is the keyword discrete. This tells 
Maxima that a set of discrete points are to be plotted. The second 
element is the list of the coordinates of the points to be plotted. 


The second argument of the wxplot2d command is 
[style, points]. 
This tells Maxima to simply plot points, and not join them with lines. 


The final two arguments specify labels for the x- and y-axes. These 
arguments were introduced in Subsection 3.3 of this Guide. 


You can also use other standard arguments of wxplot2d that you met 
| a For example, including [x,0,4] changes the range of the 


horizontal axis to 0 < x < 4. 
—_i 


The use of wxplot2d for plotting sequences is summarised as follows. 
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wxplot2d([discrete, pts], [style, points],[xlabel, "n"], [ylabel, "c[n]"]); 
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Plotting graphs of sequences 


Operation Command Example 
Plot a sequence | wxplot2d([discrete, MiStlet|eoordinatelpairs] , | wxplot2d([discrete, 
(style, points],...) ((1,5], [2,10]]], 


[style, points]); 


Graph plotting options 


Option | Argument | Example 


ae 


Plot points | [style, points] | [style, points] 


In Computer activity 58, the sequence plotted was short, so it was 
relatively easy to create the list of coordinates to plot. For a longer 
sequence this can be laborious, but Maxima can help by creating the list 
for you. This is demonstrated in the next activity. 


Plot the graph of the sequence given by the recurrence relation 


dy=0, dn=3(2—dn1) (n=2,3,4,. 


as follows. 


15) 


(a) First define the sequence (d,) in Maxima. 
(b) To plot the sequence, we need to plot the points 
(1,d1), (2,d2), (3,d3), ..-, (15, dis). 
Create a list of these points. and assign the list to the variable pts, 
using the command 
pts:makelist([n,d[n]], n, 1, 15); 


The makelist command creates a list. The first argument of the 
command is the general form of the elements of the list. In this case, 
each element has the form [n,d[n]]; that is, it is itself a list 
containing two elements: the term number and the value of the 
corresponding term. These are the coordinates of a point to be 
plotted. Here, the general form is expressed using the dummy 
variable n. 


The second argument of the makelist command is the name of the 
dummy variable used, and the third and fourth arguments are the 
minimum and maximum values taken by the dummy variable. 


When you create a long list you may not want to see all the output! 
Remember that you can end a command with a dollar sign ($) rather 
than a semi-colon (;) to suppress the display of the output. 
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(c) Plot the points in the list pts by using the command 
wxplot2d([discrete, pts], [style, points], [xlabel, "n"], [ylabel, "d[n]"]); 


as seen earlier. 


The makelist command used in the previous activity is summarised as 
follows. 


Creating lists 


| Operation Command 


Create a list | makelist (GeHeralterm, (@immylvariablé|, makelist(2*n,n,1,100); 
‘start value, end value) 


As you know from your study of Unit 10, a sequence can have infinitely 
many terms. Such a sequence is called an infinite sequence. Since 
Maxima cannot plot an infinite number of points, when plotting a graph of 
an infinite sequence, you need to restrict the number of points to a large, 
but finite, number. Try to choose a number that is large enough to reveal 
the long-term behaviour of the sequence. 


Use Maxima to plot graphs of the sequences specified below, by restricting 
the range of values of n to be n = 1,2,3,...,50 in each case. 


How do you think each sequence will behave as more and more terms are 
considered? 


(a) n=1, te=1.5rn-1—O0.1(r_-1)? (n=2,3,4,...) 
(b) sn =n®-10 (n=1,2,3,...) 
(Clit oe) 


If you began studying this section from Activity 17 of Unit 10, now return 
to the unit to continue your study, and in particular learn about the 
long-term behaviour of sequences. 
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11.2 Summing series 


You can calculate the sums of various series using the Maxima command 
sum which is demonstrated in the following activity. 


(a) Find the sum of the squares of the first 100 integers, that is 
100 
u” 
n=l 
by entering 
sum(n*2, n, 1, 100); 


The first argument of the sum command is the general term of the 
series to be summed, which in this case is expressed in terms of the 


index variable n. The second argument is the index variable itself and 
the third and fourth arguments are the lower and upper limits of the 
summation. 


(b) Write the sum 


tl fe ALN? nh 
3+(3) +) +-*(3) 
in sigma notation, and then use the Maxima sum command to 
calculate it. 
= 1 
(c) Try to calculate the sum of the infinite series ay mo by entering 


n=1 


sum(1/(4%n), n, 1, inf); 


Maxima uses inf to represent oo. 


Notice that Maxima displays the sum (unevaluated) using sigma 
notation. By default, Maxima does not attempt to calculate the sums 
of infinite series. 


You can, however, tell Maxima to simplify the result of a summation, 
which can sometimes result in it finding a value for a sum of an 
infinite series, if it exists. To do this, you change the value of the 
system variable simpsum, as shown below. (simpsum is short for 
‘simplify sums’.) 


Sequences and series 
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(d) Set the system variable simpsum to have the value true by entering 
simpsum: true; 


ox 
a! 
Now try to calculate SS me again. 


n=1 


Maxima now returns the sum of the series. 


Cae | 
(e) Find »X a 
= 
Gee Ss > ene eS eae Re 


The commands introduced in Computer activity 61 are summarised as 
follows. 


Summing series 


Sia 


simpsum: true; 


sum(i/n*2,n,1,inf) ; 


If you began studying this section from Activity 32 of Unit 10, now return 
to the unit to continue your study. 


12 Taylor polynomials 


You could use Maxima to find a Taylor polynomial for a function by using 
the diff command to calculate the required derivatives and then 
combining these to form the required polynomial. However, Maxima also 
has a command for calculating Taylor polynomials more efficiently. This is 
demonstrated in the following activity. 


12 Taylor polynomials 
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12. Taylor polynomials 


(a) Find the cubic Taylor polynomial about 7/6 for the function 
f(x) = cosa, by entering 


taylor(cos(x), x, %pi/6, 3); 


The taylor command calculates a Taylor polynomial for the function 
that is its first argument. Its second argument specifies the input 
variable of the function. 


The third argument of the command is the point about which the 
Taylor polynomial is to be found, which is also known as the centre of 
the polynomial. The final argument of the command is the degree of 
the Taylor polynomial required. 


Notice that the Maxima output ends with ..., which indicates that 
more terms could have been found. Also, the output line number is 
followed by /T/. This indicates that the output is a (truncated) 
Taylor series. 


(b) Define p to be the quintic Taylor polynomial about 0 for the function 
f(x) = tana, by entering 


p(x):=''(taylor(tan(x), x, 0, 5)); 


Here, the '' (§§) command is needed to force Maxima to calculate the 
Taylor polynomial which is then used to define p(x). 


Without the '' (§§) command, p(x) would be given by the expression 
taylor(tan(x), x, 0, 5). So if you tried to find, for example, 


p(3), then Maxima would interpret this to mean expression 
taylor(tan(3), 3, 0, 5); but this is meaningless and would 
therefore produce an error. 


You first met the ''(§§) command in Section 8 of this Guide. 


(c) Express the Taylor polynomial p(x) found above in terms of the 
variable y by typing 


p(y); 


As soon as a Taylor polynomial is used in a further calculation or 
command, the trailing ... is dropped, and the Taylor polynomial is 
converted to a rational expression. This is indicated by the marker 
/R/ following the output line number. 
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13 Complex numbers 


(d) Find the decimal value of p(x) when x = 0.1. Also find the difference 
between p(x) and tanzx at z= 0.1. 


Maxima warns you it has replaced certain decimal numbers with 
rational approximations when calculating these values. 


The taylor command is summarised as follows. 
Taylor polynomials 


Operation Maxima Command Maxima Example 


Find a Taylor polynomial | taylor (#ametion, Wariable), | taylor(sin(x),x,0,3); 
(centre), degree) 


Consider the function f(x) = e™(**). 


Plot a graph of the function and the graph of its quintic Taylor polynomial 
about x = 1 on the same axes over the range 0 < x < 2. 


(Remember, you learnt how to plot more than one function on the same 
graph in Subsection 3.3 of this Guide.) 


If you began studying this section from Activity 21 of Unit 11, now return 
to the unit to continue your study. 


13 Complex numbers 


In this section you will learn how to use Maxima to work with complex 
numbers. 


13.1 Working with complex numbers 


The imaginary number 7 is represented in Maxima by %4i. You saw this in 
Subsection 3.2 of this Guide, when you solved equations that have complex 
solutions as well as real solutions. The usual arithmetic operations work 
with complex numbers in Maxima, and there are commands for finding the 
real and imaginary parts of a complex number, its conjugate, its modulus 
and principal argument, as demonstrated in the following activity. 
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(a) Assign the complex number 3 + 4i to the variable a, by entering 


a:3t4*hi; 


When Maxima displays a complex number, it usually arranges it with 
the term in 7 as the first term, unless the imaginary part is negative 
and the real part is positive. 


(b) Similarly, assign the value 2 — 2i to b. 
(c) Calculate a+b, a-b, a*b and a/b. 


The sum and difference are displayed in simplified form, but the 
product and quotient are not. 


To simplify the result of a calculation and display it in Cartesian form 
(or rectangular form, as Maxima calls it) use the command 
rectform. 


(d) Find a*b in Cartesian form, by entering 


rectform(a*b) ; 


Find a/b in Cartesian form similarly. 
(e) Find the real and imaginary parts of a by entering 
realpart (a); 
and 


imagpart (a) ; 


(f) Find the conjugate of a by entering 


conjugate (a) ; 


(g) Find the modulus and the principal argument of b by entering 
abs(b); 
and 
carg(b) ; 


respectively. 
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13 Complex numbers 


Notice that the abs command is the same command that you used 
earlier to find the magnitude (also called the modulus, or the absolute 
value) of a real number. 


The name of the carg command is short for ‘complex argument’. The 
argument calculated is the principal argument, that is, the argument 
in the interval (—7, 7]. 


As you saw in Unit 12, any non-zero complex number can be expressed in 
polar form r(cos@ + isin @), where r is the modulus and @ is one of its 
arguments. The next activity uses this form. 


(a) Assign the complex numbers 3 (cos (57/13) + isin (57/13)) and 
4 (cos (77/13) + isin (77/13)) to the variables c and d respectively. 


(b) Use Maxima to find the moduli and principal values of c*d and ¢/d. 


Hint: use the trigreduce command to simplify the initial answers 
given. 


The commands introduced in the activities above are summarised as 
follows. 


Working with complex numbers 


Operation Maxima Command | Maxima Example 


Real part of a complex number realpart (() realpart (1+2*%i) ; 
Imaginary part of a complex number imagpart (/)) imagpart (1+2*%i) ; 
Conjugate of a complex number conjugate (jl) conjugate (1+2*%i) ; 
Modulus of a complex number abs (J) abs (1+2*%i) ; 
Principal argument of a complex number carg() carg(1+2*%i) ; 


Express a complex number in Cartesian form | rectform(\§) rectform(2/%i) ; 
If you began studying this section from Activity 27 of Unit 12, now return 


to the unit to continue your study. 


13.2 Solving polynomial equations 


In Subsection 3.2 of this Guide you met the solve command that is used to 
find solutions of equations. This command tries to find all the solutions of 
the equation (if they can be found using algebraic manipulation), including 
solutions that are non-real complex numbers. You might remember that in 
Computer activity 20(e) you solved an equation with non-real solutions. 


(a) Use the solve command to find the solutions of the equation 
H24+24+1=0, 


and assign the result to the variable solns. 


Notice that the solutions are given exactly, although in a rather 
unwieldy form! 


It is often more useful to find decimal approximations to solutions, 
expressed in Cartesian form, which you can do using the float and 
rectform commands. 


(b) Express the solutions found in part (a) in Cartesian form, with each 
real number, and the real and imaginary parts of each complex 
number, correct to two decimal places. 


(c) Find, in Cartesian form, all the solutions of 
24 — 2023 + 1712? — 6262 + 962 = 0. 


You may have noticed that in Computer activity 66 all the non-real 
complex solutions occur in complex conjugate pairs. In other words, for 
each of the equations, if a particular non-real complex number z is a 
solution, then so is its complex conjugate Z. You saw in Unit 12 that this 
property holds for any quadratic equation that has real coefficients. In 
fact, it holds for any polynomial equation that has real coefficients. It 
doesn’t hold for polynomial equations that don’t have real coefficients, 
however. For example, the solutions of the polynomial equation 


2-2 =0 


are 1+7 and —1 —i, which are not complex conjugates of each other. 


a3 


Complex numbers 
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Plotting solutions in the complex plane 


It is sometimes informative to plot the solutions of an equation in the 
complex plane. You can do this in Maxima using the methods that you 
learned in Subsection 11.1 of this Guide for plotting the discrete points 
that form the graph of a sequence. 


To plot the solutions of an equation, you first need to create a list of the 
points to be plotted, with the first coordinate of each point being the real 
part of a solution, and the second coordinate being the corresponding 
imaginary part. 


The following activity shows how you can do this. 


(a) Solve the equation 
z4+4z+5=0 


and assign the list of solutions to the variable s. 


Each element of the list s is an equation of the form z=value, where 


value is a solution of the equation. 


(b) Create a second list containing just the solutions of the equation 
(without the z=) and assign it to the variable v, by entering 


v:makelist(rhs(s[k]), k, 1, length(s)); 


The k*” element of the first list s is s[k], so rhs(s[k]) is the 
right-hand side of this element, which is the k** solution of the 
original equation. 


You met the makelist command in Subsection 11.1 of this Guide. 


The index variable k starts with the value 1 and ends with the value 
given by length(s), which is the number of elements in list s. This 
ensures that all of the the solutions of the equation are included in 
the list v. 


13 Complex numbers 
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13 Complex numbers 


(c) Create a list of the coordinates of the points to be plotted and assign 
it to the variable pts, by entering 


pts:makelist([realpart (v[k]) , imagpart (v[k])],k,1,length(s)); 


The kt" solution of the original equation is v[k]. The point in the 
complex plane corresponding to this solution has coordinates 


(Re(v 0k] ), Im(v [k])) , 
which can be expressed in Maxima syntax as 


(realpart(v[k]), imagpart(v[k])]. 


The command above uses makelist to create a list of such points, 
one for each solution of the equation. 


(d) Finally, plot the points found by entering 
wxplot2d([discrete,pts] , [style,points] , [xlabel,""], [ylabel,""]) 


This command plots a set of discrete points. The first two arguments 
of this wxplot2d command are the options that you used to plot 
graphs of sequences in Subsection 11.1 of this Guide. The final two 
arguments specify that the axis labels are to be left blank. 


Notice that since the complex solutions of polynomial equations with 
real coefficients occur in conjugate pairs, the points plotted are 
symmetric about the horizontal axis. 


Computer activity 67 introduced the length command, which is 
summarised as follows. 


Length of a list 


Operation Command Example 


Find the length of a list | length (G8) 
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(a) Find the seventh roots of unity by solving the equation 
2=1. 
Give the solutions in Cartesian form, with each real number, and the 


real and imaginary parts of each complex number, correct to two 
decimal places. 


(b) Plot the solutions in the complex plane. 
est Sa ee i eee’ 


If you began studying this section from Activity 35 of Unit 12, now return 
to the unit to continue your study. 


13 Complex numbers 
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Maxima accessibility guide 


This section outlines issues associated with the use of Maxima if you 
interact with a computer using only a keyboard, use a screenreader or 
require different colour settings. The information should help you decide 
whether the wxMaxima interface is suitable for you, or whether you should 
consider using an alternative interface, such as the command-line interface. 
It also explains how to configure Maxima to make it easier for you to use 
and contains an alternative version of Section 2 describing how to use the 
command-line interface. 


Using wxMaxima with a keyboard alone 


The majority of the wxMaxima interface is accessible using a keyboard 
alone. On a Windows computer you can access the wxMaxima menus by 
pressing and releasing the key, navigating through the menus using the 
keyboard arrow keys, and then pressing (Enter) to make your selection. 
Alternatively, you can press and hold tot key, and then press the 
keyboard key that corresponds to the underlined letter in a menu entry to 
select it. In both cases, you can leave menu navigation without making a 
selection by pressing the key. 


On an Apple Mac you can access the wxMaxima menus by pressing and 


releasing the key while holding down the and (fm) keys (or just: (ctrl) 


on some keyboards), navigating through the menus using the keyboard 


arrow keys, and then pressing to make your selection. You can 
leave menu navigation without making a selection by pressing the key. 


The wxMaxima toolbar cannot be accessed using the keyboard alone, but 
this is not used extensively in this Guide and alternative typed commands 
for accessing the functions are given as required. 


Changing wxMaxima colours and fonts 


You can change the appearance of the text used by wxMaxima as follows. 
First, click on the following toolbar icon. 


a 


Alternatively, select Configure from the Edit menu (or, on an Apple Mac, 
select Preferences from the wxMaxima menu). 


In the configuration window that opens, select the following ‘Style’ icon. 


i 


94 


You can then change the colour, style, and sometimes the font and size, of 
individual elements of the wxMaxima interface by first choosing the 
appropriate element from the ‘Styles’ list, and then selecting your 
preferred settings. 


Where a selected element has a greyed out font button this is usually 
because it shares its font with other elements of the interface. For 
example, those elements that form part of the mathematical output of the 
interface share a font that you can change by first clicking on the Math 
font ‘choose font’ button and then selecting your preferred font and size. 
Most of the remaining elements with a greyed out font button share a font 
that you can change in a similar way by clicking on the Default font 
‘choose font’ button. 


Unfortunately the Configuration window is not keyboard accessible, and so 
if you only use a keyboard you may need the assistance of a non-medical 
helper to initially configure the system. If you choose a large font size, 
some horizontal scrolling may be necessary to read all the text displayed in 
a worksheet. 


The Configuration window does not change the colours and fonts used in 
graphs produced by Maxima. The ‘Changing graph properties’ subsection 
below explains how to change the colours, fonts and line thickness used in 
graphs. You should read this subsection when you study Subsection 3.3 of 
this Guide, where plotting graphs is introduced. 


Using a screen magnifier or screenreader 


The text in wxMaxima can become pixelated when you use a screen 
magnifier. An alternative is to increase the size of the text font as 
described above and/or use the zooming facilities on the Edit menu, which 
you can also access using Alt-I (zoom in) and Alt-0 (zoom out). 


The wxMaxima interface is not accessible to a screenreader. If you use a 
screenreader you might wish to consider using the command-line interface 
to Maxima instead of the wxMaxima interface. You should ensure your 
screenreader is set to read all the punctuation. Details of using Maxima 
with the command-line interface are given below. The graphs produced by 
Maxima are not accessible to a screenreader. You might need assistance 
when using these. 
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Changing graph properties 
You can change the thickness of curves plotted in graphs by including an 
additional argument in the wxplot2d command, of the form 

(style, [lines, 3]] 


where the number indicates how thick the line should be. By default, lines 
have thickness 1. 


So, to plot a graph of the equation y = 2x with a thicker line, enter, for 
example, 


wxplot2d(2*x, [x, 0, 1], [style, [lines, 3]]); 
Typing this additional argument each time you want to plot a graph can 
be cumbersome, so you might like to assign this argument to a variable 


that you can reuse. For example, to assign your preferred style to the 
variable thick enter 


thick: [style, [lines, 3]]; 
Then you can plot graphs using commands such as 
wxplot2d(2*x, [x, 0, 1], thick); 


To change the colour of a curve, you can include an argument of the form 
[color, red] as described in Subsection 3.3. Alternatively, you can add a 
further element to the styles argument described above. For example, 
including 


(style, [lines, 3, 5]] 


produces curves with thickness 3 and colour number 5, which is black. The 
colour number codes used by the style argument are given in the table 
below. 


style argument colour numbers 


Number 


blue 
red 
green 
magenta 
black 
cyan 


Onkwnye 


Note that the value of any variable you define to help set curve properties 
(such as thick above) will be lost when you close Maxima at the end of a 
session unless you save your session. The style argument also only 
changes properties of the curve plotted, not other aspects of the graph 
such as the axes and labels. 


Changing the axes and labels, and ensuring such changes are automatically 
used each time Maxima is started is more complicated. A computer file is 
provided in the Accessibility section of the OU Maxima website 
learni.open.ac.uk/site/maxima which you can download and configure 
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to set the graph properties to your needs. These settings will then be 
automatically used for each graph in every Maxima session. Details of how 
to download and use this are given in the Accessibility section of the OU 
Maxima website. 


Getting started with Maxima using the 
command-line interface 


This section describes how to use Maxima with the command-line 
interface. If you are using this interface you should read this section 
instead of Section 2 of this Guide, which this replaces. Differences between 
the wxMaxima interface and the command-line interface in later sections 


are minor and are mentioned within those sections when necessary. 


Activities in this section are numbered to match the corresponding activity 
in Section 2. The format of the solutions given at the end of this Guide 
may not exactly match the output you obtain using the command-line 
interface. 


If you encounter unexpected problems when working through this section, 
remember to check the Frequently asked questions (FAQ) section on the 
OU Maxima website: learn1.open.ac.uk/site/maxima. 


Using Maxima 


Computer activity 2 Starting Maxima 


Start the command-line interface to Maxima on your computer and keep it 
open to work with as you study this section. 


If you cannot remember how to do this on your computer, see the 
OU Maxima website: learni.open.ac.uk/site/maxima 


The command-line interface to Maxima is illustrated in Figure 12. 


(MX Command tine Masima 


cct) 
See the file copy: 

Dedicated to the memory of William Schelter 

The function bug report() provides bug reporting information 

(441) 


Figure 12 The initial command-line interface window 


Maxima accessibi 
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At the top of the window is some introductory text, followed by 
(hil) 
This is how Maxima labels lines. The i stands for ‘input’: this is input line 


number 1 and is where you can type your first command. These commands 
either instruct Maxima to perform calculations or change system settings. 


The percentage symbol (%4) is used by Maxima to indicate objects built 
into the system. You will learn more about this later. 


The following activity shows you how to use Maxima to perform a simple 
calculation. Follow the steps using Maxima as you read the activity. 


Follow these steps to use Maxima to calculate 2+ 3 x 3. 


1. Click anywhere on the Maxima window, to ensure that the text you 
type next is directed to the correct place. 


2. Type the following exactly as it appears here: 
2+3%4/5; 


including the semicolon (;) at the end. 


The text you type will appear next to the input line number (%i1), as 
indicated by a flashing horizontal line known as the editing cursor. 


3. Press : 
If you forget to type the semicolon (;) at the end of the command, 


the editing cursor simply moves to the next line, waiting for you to 
type more. Maxima will continue to do this until you indicate you 

have reached the end of the command by typing a semicolon. So if 
you forgot the semicolon, then type it now and press again. 


Pressing (_Enter_) instructs Maxima to calculate the expression. This 
is known as evaluating the expression. 


The result is displayed under the input line, preceded by the label 
(401). The o stands for ‘output’, so this label identifies the line as 
output line 1. It gives the output corresponding to input line 1. 
Throughout each Maxima session the line numbering starts at 1 and 
then increases by one for each new input (and output). 
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Maxima then displays another input line number, ready for you to 
enter your next command. 


BB Command tine Mauma = tee 


Maxima 5.30.0 http: //maxima.sourceforge.net 
jusing Lisp GNU Common Lisp (GCL) GCL 2.6.8 (a.k.a. GCL) 
Distributed under the GNU Public License. See the file COPYING 
Dedicated to the memory of William Schelter 
ithe function bug report() provides bug reporting information 
(842) 24384/5 
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(¥o1) 


(a2) 


You should end all Maxima commands with either a semicolon (;) (as you 
saw in Computer activity 3) or a dollar sign ($). Ending a command with 
a dollar sign tells Maxima to perform the instruction, but not to display 
the answer. 


The result of the calculation in Computer activity 3 was displayed as 


formatted mathematics over several lines, replicating how mathematics is 
usually written. This can, however, cause problems if you are using a 
screenreader. The display behaviour can be changed for every Maxima 
session as described in the following activity. 


Computer activity 5 Configuring the Maxima display of 
mathematics 


You should only do this activity if you wish to change how Maxima 
displays mathematics. 
(a) First exit Maxima by entering the command 
quit (); 
(on some installations you may also have to close the window). 


(b) Turn off the two-dimensional display of mathematics by adding the 
Maxima command display2d:false; to your maxima-init .mac 
configuration file, as follows. 


(i) Follow the instructions in the Accessibility section of the OU 
Maxima website to locate your maxima-init.mac configuration 
file, or create one in the correct location as appropriate. 
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(c) 


(ii) Open your maxima-init.mac file in a plain text editor (such as 


Notepad on a Windows computer, or TextEdit on an Apple Mac), 


and add the text 
display2d:false; 


on a new line. There are further details on how to do this in the 
Accessibility section of the OU Maxima website. 


(iii) Save your modified maxima-init.mac, and close your text editor. 


Restart Maxima and enter the calculation from Computer activity 3 
again, that is, enter 


2+3%4/5; 


Notice that Maxima now displays the result on a single line: 


22/5 


Troubleshooting Maxima 

If you are using Maxima and it does not seem to be responding, try the 
following suggestions. If you encounter other problems, check the 
Frequently asked questions on the OU Maxima website: 
learni.open.ac.uk/site/maxima. 


Troubleshooting Maxima 


If Maxima does not seem to be responding, the following might be 
helpful. 


Check you have ended your previous command with either a 
semicolon (;) or a dollar ($). Maxima will not do anything 
without them! 


Check that Maxima is not waiting for you to type something. Is 
a question being displayed that you have not yet answered? 


If Maxima is taking too long to evaluate a command, you can 
interrupt it by typing Ctrl-C. 


If you do this, you may see the following output, or something 
similar, which you can safely ignore. 


Maxima encountered a Lisp error: 
Console interrupt. 
Automatically continuing. 


To enable the Lisp debugger set *debugger-hook* to nil. 
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Closing wxMaxima 

When you have finished using Maxima, close it by using one of the 
following methods. 

e = =©Enter quitQ; 


e Ona Windows computer, click the usual small cross button at the top 
right-hand corner of the Maxima window. 


Basic calculations using Maxima 


In Computer activity 3 you were asked to enter a numerical expression into 
Maxima using the symbols + for addition, * for multiplication and / for 
division. Maxima evaluated it using the rules for the correct order of 
mathematical operations; that is, the BIDMAS rules. The way in which 
calculations and commands have to be entered in Maxima is known as its 
syntax. 


The Maxima syntax for basic mathematical operations is summarised 
below. 


Basic mathematical operations 


Operation Syntax | Example 
Addition a 443; 
Subtraction = 4-3; 
Multiplication * 4*3; 
Division i 4/3; 
Brackets Cand) | 2*(3+4); 
Powers, for example 2° “ or ** = | 243; 
2**3; 
Square root, for example /5 sqrt (5); 
Magnitude (absolute value) abs (-3) ; 


In this Guide, we will use “ for powers. Note that many Maxima 
commands take the form of a command name, such as sqrt or abs, 
followed by a pair of round brackets containing one or more objects, such 
as numbers, that the command operates on. Such an object is called an 
argument of the command. 


In the summary tables in this Guide, arguments of commands are often 
represented by ff, as illustrated in the table above. 
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Warning: Do not omit multiplication signs! 

A common error when first using Maxima is to omit multiplication 
signs. For example, you might type 

e = =2sqrt(2) rather than 2*sqrt (2), or 

e 2(3+4) rather than 2*(3+4). 


Doing this will result in an error. The displayed error message might 
include one of the following statements: 


e parser: incomplete number; missing exponent? 


e incorrect syntax: Syntax error 
incorrect syntax: Too many )'s 
incorrect syntax: Premature termination of input at ;. 


e ... is not an infix operator 


If you receive one of these error messages, first check for missing 
multiplication signs! 


Warning: Take care with question marks (?) 


Do not type question marks when asking Maxima to evaluate an 
expression. The symbol ? is a special symbol in Maxima which 
provides access to underlying systems. 


The one exception to this is when you are using the commands for 
obtaining help from the system, which you will see later. 


Maxima generally ignores all spaces that you type while entering a 
calculation or command, so you can use spaces to make a command easier 
to read. You can calculate roots other than square roots, such as cube 
roots, by using the index laws. 


Computer activity 6 Calculating cube roots 


(a) Remembering that 7 = 278, calculate 9/27 by entering the 
following line into Maxima. 


27*(1/3) ; 


Remember to press to evaluate the command. 


(b) What happens if you enter 
ZOU Ss 


instead? Can you explain this? 
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Maxima usually tries to calculate quantities exactly rather than 
approximating them by decimal numbers. It does this by manipulating the 
symbols in the expression according to the rules of mathematics. This is 
sometimes known as symbolic computation. You can see an example of 
this in the next activity, which also shows you two different ways to 
instruct Maxima to output a decimal answer instead. 


In Maxima, calculate 1309/3451 in three different ways, by entering each 
of the following expressions in turn. 


(a) 1309/3451; (b) 1309.0/3451; (c) float (1309/3451) ; 


Maxima returns an exact answer for part (a). 


Since part (b) involves a decimal number 1309.0, Maxima returns a 
decimal approximation to the answer. 


In part (c), the command float instructs Maxima to return a 
decimal answer. 


Converting to decimal numbers 


Operation Command | Example 


Convert to a decimal number | float (§) | float (sqrt (2)); 
SS See 


The name of the command float arises from the fact that the method 
used by computers to store decimal numbers internally using the binary 
digits 0 and 1 is called a floating-point representation. The float 
command instructs Maxima to convert a number stored symbolically as a 
mathematical expression to one stored using a floating-point 
representation. 


Maxima usually displays decimal numbers to 16 significant figures (though 
it suppresses trailing zeros at the end of the decimal part of a number; for 
example 0.125 is displayed simply as 0.125 rather than as 
0.1250000000000000). This is also the precision to which Maxima 
performs decimal calculations. You will see later how to change the 
number of significant figures displayed. 
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For each of the following expressions, use Maxima to simplify it and then 
find a decimal approximation for it. 


2/3 
(a) VF (b) (3) (o) 2 


The expression in part (c) doesn’t involve decimal numbers and so the 
result is evaluated exactly. The decimal form of the result 
3.2733906078961419E+150, is given in the command-line interface’s 


way of displaying scientific notation. It means 
3.2733906078961419 x 10!°°. 


Maxima knows the values of standard mathematical constants such as the 
irrational numbers e and 7. It treats them symbolically unless you 
explicitly ask for them in floating point form. 


The constants e and 7 


As with the input and output line labels, the ‘%’ symbol indicates the 
name of a quantity built into Maxima. 


Use Maxima to calculate the value of the expression 7+ ; +e, both 
exactly and as a decimal number. 


Notice that in the exact answer, Maxima displays 7 and e in the form 


they were entered: %pi and Ze. 


Reusing and editing commands 


In Computer activity 9, you may have entered the expression twice, once 
to evaluate it exactly and once for a decimal answer. This is, however, not 
necessary. When you have entered and evaluated an expression, you can 
edit it and then evaluate it again, as shown in the following activity. 
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(a) Calculate 2/3 + 5/27 by entering the following line into Maxima. 
2*sqrt (3)+5*sqrt (27) ; 


(b) Edit your entry in part (a) to calculate 2/3 — 5/27 as follows. 
1. Press the keyboard up arrow key once, so that the previously 
entered command appears against the current input line number 
(on an Apple Mac you may need to hold down the (Alt) key while 
pressing the keyboard up arrow key once). 


You can access all the commands you have previously entered in your 


Maxima session using the up and down keyboard arrow keys. 


2. Edit the expression by moving the editing cursor using the left 
and right keyboard arrow keys, typing new text and using the 
(backspace) or key to make deletions. 


3. Press (Enter). 
SS ES ee 


Maxima also permits the result of one calculation to be used in another 
calculation. The symbol % represents the result of the last evaluated 
command. So, for example, if you enter sqrt (%) ; Maxima calculates the 
square root of the last evaluated output. You can also use the output (or 
input) of other cells, by typing the cell line label in a calculation. For 
example, if you enter 405*3; Maxima calculates the cube of the expression 
in output line 5. 


(a) Enter 6*2; into Maxima and evaluate it. 
(b) On the following input line type 3+%; and enter it. 


Maxima returns 15, the sum of 3 and the previous answer. 


(c) Enter Zoff+5; with §§ replaced by the output line number 
corresponding to the result of 6*2 entered earlier. 


This command adds 5 to the result of the line referenced. 
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Assigning values to variables 


Maxima allows you to assign a value to a variable, and then use the 
variable in further calculations. 


For example, you can assign the value 23 to the variable a, and then 
evaluate an expression in a, such as a” — 3a+ 2. You can also assign 
expressions, both numerical and (as you will see later) algebraic, and 
indeed any other Maxima object, to a variable, and use that variable in 
subsequent commands. 


To assign a value to a variable you use a colon (:). For example, the 
command 


a233 


assigns the value 23 to the variable a. The colon can be read as ‘is 
assigned the value’. 


Variable names can be any combination of letters and numbers that begin 
with a letter. For example, Maxima will accept any of the following 
variable names. 


a A solution solution2 x2b 


Variable names are case-sensitive; for example, the variable x is different 
from the variable X. 


Once a variable has been assigned a value, Maxima will remember this 
value for the rest of your session. This may sometimes surprise you. It is 
easy to forget that earlier in the day you assigned a value to, say, x; you 
may then be confused later on when a calculation involving x gives an 
unexpected result. You can, however, tell Maxima to ‘forget’ about a 
variable to which you have previously assigned a value. You can also ask it 
to list all the variables that you have assigned values to. These commands 
are demonstrated in the following activity. 
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(a) Assign the value of V8 to the variable a by entering a: _sqrt(8); 


The value of a is displayed as output, in an equivalent form. | 
L 


(b) Enter a; 


This displays the current value of a. 


(c) Assign the value of a2 to b. 


Don’t forget to include a multiplication sign when typing ay2. 
The value of b is simplified when displayed. 


(d) Edit the line where you assigned 8 to a, to assign the value V7 to a. 
(e) Enter b; to display the value of b. 


Notice that the value of b has not changed. The variable b was 
defined when a had its original value. 


(f) Enter values; to list the names of all the variables that are currently 
assigned values. 


The variable names are given within square brackets, separated by 
commas. This is how Maxima displays lists. 


g) Enter kill(a); to remove the variable a from the system. 


You should obtain the output done 


h) Enter a; to display the value of a. 


The variable name is output, as a no longer has a value. 


i) Enter b; to display the value of b. 


| b is still defined. | 
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Notice that in Computer activity 12(e), when the value of a was changed, 
the value of b was not affected, since it was defined when a had its original 
value. If you want to update the value of b using a different value of a, 
then you should enter a new value of a then use the up and down 
keyboard arrow keys to obtain the line defining b and then press 
to re-evaluate it. 


Working with variables 


Operation Command 


Assign a value to a variable 
Display the value of a variable 
List all user assigned variables 
Remove an assigned variable 

Remove all assigned variables 


values values; 


kill (WSETEBIE) | kill (a); 
kill(al1) kill (all); 


Note: here the placeholder Wa##ab) represents any variable name. 


You can practise using variables in the following activity. 


(a) Define the variable a to have the value 42. 
(b) Define the variable b to be equal to a® — 3a + 4. 


(c) Define the variable c to be the decimal approximation of the square 
root of b. 


System variables 


There are some variables built into Maxima whose values affect the 
behaviour of the system. These are called system variables. You met 
one system variable in Computer activity 12; the variable values holds a 
list of the names of all the variables you have defined. 


Another system variable, fpprintprec, specifies the number of significant 
figures of decimal numbers that are displayed. The name fpprintprec is 
an abbreviation of ‘floating point print precision’. 


To change the system behaviour so that, for example, only 4 significant 
figures are displayed, use the command fpprintprec:4; to assign the 
value 4 to the variable fpprintprec. 


You can set the variable fpprintprec to any value between 2 and 16. 
Also, setting it to 0 restores the default behaviour of displaying 16 
significant figures. The value of fpprintprec sets the requested number of 
significant figures, but Maxima does not always exactly meet the request. 
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Display the decimal approximation of x to 8 significant figures by doing 
the following. * 


(a) Set the system variable fpprintprec to be 8. 
(b) Display the decimal approximation of 7. 


You can reset the values of fpprintprec and all other system variables to 
their original values by using the reset() command. Try entering 
reset(); followed by float (%pi) ; and check that the value of 7 is 
displayed to 16 significant figures. 


Resetting system variables 


Reset all system variables 


Another useful system variable is leftjust. Setting this to be equal to 
true using 


leftjust:true; 


makes Maxima display the result of a command left-justified against the 
output line number. The default value of leftjust is false, which 
displays output centred in the line. 


Note that the list of variables displayed when typing values; does not 
include any system variables. If one of your variables does not appear in 
the list displayed by values;, then you have probably used a variable 
name that is also the name of a system variable and hence changed the 
value of that system variable. In extreme circumstances this may change 
the behaviour of Maxima, which can be restored by resetting all system 
variables as described above. 
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Saving and printing your work 


After working with Maxima you will probably want to save your 
calculations so that you can reuse or read them at a later date. 


When you are using the command-line interface you cannot save the 
session directly, as you do when saving other files. Instead there are two 
stages in saving your work: 


e If you want to keep a record of the contents of the command-line 
window, then you need to save a transcript of the session. 


e If you want to save the values of all the variables, so that you can 
continue working with them in a new Maxima session, then you will 
need to save the current state of the system. Note that this does not 
give you access to the commands used to define the variables, so you 
cannot edit them in the new session. 


These two processes are illustrated in Computer activity 16. 


When saving your work, you need to tell Maxima exactly where to save 
the file containing the session state or your transcript by giving the 
complete pathname of the file, which includes not only the name of the file 
but also the hierarchy of folders in which the file is located. Unless you do 
this, Maxima will try to save the file in the location in which the system is 
installed on your computer. Not only is this not recommended, but as a 
computer user, you may not have sufficient permission to do this. 


The pathname of a file depends on its location, the type of computer you 
are using, and your computer username. In the following activity, we 
assume that a user with the computer username ‘bill’ wishes to save a 
transcript of his Maxima session in a file called transcript .txt within a 
folder called maxima which he created in Computer activity 1. The maxima 
folder is within Bill’s ‘Documents’ area on a recent Windows computer 
(that is, Windows vista, 7 or 8). In this case, the full pathname of his file is 


C:/Users/bill/Documents/maxima/transcript.txt 


(Notice that Maxima separates parts of the pathname with a forwards 
slash (/) rather than the usual Windows backwards slash (\).) 


If Bill were using a older Windows computer (such as Windows XP), the 
equivalent pathname would be 


C:/Documents and Settings/bill/Documents/maxima/transcript.txt 
or using an Apple Mac the equivalent is 


/Users/bill1/Documents/maxima/transcript.txt 
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(a) Create a transcript of your current session by following these steps. 
(i) Enter 7 
writefile("C:/Users/bill/Documents/maxima/transcript.txt") ; 


replacing the pathname with one appropriate to you and your 
system. 


The path of file is given in double-quote marks. 


This tells Maxima to save a copy of all the subsequent input and 
outline lines to the file named transcript .txt in your chosen folder. 


The output 
Starts dribbling to 
C:/Users/bill/Documents/maxima/transcript .txt 


is given, where the pathname is the one you gave in the command, 
followed by the current date and time (in a probably incorrect 
timezone!). Some computer implementations may produce a different 
output. 


(ii) Enter 
playback () ; 


This command redisplays each of the input and output lines of your 
session, and when used in conjunction with the writefile command 
results in a copy of these lines being stored in your chosen file. 


To playback only a selection of the lines from your session, you can 
use a command of the form 


playback ( (iS Lne) URGE , Last Ant mmber) ); 


which plays back only those lines within the given range. 


(iii) Enter 
closefile(); 


This stops the recording of the session into the file. 


The output Finished dribbling to 
C:/Users/bill/Documents/maxima/transcript .txt is given. 
Again, some computer implementations may produce a different 
output. 
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(iv) Find the file transcript.txt on your computer and open it 
using a text editor, for example, Microsoft Word. 


The file contains a copy of the input and output lines from your 
session. (When opened in some text editors the file may appear all on 
one long line.) 


(b) Save your current session then reload it, by following these steps. 
(i) Enter 
save("C:/Users/bill/Documents/maxima/mysession.txt", all); 
replacing the pathname with one appropriate to you and your 
system. 


This tells Maxima to save your current session to the file 
mysession.txt in your chosen folder. 


The path of the file is output. 


The file is stored using a form that Maxima can reprocess. 


(ii) Close Maxima by entering 
quit(); 
(iii) Restart Maxima in the usual way. 
(iv) Load your previously saved session by entering 
load("C: /Users/bill/Documents/maxima/mysession.txt") ; 


using the same pathname you used in part (b)(i). 


This tells Maxima to restore the state of the system to that saved in 
the specified file. 


The error 


assignment: cannot assign to gf_data(characteristic, 
exponent, reduction, primitive, cardinality, order, 
factors_of_order) 


-- an error. To debug this try: debugmode(true) ; 
might be given, which can be ignored. 
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(v) Enter 
values; 


to display the names of all known variables. 


You should see that all the variables from your previous session are 
known. You might like to check their values are as expected too. 


An alternative way of using the writefile and closefile commands is to 
enter the first at the very start of your session and the second at the end, 
so that everything you do is recorded. 


Saving your work 


Operation 


Start writing a transcript writefile("#ilelpath") ; 
Stop writing a transcript closefile(); 
Playback all the input and playback() ; 


output in a session 

Playback a range of input and | playback ( [Btarelane"Himber, Sndaneynimber) ) ; 
output 

Save the state of a session save ("(filepath", all); 

Reload a session load ("fale path) 


You can also obtain an image of the Maxima window by taking a 
‘screenshot’. To do this on a Microsoft Windows computer, first click on 
the Maxima window, then press Alt-PrintScreen (that is, hold down the 
keyboard key while pressing the PrintScreen key, which may be 
labelled * LES something similar). This stores the image of 


the Maxima window, which you can then paste into a suitable application. 


If you are using Windows Vista, Windows 7 or Windows 8, you might like 
to use the ‘Snipping Tool’, available by typing ‘snip’ into the Start menu. 


To take a screenshot on an Apple Mac computer, press Cmd-Ctr1-Shift-4 
(that is, hold down the Grae)» (ctrl) and keyboard keys while pressing 
the key). Next press the space bar, move the camera pointer over the 
Maxima window and click. This stores the image of the Maxima window, 
which you can then paste into a suitable application. You may also like to 


use the Grab application from the Utilities folder. 
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Getting help 


Maxima includes commands that you can use to obtain help. These are 
demonstrated in the following activity and summarised in the table that 
follows. 


(a) Enter ? float; to show the help for the float command. 


Note the space between the ? and float, which is needed. This is one 
of the few times in Maxima when spaces matter! 


The help information for the command may include more detail than 
you need! 


After giving the help information, Maxima displays the output true. 
This means that the command you entered was successful. 


(b) Suppose that you cannot remember the Maxima command for square 
root, but you do remember that it was something like sqr. 


Type ?? sqr; to list all the help information titles that contain the 
letter sequence ‘sqr’. 


A list of possible titles is displayed. 


Title number 1 is what you were looking for, so enter 1 to display the 
relevant information. 


If more than one title is displayed in response to a ?? command, then 
Maxima will not continue until you enter one of the following: 


the appropriate number, to display your chosen information 


several numbers, separated by spaces, to see help information 
on several topics 


all to display the information on all topics listed, or 


none to see no information. 
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Note 1: the space after the ?, which is needed. 


Note 2: the second pair of commands displays a list of the titles of all 
help information pages whose title contains the given text. To obtain 
help on a particular topic listed, type the number corresponding to the 
required title in the list. 


There is further help with Maxima in the Frequently asked questions 
section of the OU Maxima website, at learni.open.ac.uk/site/maxima. 


If you are reading this section as an alternative to Section 2, you should 
now continue with Section 3. Note, however, the information below 
regarding the plotting of graphs, which you will need when studying 
Subsection 3.3. 


Plotting graphs using the command-line interface 


When using the command-line interface, you should use the plot2d 
command rather than the wxplot2d command described in Subsection 3.3. 
These commands work similarly, and have the same arguments. The main 
difference is that plot2d opens a new computer window containing the 
graph. 


On some computer systems, the plot2d command opens two windows, one 
containing the graph and a second containing a command-line interface to 
gnuplot, which is the system used by Maxima to plot graphs. To prevent 
this second window opening each time, enter the following command. 


gnuplot_viewargs: "-persist ~s"; 

The Accessibility section of the OU Maxima website contains details of 
how to configure Maxima to automatically prevent the display of this 
second window. 


The graph window can be closed by pressing the usual cross at the top 
right-hand corner, or by pressing Alt-F4. You will not be able to enter 
any commands into the Maxima command-line interface until the graphing 
window is closed. 


The graph can be printed using the print button on the graph window, or 
saved as an enhanced metafile (.emf) file using the save button. These files 
can be opened with many image viewing programs and can be inserted 
into word-processed documents. 


Operation Command | Example 
Get help on a command ? (Gommand ? float; 
or, describe ((Gemmand)) describe (float) ; 
Find help information whose | 7? J ?? flo; 
title contains the given text | or, describe(J, inexact) | describe(flo, inexact); 
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Computer methods for CAS Activities 
in Books A-D 


Computer solution to Unit 5, Activity 17 


First, load the implicit_plot package and set the default plotting options 

to use equal scales on the «- and y-axes. 

(%il) load(implicit_plot) ; 

(4o1) C: /PROGRA~1/MAXIMA~1 .0/share/maxima/5.30.0/share/contrib/implicit_plot.lisp 


(%i2) set_plot_option([gnuplot_preamble, "set size ratio -1"])$ 
P. plot_p: 


This line is ended with $ to prevent the display of the output of the 
command. 


(a) Assign the equation representing the parabola and the equation 
representing the circle to the variables p and c respectively. 
(4i3) pry=2* x*2-10*x+14; 

(403) y=2x?-10 x+14 
(4i4) c:x2-6¥xty*2-8*y+22=0; 
(404) y?-8 y+x?-6x+22=0 


Now plot the two curves. Appropriate axis ranges are chosen to ensure 
the whole of the circle is shown. 


(4i5) wximplicit_plot({p, c], [x,0,7], [y,0,7]); 


T T i 
= 2°x42-10"x+14. —— 


7 r = 
yreelyexr2 6x22 =0 Sarre) 


(nb) ee | 


(405) 
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(b) To find the points of intersection of the parabola and circle, we solve 
the equations representing each simultaneously, for 2 and y. The 
equations have previously been assigned to the variables p and c 
respectively. 


(416) solve([p,c], [x,y]); 

(406) [[x=1.304970075461879, y=4.356193142057383] , 
[x=1.664051355206847 , y=2.897620365246265] , 
[x=3 . 121372031662269 , y=2.272207016670814] , 
[x=3 .909606299212598,, y=5 .473980309423348] ] 


So, to 2 d.p., the four points of intersection are 
(1.30, 4.36), (1.66, 2.90), (3.12, 2.27) and (3.91, 5.47). 


(c) The maximum number of intersections of a circle and a parabola is 
four. This fact can be justified by the geometric argument that each 
‘half’ of the parabola (each part on one side of the vertex) can meet 
the circle at most twice. 


Alternatively, if the equation of the parabola is used to substitute for 
y in the equation of the circle, then an equation for x containing terms 
in «4, 2°, 2, 2 and a constant term is obtained, and such an equation 
can have at most four solutions (by a result that you will meet in 
Unit 12). 


Computer solution to Unit 7, Activity 23 
(a) The function f is defined as follows. 


(hit) £(x) += (sqrt (14+x%2) ) /8+(2-x)/16; 
2 = 
(ot) £G) =" tee io 


16” 


(b) The graph of f(a) is plotted below. 
(4i2) wxplot2d(f (x), [x,0,2]); 


0.28 : + . - 
© 0.275 + A. 
= 027+ / 
a 0.265 ye 
gS 0.26} A 

‘ S 0.255 + 

(42) 0.25 & me 
0.245 5 
= 024+ \ 2 
3 0285} 9 

0.23 . 1 . 
0 05 1 1.5 2 
x 
(ho2) 
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(c) The value of x that gives the minimum value of f(x) seems to be 
about 0.6. 


(d) The derivative of f can be found and assigned to df (x) as follows. 
(%i3) af (x) :='' (diff (f(x) ,x)); 


5 poe 
(03) df (x): WER oC 


eee ee 
— 8¥224+1 16) 
(e) First, we try to use the solve command to find the solution of 
df (x)=0. 
(%i4) solve (df (x)=0) ; 
ee 
(hoa) (x=V%**y 


So the derivative of f is f'(x 


Maxima returns a rearranged version of the equation, but does not 
solve it. The solve command cannot find an exact solution of 
equation df (x)=0, so we look for an approximation to the solution 
instead. 


To do this we need to find an interval containing the solution. We 
know from the graph of f(z) that its gradient is negative to the left of 
the minimum, and in particular at x = 0, so df(0)< 0. Also, the 
gradient of f(x) is positive to the right of the minimum, in particular 
at 2 = 1, so df(1)> 0. So a solution of df(x)=0 must lie in the 
interval (0,1). An alternative way to find this interval is to plot a 
graph of df (x). 


(4i5) find_root (df (x) ,x,0,1); 
(405) 0.57735026918963 
The solution of the equation f/(x) = 0 is = 0.577 (to 3 s.f.). 
(f) The man should join the road 0.58 km (to the nearest 10 metres) along 
from the point that is closest to his initial position. 


(Because the graph of f indicates that the minimum value of f(x) occurs 
at a local minimum rather than at one of the endpoints of the interval 
(0, 2], in the solution above we don’t consider the values of f(0) and f(2).) 
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Computer solution to Unit 8, Activity 48 
(a) The domain of f is the interval [—1, 00). 
(b) The graph of f can be plotted as follows. 


(hil) £(x) :=(x%2-10*x+15) / (1+sqrt (x+1)) ; 
x?-10 x+15 
Rol) £ :=—— ; 
(hot) £(x) REST 
(%i2) wxplot2d(£(x), [x,-1,15]); 


30 


(42) 


(x42-10"x+15)/(sqrt(x+1)+1) 
a 


(402) 


The z-intercepts of f can be found by solving f(x) = 0. 
First, solve the equation, assigning the result to the variable solns. 


(4i3) solns:solve(f(x)=0) ; 

(ho3) [x=5-V10, x=V/10+5] 

Next, extract the two solutions from the list and assign them to the 
variables soli and sol2 respectively. 


(%i4) soli:rhs(solns[1]); 

(404) 5-10 

(4iS) sol2:rhs(solns[2]); 

(405) V/10+5 

Finally, express the solutions as decimal approximations. 
(%i6) float (sol1) ; 

(406) 1.837722339831621 

(4i7) float (sol2) ; 

(407) 8.16227766016838 

So the two x-intercepts of f are 1.84 and 8.16 (to 3 s.f.). 


119 


Computer methods for CAS Activities in Books A-D 


(d) The area between the graph of f and the z-axis, between the two 
x-intercepts is 


i 22 —102-+15 . 
1sa7.. l+vz+1 ~" 
We shall concentrate on evaluating 


Bo a? — 102 +15 ie 
1as7.. l+Ve+l ” 


that is, 
8.162... 
[foe 
1.837... 
remembering the required answer is the negative of this integral. 


The integrate command cannot evaluate this definite integral, so we 
find an approximate value instead. 


(4i8) quad_qags (f(x) ,x,sol1,s012) ; 
(408) [-12.37959499949233, 1.3744111403970839 107!9 , 21,0] 


So 


8.162... 
f f(x) da = —12.4 (to 3 s.f.). 
1.837... 


Remembering the required area is the negative of this, the required 
area is 12.4 (to 3 s.f.). 


Computer solution to Unit 8, Activity 49 


The area of the cross-section, in square centimetres, is 
15 
3 Va + 1dz. 


The integrate command cannot evaluate this definite integral, so we find 
an approximate value instead. First, assign the integrand to the variable y. 


(hil) y: (1/3) *sqrt (x43+1) ; 


vxe+1 
3 


(hot) 


Then find an approximate value for the definite integral. 


(%i2) quad_qags(y,x,-1,1.5); 
(%o2) [0.93924271889094 ,7.1021410974481114 107°, 231,0] 


The value of the definite integral is 0.939 (to 3 s.f.). So the volume of 
plastic required to make one metre of the edging is approximately 93.9 cm’. 
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Computer solution to Unit 8, Activity 50 


The change in the displacement of the object from time t = 0 to t = 10 is 


10 2 
sin {| — } dt. 
h(i) 
The integrate command calculates this definite integral in terms of 


values of the Error function, erf. Instead, we shall find an approximation 
to its value using the quad_qags command. 


(%i1) quad_gags(sin(t*2/150) ,t,0,10); 

(401) [2.15266541148388 , 2.3899387041435914 10°**, 21,0] 

So the value of the definite integral is 2.15 (to 2 d.p.). 

The displacement of the object at time t = 10 is (6 + 2.15)m = 8.15m, to 
the nearest centimetre. 

Computer solution to Unit 9, Activity 28 


(a) The matrix form of the equations is 


ere) 


Assign the coefficient matrix to the variable A, and the vector formed 
by the right-hand sides of the equations to b. 


(Ail) A:matrix([2,2,3],[2,3,4],[4,7,10]); 
2 2 3 

(401) }2 3 4 
4 7 10 


(412) b:matrix([2] , [-1],[3]); 
2 

(ho2) |-1 
3 


The inverse of A can be found as follows. 


(413) invert (A) a 


sie ie! 
Ce age 
oS) |=2: <@) “=n 
a ae | 
So we have 


Do sNee 
23 A = le 
47 10 


RNR 
| 
2 Rie 
| 
ero 
or 
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Therefore 


Ga 


This can be calculated using 


(Hid) invert (A) .b; 
0 

(%o4) |-11 
8 


tole 


The solution is z = 0, y = —11, z = 8. 


Note that an alternative way to solve these equations using Maxima is 
to use the solve command, as follows. 


(435) solve ((2*x+2*y+3*z=2, Qex+3ey+4ez=-1, 4x+7*y+10*z=3], [x,y,z]); 
(%05) ([[x=0,y=-11,2z=8]] 


The matrix form of the equations is 


(0-0) 


Assign the coefficient matrix to the variable A, and the vector formed 
by the right-hand sides of the equations to b. 


(416) A:matrix([3,2,1],[4,3,1],[7,5,1]); 
bate Be! 

(406) |4 3 1 

7 5 

(4i7) b:matrix([1] , [2], [4]); 

1 

(ho7) |2 

z 


The inverse of A can be found as follows. 


(/i8) invert (A) ; 


ay eS. oh 
(408) |-3 4 =1 
L 2H i er | 

So we have 
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Therefore 


Ceres: 


This can be calculated using 
(4i9) invert (A) .b; 

| 
(%o9) | 4 

2 


The solution is x = —3, y= 4, z= 2. 


Note that an alternative way to solve these equations using Maxima is 
to use the solve command, as follows. 


(%i10) solve([3*x+2*y+z=1, 4*x+3ey+z=2, 7#x+5+*y+z=1], [x,y,z]); 
(4010) C[x=-3,y=4,z=2]] 


Computer solution to Unit 9, Activity 29 


(a) The matrix form of the equations is 


ala, 


12:8 

The coefficient matrix is (: 3 ) . Its determinant can be found 
aa gag 

using Maxima as follows. 


(4i1) determinant (matrix([1,2,3],[2,3,1],[4,7,7])); 
(hol) 0 


The determinant of the coefficient matrix is zero, so the coefficient 
matrix is non-invertible and the equations do not have a unique 
solution. 


The matrix form of the equations is 
2.2 3\ (xz 5 
23 1) /y)=([6). 
477 z a 


2 
The coefficient matrix is | 2 
4 


os cd 
“1H 0 


) . Its determinant can be found 


using Maxima as follows. 


(%i2) determinant (matrix([2,2,3],[2,3,1],[4,7,7])); 
(702) 14 


The determinant of the coefficient matrix is 14. This is non-zero, so 
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the coefficient matrix is invertible and the equations have a unique 
solution. 


(c) The matrix form of the equations is 


1 Lot 1 = 0 
res Tee lggena| | |! e5)) al |e 
2 0 3-1 7) (imma (a! 
6 2 @ 0 z 0 
1 1 -1 1 
. ae t =f BE 1 
The coefficient matrix is Bip) gt at 
G2 ae 


Its determinant can be found using Maxima as follows. 


(%i3) determinant (matrix([1,1,-1,1],[1,-1,1,1],[2,0,3,-1],[6,2,4,0])); 
(403) 0 
The determinant of the coefficient matrix is zero, so the coefficient 


matrix is non-invertible and the equations do not have a unique 
solution. 


Computer solution to Unit 10, Activity 20 
(8) p= —$(1.2)"" (=1,2,3,...): 


The first 30 points on the graph of the sequence (z,,) can be plotted as 
follows. 


First, define the sequence. 
(fit) x[n] :=-(1/3)*(1.2)%n; 
(hot) m:=(-3) cE 


Then generate the points to be plotted and plot them. 


(4i2) xpts:makelist ({n,x[n]] ,n,1,30)$ 
(4i3) wxplot2d([discrete,xpts], [style,points], [xlabel,"n"], [ylabel,"x[n]"]); 


x{n] 


(4t3) 


(403) 


The sequence (a,,) is decreasing and x, —+ —0o as n — 00. 
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(b) yn =5(-0.9)" (n =1,2,3,...). 


The first 30 points on the graph of the sequence (y,,) can be plotted as 
follows. 


First, define the sequence. 


(414) y[n] :=5*(-0.9)n; 
(404) yn:=5(-0.9)" 


Then generate the points to be plotted and plot them. 


(%i5) ypts:makelist([n,y[n]] ,n,1,30)$ 
(4i6) wxplot2d([discrete,ypts], [style,points], [xlabel,"n"], [ylabel,"y[n]"]); 


5 r r r T T 
4-e 4 
3 2 4 
2 Cine 4 
1 Se ees 5 
cS oe © 0.9 
cues) = ae aie eee 

e 

2 e 4 
3 on 4 
-4 5 4 

pei - n ——— oe 


o 
oa 
3 
a 
8 
& 
8 


(406) 


The sequence (y,,) alternates in sign and y, — 0 as n — oo. 

(c) 2 =—7(-2.5)" (n= 1,2,3,...). 
The first 15 points on the graph of the sequence (z,) can be plotted as 
follows. 


First, define the sequence. 

(4i7) z[n] :=- (1/10) *(-2.5)%n; 
1 n 

(ho7) zi=(-5) (-2.5) 
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Then generate the points to be plotted and plot them. 
(418) zpts:makelist([n,z[n]],n,1,15)$ 
(%i9) wxplot2d([discrete,zpts], [style,points], [xlabel,"n"], [ylabel,"z[n]"]); 
100000 T T T T T T T 
80000 
60000 
40000 
20000 
0 
-20000 
-40000 


z{[n] 


Cht9) 


(Ho09) 


The sequence (z,) alternates in sign and is unbounded. 


Computer solution to Unit 10, Activity 21 
(a) an =17-—3(1.2)" (n=1,2,3,...). 


The first 30 points on the graph of the sequence (a,,) can be plotted as 
follows. 


First, define the sequence. 

(441) aly) Ts 

(hol) ag Balto: Reon 

Then generate the points to be plotted and plot them. 


(412) apts:makelist([n,a[n]] ,n,1,30)$ 
(4i3) wxplot2d([discrete,apts], [style,points], [xlabel,"n"], [ylabel,"a[n]"]); 


Ceecdccccsos 


cc 
(43) % -30 


(h03) 


The sequence (a,,) is decreasing and a, —+ —co as n + oo. 
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(b) bn = 5(—0.9)" +45 (n =1,2,3,...). 


The first 30 points on the graph of the sequence (b,,) can be plotted as 
follows. 


First, define the sequence. 


(4i4) bn] :=5*(-0.9)n+45; 
(h04) by:=5(-0.9)"+45 


Then generate the points to be plotted and plot them. 


(%i5) bpts:makelist([n,b[n]] ,n,1,30)$ 
(4i6) wxplot2d([discrete,bpts], [style,points], [xlabel,"n"], [ylabel,"b[n]"]); 


50 
49 
48 
47 
46 
45 
(ht6) 44 
43 
42 
4 
40 


bin] 


(%06) 


The terms of (b,,) alternate either side of 45 and b,, + 45 as n + oo. 
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Solution to Computer activity 3 
(hit) BtSA75; 
hot) — 
(hot) = au 
So2+3x 5=%. 


Solution to Computer activity 6 


(a) (4i1) 27%(1/3); 
(fo1) 3 


(b) (4i2) 27%1/3; 
(402) 9 
Without any brackets, 27*1/3 is interpreted according to the 
BIDMAS rules, as (27')/3 which is equal to 9. 
Do not be surprised if your line numbering differs from ours. Each of 
our solutions is prepared as a separate Maxima session with the line 
numbering starting at 1 (unless it is a continuation of a previous 
activity); your numbering will continue to increase for the duration of 


your session. 


Solution to Computer activity 7 
(a) (%i1) 1309/3451; 

Oeeue 

(hot) on 

So 1309/3451 = 
(b) (4i2) 1309.0/3451; 

(402) 0.37931034482759 

So 1309.0/3451 ~ 0.37931034482759. 
(c) (4i3) float (1309/3451) ; 

(%03) 0.37931034482759 

So 1309/3451 ~ 0.37931034482759. 
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Solution to Computer activity 8 
(a) (4i1) sqrt (147); 

(hot) TV3 

(%i2) float (sqrt (147)); 

(402) 12.12435565298214 

So /147 = 7/3 © 12.12435565298214. 
(b) (%i3) (8/5)*(2/3) ; 


4 
(ho8) 


(hi4) float ((8/5)*(2/3)) ; 
(f04) 1.367980757341358 


8 2/3 vi 
So (3) = aya © 1.367980757341358. 


(c) (4i5) 2500; 
(405) 327339060789614187001318969682[91 digits] 217256545885393053328527589376 


(%i6) float (2500) ; 
(406) 3.2733906078961419 10'5° 
So 25° = 3,2733906078961419 x 10!°°. 


Solution to Computer activity 9 
(hil) tpi + %pi/3 + %e%2; 

hot) ne? 

(4i2) float (%pi + %pi/3 + %e%2); 
(402) 11.57784630371704 


4 
Somnte ; += > +e? = 11.57784630371704. 


Solution to Computer activity 10 
(a) (411) 2*sqrt(3)+5*sqrt (27) ; 
(hol) 17/3 
So 2V3 + 5V27 = 17V3. 
(b) (412) 2*sqrt(3)-5*sqrt (27) ; 
(402) -13/3 
So 23 — 5/27 = -13V3. 
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Solution to Computer activity 11 


At the end of this activity, your worksheet should look like the following, 
although line numbers may differ. 

(hid) 6*2; 

(404) 12 

(4i5) 3+%; 

(%05) 15 

(hi6) %04+5; 

(406) 17 


Solution to Computer activity 12 


(a) (Ail) a: sqrt (8); 
(hor) 29/2 


(b) (442) a; 
(Ho2) 29/2 


(c) (4i3) b: axsqrt(2); 
(403) 4 


(d) (414) a: sqrt(7); 
(404) V7 


(e) (hi5) b; 
(405) 4 


(f) (4i6) values; 
(406) [a,b] 


(g) (hi7) kill (a); 
(h07) done 


(h) (418) a; 
(408) a 


(i) (hid) b; 
(ho9) 4 
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Solution to Computer activity 13 
(hil) a:42; 

(hoi) 42 

(%i2) b:a*3-3+at4; 

(402) 73966 

(413) c:float(sqrt(b)); 

(%03) 271.9669097519035 


Solution to Computer activity 14 


(fil) £pprintprec:8; 

(hoi) 8 

(This sets the number of significant figures to be shown in the display of 
decimal numbers.) 


(4i2) float (%pi) ; 
(402) 3.1415927 
So the value of 7 is 3.1415927 (to 8 s.f.). 


Solution to Computer activity 15 


The circumference of a circle of radius 5 is 31.42 (to 2 d.p.) and the area is 
78.54 (to 2 d.p.). 


Solution to Computer activity 17 


(a) 


(Ail) ? float; 

-- Function: float (<expr>) 

Converts integers, rational numbers and bigfloats in <expr> to 
floating point numbers. It is also an ‘evflag’, ‘float’ causes 
non-integral rational numbers and bigfloat numbers to be converted 
to floating point. 

There are also some inexact matches for ‘float’. 

Try ‘?? float’ to see them. 


(hol) true 


(4i2) ?? sqr; 

0: isqrt (Functions and Variables for Number Theory) 

1: sqrt (Root Exponential and Logarithmic Functions) 

2: sqrtdenest (Package sqdnst) 

3: sqrtdispflag (Functions and Variables for Display) 

Enter space-separated numbers, ‘all’ or ‘none’: le 

-- Function: sqrt (<x>) 

The square root of <x>. It is represented internally by 
“<x>%(1/2)?. See also ‘rootscontract’. 

‘radexpand’ if ‘true’ will cause nth roots of factors of a product 
which are powers of n to be pulled outside of the radical, e.g. 
‘sqrt (16*x*2)’ will become ‘4*x’ only if ‘radexpand’ is ‘true’. 


(402) true 
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Solution to Computer activity 18 


(a) (hil) expand( (2+x)*(3*x+4)%3 ); 
(ot) 27x*+162x5+360x?+352x+128 
So (2 +2) (3a + 4)9 = 27x + 16223 + 3602? + 352x + 128. 


(b) (4i2) factor( 2#x*3-7*x*2-10*x+24 ); 

(%o2) (x-4) (x+2) (2x-3) 

So 2a° — 7x? — 10x + 24 = (x — 4)(x + 2)(2r — 3). 
(c) (4i8) p: (e341) /(x+1) ; 


x8+1 
3) — 
(ho3) Sy 


(d) (414) fullratsimp(p) ; 
(ho4) x?-x+1 


(e) (%i5) subst (3,x,p); 
(405) 7 


Solution to Computer activity 19 


(a) (Ail) vr: x°2+2*x=4; 
(Hor) x?+2x=4 


(b) (442) lhs(r); 
(%o2) x?+2x 


(c) (413) rhs(r); 
(403) 4 


Solution to Computer activity 20 


(a) (Ail) solve( x2+2+x=1 ); 
(fol) [x=-V2-1, x=V2-1] 
The solutions of the equation are + = -V2-land r= v2-1. 


(b) (412) solve( 4*x*2-12*x+9=0 ); 
(02) xe5] 


The solution of the equation is x = 3. 
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(c) (4i3) solve( 2*x*3+x%2-5*x+2=0 ); 
a: 
(03) ies » x=-2, x=1] 
The solutions of the equation are z = 3, @=-—2andz=1. 


(d) (414) solve( x“9+2*x-4=0 ); 
(404) [0=x°+2x-4] 
Maxima was unable to solve this equation. 
(e) (4i5) solve( x%2-4*x+5=0 ); 
(405) [x=2-%i, x=%i+2] 
The two complex solutions are 2 = 2—i and x = 2+i. There are no 
real solutions. 


Solution to Computer activity 21 


(a) (fit) A: [2,3,5,7,11]; 
(401) [2,3,5,7,11] 


(b) (442) AL4]; 
(h02) 7 


So the fourth element of A is 4. 
(c) (i) (%i3) solns:solve(x*2-x-1=0) ; 


5-1 S+1 
(403) [xe YE , wv ] 
So the solutions are x = — = i and 2 = - ve 


(ii) (%i4) float (solns) ; 
(404) [x=-0.61803398874989, x=1.618033988749895] 
So decimal approximations to the solutions are 
«x = —0.61803398874989 and x = 1.618033988749895. 
(iii) (415) rhs(solns[1]); 
v5-1 


(Ho8) a 


Solution to Computer activity 22 


(a) (Ail) eq: 3*y+4=34x%2+9%x; 
(Hol) 3y+4=3x?+9x; 


(b) (442) eq-4; 
(02) 3y=3x?+9x-4; 


(c) (hi3) %/3; 
3x?+9x-4 


(%03) y= 3 
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Solution to Computer activity 23 
(a) (%i1) eqn: c=(2*atb+3*c)/(a-b) ; 
2atbt+3c 
a-b 


(hot) c= 


(b) (%i2) solve(eqn, a); 
_ (b+3) +b } 


(ho2) tax 


(b+3)e+b 


Soa= = 


Solution to Computer activity 24 
(a) (Ail) eql:3*x*2+2*y=20; 
(hol) 2y+3x?=20 
(Wi2) eq2:4*x-3*y=-4; 
(ho2) 4x-3y=-4 


(b) (%i3) solve( [eqi,eq2], [x,y] ); 
(403) [Ex=2, y=4] , [ x=-% , y=-8 1] 
The two solutions are x = 2, y=4 and r= -%, y=- 


Siz 
ae 


Solution to Computer activity 25 


(a) (4i1) wxplot2d( 4*x*2-8*x+2, [x,-1,2]); 
14 eet T 


2 
x 
9 
(itt) 
+ 
x 
(hol) 
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(b) (4i2) wxplot2d( 4*x*2-8*x+2, [x,-1,2], [y,-5,5]); 
plot2d: some values were clipped. 


4 
x? 
256 . 
(hare h eS , 
= -2 ——— 
-4 
So a a 
LO SEOUN OSS Tt aise ce 
x 
(02) 


(c) 
(4i3) wxplot2d( 4*x*2-8*x+2, [x,-1,2], [y,-5,5], [color, green]); 
plot2d: some values were clipped. 

r T Tr T 


4 4 
:: 2 h, 4 
PO 4 
tsa ues 
& 2 — | 
4 4 
‘ H ‘ . . 
A -0.5 0 0.5 1 1.5 2 
x 
(403) 


Solution to Computer activity 26 
(a) (4i1) wxplot2d( [15*t%3+10*t, 60*t], [t,0,2] ); 
140 


15°t43+10°t —— 
120 60't —/ 
100 
80 
(41) 60 
40 
20 
0 
t 
(hot) 
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(b) 


(c) 


Solutions to Computer activities 


(4i2) wxplot2d( [15*t*3+10*t, 60*t], [t,0,2], [color, green, magenta] ); 
140 T 
120 
100 


“SSO — 
60't 


y. 


(%t2) 60 


(ho2) 


(4i3) wxplot2d( [15*t*3+10*t, 60*t], [t,0,2], [color, green, magenta], 
(legend, "si", "s2"] ); 
140 a T T 
120 
100 
80 


(U3) 60 


(403) 
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(4i4) wxplot2d( [15*t*3+10*t, 60*t], [t,0,2], [color, green, magenta], 


[legend, "si", "s2"], [xlabel, “time, t (h)"], [ylabel, 


140 =—— T 
Ee 120 2] 
g — 
> 100 
= 80 
(ht4) § 60 
4 40 
3% 20 = 
0 —j_ . 1 
0 0.5 1 15 2 
time, t (h) 
(h04) 


Solution to Computer activity 27 


(a) 


(b) 


(c) 


(d) 


(f) 


(hil) £ (x) :=x*5-3¥x%2+4; 
(fot) £ (x) :=x5-3x2+4; 


(%i2) £(2); 

(02) 24 

So f(2) = 24. 

(hi3) g(x) := Cea 
x 

(ho3) g(x) := aoe 

(hi4) g(3); 

(h04) 2 

So g(3) = 2. 

(hi) g(2); 

expt: undefined: 

#0: g(x=2) 

-- an error. 


0 to a negative exponent. 


To debug this try: debugmode(true) ; 


(%i6) h(1); 
(406) h(1) ; 


(4%i7) functions; 
(Ho?) [£(x), g(x)] 


"displacement, s (km)"] ); 
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Solution to Computer activity 28 
(a) (418) wxplot2d(£(x), [x,-2,2]); 


+ 
+ 
s SS 
* 
ts) 
1/05 0° 05 1 LS “2 
x 
(408) 


(b) (4i9) wxplot2d(g(x), [x,-5,5]); 
800 T T T r 
600 ; 4 


200 7 


(U9) 


(x-1)/(x-2) 


3 
S 
Ss 


(h09) 
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(c) (4i10) wxplot2d(g(x), [x,-5,5],[y,-5,5]); 
plot2d: some values were clipped. 


: 1 = . 
4 \ 4 
z 2 es} 
=o ata 4 
(ut10) -& 
“5 | 
-4 | 4 
: : i ki 1 
-4 -2 0 2 4 
x 
(%010) 


Solution to Computer activity 29 


(a) (hil) he*5; 
(hot) he®; 


(b) (4i2) float (%e%S5) ; 
(402) 148.4131591025766 
So e® = 148.41 (to 5 s.f.). 
(c) (4i3) exp(5); 
(h03) he®; 


(d) (4i4) float( log(8) ); 
(04) 2.079441541679836 
So log(8) = 2.079 (to 4 s.f.). 
(e) (4i5) solve( pe es 
log(-V5) log (5) 
pacers mee ; eer 
log(5) 
2log(3)" 


] 


The only real solution is 


Solutions to Computer activities 
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Solution to Computer activity 30 


(a) (hil) u:( log(x*3+2*x*2+x) - log(x) )/log(x+1) ); 
(401) log (x3+2x7+x)-log (x) 


log (x+1) 
(b) (412) radcan(u) ; 
(402) 2 
Notice, 
In(x? + 22? + x) — In(x) 


Ine +1) oo +1) (In(a? + 22? + x) — In(z)) 


oa (in (ce ae = 


(In (a? + 22 + a 


~ ae +1) 


- ne. In(w + 1) (in Gr 1? 


_ 2In(x +1) 
~ In(e +1) 
=2. 
(c) (4i3) v:log(a)+2*log(b)-log(c) ; 
(403) -log(c) + 2log(b) + log(a) 


(d) (4i4) logcontract(v) ; 
a b? 
(404) log “a 


Solution to Computer activity 31 
(a) (4i1) cos (%pi/6) ; 

‘ 3 

(hot) rt 


Cy 
So, cos (=) = co 
(b) (%i2) sin(45*%pi/180) ; 
402) 2 
(%02) Ja 


eee gee 
So, sin 45°= Vi 
(c) (413) sec(3.4); 
(403) -1.034342024711439 
So, sec(3.4) + —1.034342024711439. 
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(d) (414) asin(sqrt (3)/2); 


(h) 


Solution to Computer activity 32 


m Tv 
(h04) = 


So, sin! ¥3 =e 
2 3 


(4i8) csc (%pi/7); 
a qT 
(405) csc (5) 


(416) tan(-%pi/12) ; 
PEC ee Ce 
(408) tan(4) 


(4i7) acos(2.0); 


(407) 1.316957896924817 %i 


So, cos~! 2.0  1.316957896924817i. 


(4i8) acos(2); 
(408) acos(2); 


(a) (hil) solve (cos (x)-0.6=0) ; 


rat: replaced -0.6 by -3/5 = -0.6 
solve: using arc-trig functions to get a solution. 
Some solutions will be lost. 


). 


(hot) [x=acos (2); 


So, a solution is 2 = cos! ( 


(4i2) float (%); 


3 


5 


(402) [x=0.92729521800161] 


So a solution is 2 = 0.927 (to 3 s.f.). 


(4i3) solve(sin(3*x+4)=1/2) ; 
solve: using arc-trig functions to get a solution. 
Some solutions will be lost. 


7-24 
403 a 
(h03) [x is ] 
So a solution is 2 = le 
18 


Solutions to Computer activities 
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(d) (414) solve(cos(2*x+4)-sin(x/3)=0) ; 
(H04) (sin(=)=cos (2x+4)] 


Maxima cannot find an exact solution of this equation, so returns the 
equation slightly rearranged. 
Solution to Computer activity 33 
(a) (i) (Ail) £x) :=cos(2*x+4)-sin(x/3) ; 
(hot) £(x) :=cos (2x+4)-sin( =) 


(ii) (Ai2) wxplot2d(f(x), [x,-2,2]); 
2 


(4t2) 


cos(2*x+4)-sin(x/3) 
Oo 
a 


(02) 


Since £(0)< 0 and f£(1)> 0, there must be a solution in the 
interval (0,1). 

(iii) (413) find_root (f(x), x, 0, 1); 
(%03) 0.42743338823081 
So a solution of cos(2x + 4) — sin(x/3) within (0,1) is 0.427 
(to 3 s.f.). 

(b) (444) g(x) := (sin(x+2))%2-3+x; 
(ho4) g(x) := sin(x+2)?-3x 
(hid) wxplot2d(g(x), [x,-2,2]); 


= 


4 ee 
x ae 
m 2 
z 0 
(yes) 
= 2 
& 
4 
6 a 4 4 4 
245 41 05 0 05 
x 
(405) 
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Since g(0)> 0 and g(1)< 0, there must be a solution in the 
interval (0, 1). 


(If you chose an interval other than (—2,2) then your graph will look 
different from the one given above. If, in particular, your graph does 
not cross the z-axis, then you will need to plot the graph of g again 
using a larger interval.) 

(416) find_root (g(x), x, 0, 1); 

(406) 0.21356701730612 

So a solution of sin?(« + 2) = 3x is z = 0.214 (to 3 s.f.). 


Solution to Computer activity 34 
(a) (4i1) trigreduce( sin(x)*2 ); 
(Hot) scene (a 


(b) (4i2) trigexpand( % ); 
sin(x)?-cos(x)?+1 
(ho02) —— 


(c) (4%i3) trigsimp( % ); 
(%03) sin(x)? 


Solution to Computer activity 35 


(Ail) trigrat( sin(4*x)/sin(x) ); 
(fot) 2cos(3x)+2cos (x) 


sin(4a) 


ae) = 2cos(3x) + 2cos(z). 


0, 
Solution to Computer activity 36 
(a) (Ail) £:cos(5*x) ; 
(hot) cos (5x) 
(/i2) g:trigexpand(f); 
(402) Scos(x)sin(x)*-10cos(x)%sin(x)?+cos(x)® 
So, cos(5a) = 5cos(x) sin*(x) — 10 cos? (a) sin?(«) + cos®(xr). 
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Solutions to Computer activities 


(b) Choosing the range —27 < x < 27 to plot the graphs, we obtain 
(4i3) Bihari 2 gen 2*%pil); 


ieinaiiaiia an ze 


i M 


|| 


| 


(ho3) 


The two graphs appear identical, which suggests the two expressions 
are equal. 
Solution to Computer activity 37 


(a) (4i1) load(implicit_plot) ; 
(hol) C: /PROGRA~1/MAXIMA~ 1. 0/share/maxima/5.30.0/share/contrib/implicit_plot.lisp 


(The output from this command is the location of the package on your 
computer system, and may differ from that shown here.) 


(b) (4i2) wximplicit_plot (x*2+y*2=1, [x,-2,2], [y,-2,2]); 
2 


— 7 T ye 


1.5 


° 
a 
T 


6 
a 
T 


a 
T 


(h02) 
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Solution to Computer activity 38 
(It is assumed the implicit-plot package has already been loaded.) 


(fil) wximplicit_plot (x*2+y*2=1, [x,-2,2], [y,-2,2], [gnuplot_preamble,"set size ratio -1"]); 
2 aa at fed ee 


> 


(htt) 


205i 10s OmOSma neste 
x 


(hol) 


Solution to Computer activity 39 
(It is assumed the implicit-plot package has already been loaded.) 
(a) 


(/i1) set_plot_option([gnuplot_preamble, "set size ratio -1"])$ 


(b) (4i2) wximplicit_plot((x-1)*2+y*2=3, [x,-1,3], [y,-2,2]); 


T T ae Se ia T T 
15 — ‘ 
1h, , 
f \ 
O5- / Nes 
\ 
| \ 
ee a } 
05} \ / 
4 i 
A5 bo = 4 
2 a Sa a 
105 005 1 15 2 25 3 
x 
(ho2) 
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Solution to Computer activity 40 


(It is assumed the implicit_plot package has already been loaded, and 
the plot settings changed to use equal scales for the x- and y-axes.) 
(a) (4i1) line:y=3*x-5; 

(hot) y=3x-5 

(%i2) circle: (x-2)%2+(y-3)*2=2; 

(%o2) (y-3)?+(x-2)?=2 


(b) (%i3) wximplicit-plot([line, circle], [x,0,5], [y,0,5]); 
5 


T =e 1 a ae 
y=3/x-5 —— 
(y-3)924(-2)'2 = 2 —— 
a er ‘ ; 
aie ul / 4 
\ a 
/ 
(43) \ 
Ait / | 
0 == 1 —-= 1 
0 1 2 3 4 5 
x 
(403) 
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Solution to Computer activity 41 


Solutions to Computer activities 


(It is assumed the implicit_plot package has previously been loaded and 


the plotting options set for equal-scale axes.) 


(a) (hil) wximplicit-plot(x*2+4+#y*2=1, [x,-2,2], [y,-2,2]); 
es 


T T T | aces T T 
155 
mall 
05+ — 
0 
(Ut) 
0.5 - —— 
| 
-1.5 = 
2 Ll Ll n a 
2-15 -1 05 0 05 1 15 
x 
(hot) 


(The curve is an ellipse.) 


(b) (4i2) wximplicit_plot (x*2-4+*y*2=1, [x,-4,4], [y,-4,4]); 
4 


eer aan at a a ca 
3b 
ar 


i- 


(402) 


(The curve is a hyperbola.) 


2 
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Solutions to Computer activities 


(c) (4i3) wximplicit_plot (y*2=x*3-2ex+1, [x,-4,4], [y,-4,4]); 
4 


een ee Se Se LT 
st a 


(43) 


(h03) 


(d) (%i4) wximplicit-plot (y*2=x*3-2*x+5, [x,-4,4], [y,-4,4]); 
4 


T a ioe | aes T Fs ja | 
3b ise 
ule | 
> OF 4 

(44) 
atte \ a 
a | 
3b " 4 
pe et 4 ; a 2 1 S| 
Vile eerie eae Yo | 

x 
(ho4) 
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(e) (445) wximplicit_plot(y*6=x*2-x*6, [x,-2,2], [y,-2,2]); 
2 


(f) 


i ee 
1.5 - 
iE 
pes a ee: —* 
0.5 + | / 
oF 
(45) | | | 
0.5 + | a \. } 
4 aw 
ae coon 
ave 
1.5 
2 Se 
2-15 -1 05 0 05 1 15 2 
x 
(405) 


(4i6) wximplicit_plot 
ie 


1.5 


(ht6) 


(406) 


‘15 -1-05 0 05 1 15 2 
x 


Solutions to Computer activities 


((x*2+y*2-1)*3=x"*2ey*3, [x,-2,2], [y,-2,2]); 
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Solution to Computer activity 42 


(a) 


(b) 


(e) 


(%il) diff (x*4+5*x*2+15, x); 
(hol) 4x3+10x 
So the derivative of 2* + 5x? + 15 is 4x7 + 102. 


(hi2) psti2/ain (ee) ; 
(02) aD 
(443) diff (p,t); 
(403) uy 2 Se 
re) 2t —_4t cos (4t) 


iatidw ok i = Ss 
So the derivative o sn) is sin(4#) sin? (4t) 


(414) gu) :=(u%2+3) *1log(u%2) ; 
(%04) g(a) :=(u?+3)1og(u?) 


(4iS) diff (g(a) ,u); 


2. 
(/05) 4ulog (a+? oe 


F 2 (uw? +3 
So the derivative of (u? +3) In(u?) is 4uIn (u) + Ss 


(4i6) diff (p,t,2); 
16 t? 2 16tcos(4t) 32t?cos(4t)? 
Caos? sin(4t) ss sin(4t) | sin(4t)? sin(4t)? 
ese ae 
So the second derivative of sind) is 
16¢? a 2 _ 16tcos (4t) 32t?cos? (4t) 
sin (4t) ' sin(4t) sin? (4t) sin? (4t) ~ 
(%i7) diff (g(u) ,u,3); 
2. 
on = 
u 


4 (u? + 3) 


So the third derivative of (u? + 3) In(u?) is a 
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Solution to Computer activity 43 
(a) (hil) diff (2*x*6+5/x, x); 
(hot) 12x6-5 
x 


ee 5 
The derivative of 2a° + — is 122° — —. 
x x 


(b) (4i2) diff (log(exp(x)+x), x); 


5 he*+1 
(402) vera 
x 
The derivative of In(e” + 2) is en y 
Cara 


(c) (443) diff ((tan(t)-t*2)/exp(t), t); 
(403) he (sec(t)?-2t)-%e* (tan(t)-t?) 
This can be simplified as follows. 
(444) fullratsimp(%) ; 
(404) ~he* (tan(t)-sec(t)?-t?+2t) 
(Remember, % in input line 4 refers to the result of the last evaluated 
command.) 


an(t) — t? 


Thetdenvartrelon -. is —e~* (tan (t) — sec (t)? — t? + 2). 


Solution to Computer activity 44 
(a) (hil) £(x) :=sin(x)/ (exp (x) +1); 


‘ ___ Sin(x) 
(hot) BS aGoEL 


(b) (4i2) df (x) :=diff(£(x),x); 
(402) df (x) :=diff (f(x) ,x) 


(ec) C4aS) deca) 
diff: second argument must be a variable; found 3 
#0: df (x=3) 
-- an error. To debug this try: debugmode(true) ; 
(d) (4i4) df(x):='' (diff (f(x) ,x)); 
Yorsi 
(od) af (x) gecesi) _he sinks) 
het+1 — (e*+1) 
_ cos(x) _ e* sin(x) 
e+1 (e= +1)? 


So f(x) 
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(e) 


(f) 


(%i5) df (3); 3 
7 cos(3) _e*sin(3) 
(h05) Gert1 (he2+1)° 
(416) float (%) ; 
(%06) -0.05332658917586 
So the value of f’ at « = 3 is —0.053 (to 3 s.f.). 


First, define the function. 


(hiT) £(x) :=sin(x) *exp(cos(x)) ; 
(ho7) £ (x) :=sin(x) exp(cos(x)); 
Then calculate the derivative, defining df(x) to be the result. 


(4i8) df (x) :='' (diff (f(x) ,x)); 
(408) df (x) :=%e®™ cos (x)-Keo*™ sin (x)? 


Solutions to Computer activities 


Finally evaluate the derivative at « = 1, giving the result as a decimal 


approximation. 

(4i9) float (df (1)); 

(09) -0.28798342608583 

So, the derivative of sin(x)e°*) at a = 1 is —0.288 (to 3 s.f.). 


Solution to Computer activity 45 


(a) 


(b) 


(Ail) integrate (2*x*2+exp(3*x), x); 


Ye 2x7 
hor + 
(hot) aaa : . 
pene! 
Sof 22?+e% de = +p te 


(4i2) integrate(1/x, x, 2, 3); 
(ho2) log (3)-log (2) 

(4A3) float (%) ; 

(403) 0.40546510810816 


3 
So if fae = In(3) — In(2) = 0.405 (to 3 s.f.). 
2 


(hi 4) pCa) 1=x/ (e241) ; 
y ——— 
(hod) pbx) = 


(%i5) q(x) :='' Cintegrate (p(x) ,x)); 


2. 
(Ho5) q(x) p= Doge st) 


2 
so f par= MFO 46, 
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Solution to Computer activity 46 


(a) (4i1) integrate ((x-3)*(x-1), x); 
x5-6 x74+9 x 
Ce ———— 


cays 
80 f(x —3)(«-1) ax = eae ae. 
(b) (412) integrate(1/sqrt(9-t%2), t); 


(ho2) asin(*) 


So ee =sin™ (5) +e. 


Vv9-# 3 
(c) (4i3) integrate(exp(t)*sin(t), t, 0, %pi); 
he™ 1 
(413) io 


Dae 
(hid) float (%); 
(404) 12.07034631638963 


7 


es e 1 
So f e'sintdt = — + = = 12.1 (to3s.f.). 
0 PAE 


Solution to Computer activity 47 
(a) (4il) integrate(exp(sin(x)), x); 
(hot) [tePax 
Maxima was unable to calculate this integral. 
(b) (4i2) integrate(exp(-x%2), x); 
ea) SS 
Maxima has calculated this integral, in terms of a function called the 
Error function and denoted erf. 


(c) (4i3) integrate(exp(x)/x, x); 
(403) -gamma_incomplete(0,-x) 
Maxima has calculated this integral, in terms of a function called the 
incomplete Gamma function. 
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Solution to Computer activity 48 


(a) (i) (hat) kill @); 
(fol) done 


(ii) (412) integrate(x*n, x); 
Is nt1 zero or nonzero? n; 
ett 
(402) — 
nt+1 


(b) (i) (443) assume (notequal (n,-1)) ; 
(403) [notequal(n,-1)] 


(ii) (4i4) integrate(x*n, x); 


n+l 
(404) ~ 


nt+1 


(iii) (415) factsQ; 
(405) [notequal(n,-1)] 
So Maxima knows n # —1. 


(iv) (416) forget (motequal(n,-1)); 
(406) [notequal(n,-1)] 


(v) (4i7) facts(; 
(it) O 
No facts are now known. 
Solution to Computer activity 49 


First, tell Maxima what we know about a. 


(Ail) assume(a>0) ; 
(401) [a>0] 
Then do the integration. 


(412) integrate (1/(x*2ta) ,x); 


where a > 0. 


Solutions to Computer activities 
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Solution to Computer activity 50 
(a) (4i1) integrate (log(cos(x*2)), x, 0, 1); 
(hot) i log (cos (x?) )dx 
Maxima cannot evaluate this integral exactly. It returns the integral 
unevaluated. - 
(b) (412) quad_qags(log(cos(x*2)), x, 0, 1); 
(402) [-0.11151210243615 ,3.3310015826528369 10°" 21,0] 


(Remember 3.3310015826528369 107'? means 
3.3310015826528369 x 10-!*.) 


1 
So | In(cos(a”)) da = —0.112 (to 3 s.f.). 
0 


(c) (4i3) quad_qags(exp(-t*2), t, 0, 1); 
(403) [0.74682413281243 ,8.2914134759407266 10°*®, 21,0] 


1 
so | ce” dt = 0.747 (to 3 sf). 
0 


Solution to Computer activity 51 
(a) (4i1) A: matrix( [1,3,5], [2,4,6] ); 


to ac 
(hot) iF 4 A 


(b) (412) B: matrix( 
Ei 2 735495. 
[4,1,2,3], 
(2,4,1,3], 
) 


1 
(h02) |4 
2 


Solution to Computer activity 52 
(a) (4i3) matrix_size(A) ; 
(403) [2,3] 
So the matrix A has size 2 x 3. 
(b) (4i4) BIS,2]; 
(ho4) 4 
The element in the third row and second column of B is 4. 
(c) (415) BI3]; 
(405) [2,4,1,3] 
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(4) 


(%i6) B[3,2]:0; 
(%06) 0 
CHLv Bs 

r 


1 23 4 
(4o7) |4 1 2 3 
[20 1 


Solution to Computer activity 53 


(a) 


(c) 


(fit) P:matrix([2,0] , [3,-1], [-4,6]); 


2 0 
(401) | 3-1 
|-4 6 


(4i2) Q:matrix((0,3],[5,-4],[1,3]); 
0 3 
5 -4 
ae 


(ho2) 


(4i3) Rimatrix([1,3], (2,41); 
1 3 
(ho3) |n 4 


(#44) S:matrix([1], [2]); 


(Hoa) |3 


(4i5) T:matrix([3,4]); 
(405) [3 4] 


(4i6) P+Q; 


2ug3 
(406) | 8-5 
[-3 9 


2 3 
SoP+Q= 8 -5]. 
—3 9 


(hiT) P-R; 
fullmap: 
-- an error. 


arguments must have same formal structure. 
To debug this try: debugmode(true) ; 


P —R cannot be calculated, since the matrices are of different size. 


An error is given. 


Solutions to Computer activities 
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(d) (4i8) 2*P; 


4 0 
(os) 6 =2 
ee de 


4 0 
so2P= 6 —2'}. 
—8 12 


(e) (4i9) P-R; 


(f) (4i10) P.Q; 
MULTIPLYMATRICES: attempt to multiply nonconformable 
matrices. 
-- an error. To debug this try: debugmode(true) ; 


PQ cannot be calculated, since the matrices are not of the 
appropriate sizes. An error is given. 


(g) (4i11) P.S; 


(%o11) 


(h 


(4112) P.T; 

6 
(4012) | 5 
12 


6 
Maxima interprets the row vector T as a column vector and obtains ( ) . 
12 


(i) (4413) R42; 


mi q 435 
(%013) B A 


7 15 
rie 
So R* = i = i 
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(j) (#i14) R42; 


F rere: 
(4014) Fe a 


Maxima returns the matrix (; a) whose elements are the squares 


of corresponding elements of R. 
(k) (4i15) P*Q; 


0 © 
(4015) |15 4 


-4 18 
0 0 
Maxima returns the matrix | 15 4] whose elements are the 
—4 18 


products of corresponding elements of P and Q. 


Solution to Computer activity 54 
(a) (441) A:matrix([1,2,0], [-4,4,2],[-1,2,1]); 


pL ey ty 
(hol) |-4 4 2 


(b) (4i3) determinant (A) ; 
(403) 4 
So det A = 4. 


(c) (4i4) C:invert (A); 
o -i 1 


2 
(oa) | 2-2 
Sele es 


So AS *— a 


Boar ie 
oo hi 


(d) (4i5) A.C; 


100 
(405) JO 1 0 
oo 1 


So AC = CA =I, and C is the inverse of A. 
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(e) (4i7) invert (B) ; 
expt: undefined: 0 to a negative exponent. —- an error. 
To debug this try: debugmode(true) ; 


B does not have an inverse. 


Solution to Computer activity 55 


(a) First, enter the matrix into Maxima, and assign it to a variable. 
(hit) A:matrix([2,-1,1],[2,0,1],[4,2,1]); 


ie al ak 
(Ket) |2 0° 4 
A 2 tL 


Then, check if the matrix is invertible by finding its determinant. 


(4i2) determinant (A); 
(ho2) -2 


We have 


Pie aE 
det}2 0 1) =-2, 
oer aa 


so this matrix is invertible. 
Now find the inverse. 
(443) invert (A) ; 


{oon 

a 2 

(403) |-1 1 0 
<2; 21 


Dy Siti SS 1-2 2 
Soil|ioe i070) ==||/=shoee sil: 
ae Soe eee 


(b) Again, first enter the matrix and calculate the determinant. 
(4i4) B: matrix( [-1,0,3,-2], [4, 1/2, -1 ,0], [2,0,1,-1],[-1,2,1,-1]); 


Oh a) 2 
; ae ert a 
ho) 15 6 4 1 
Sal tol a 


(4i5) SOLED, 
(405) ae 
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We have 
-10 3 -2 
44-1 0 1 
LA Sy tis aieeest|| eh 
—1 2 i -1 


so this matrix is invertible. 
Now find the inverse. 
(416) invert (B) ; 


ie ot 
6 -6 10 2 
(h08) | 19 -00 33 5 
-27 -28 46 7 
Stay 8 Sy ee ea, ee 
3 43-1 0) _|{-6 -6 10 2 
ase ail) ow [lay 22h) G5 
=12) el =oF 2846 7 


Solution to Computer activity 56 


(a) 


(b) 


(4i1) aln) :=0.8%n; 
(hol) ay:=0.8" 


(%i2) a[100] ; 
(402) 2.0370359763345081 10°*° 
(Remember, this means 2.0370359763345081 x 10~1°.) 


So a1o9 = 2.037 x 10-19 (to 4 s.f.). 


Solution to Computer activity 57 


(a) 


(b) 


(hii) bL1] 21; 

(hot) 1 

(4i2) b[n] :=2*b[n-1]+1; 
(402) by:=2 by-1+1; 


(%i3) b[5]; 

(%o3) 31 

(%i4) b[100] ; 

(%04) 1267650600228229401496703205375 

So, bs = 31 and bio9 = 1267650600228229401496703205375 
(%i5) b[1] :10; 

(%05) 10 

(%i6) b[5]; 

(406) 31 
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(d) (4i7) kil1(b) ; 
(h07) done 
(448) b[1]:10; 
(%08) 10 
(%i9) bn] :=2*b [n-1] +1; 
(%09) by:=2bp-1 +1; 
(%i10) bf5]; 
(4010) 175 
So bs is now 175. 


Solution to Computer activity 58 


(a) (4i1) c{n) :=n*2; 
(401) cy:=n? 


(b) (442) pts:( [1,cf1]], (2,cl2]], (3,c(3]] 1; 
(02) [f1,1], [2,4], [3,91] 


(c) 
(/i3) wxplot2d([discrete, pts], [style, points], [xlabel, "n"], [ylabel, "c{n]"]); 
9 T T T 
8 
if 
a 6 
& 5 
Cht3) 9 4 é 
3 
2 
1 1 1 1 
1 15) 2 25 3 
n 
(403) 


Solution to Computer activity 59 


(a) (hit) a1] :0; 
(hol) 0 
(442) d[n] :=1/2*(2-d[n-1]) ; 


(%02) dy: = (2-dy-1) 


(b) (443) pts:makelist([n,d[n]], n, 1, 15); 


i 1 3 5 11 21 43 85 171 
(03) ((1,0], [2,1], (3,51, (4,71, [8,5], (6,761. (7,55) (8,71, (9, 7581) (10, 555] 
341 683 1365 2731 5461 
115) (12.5554) (13,5048! (14.7096 (15,7591 
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(c) 
(4i4) wxplot2d([discrete, pts], [style, points], [xlabel, "n"], [ylabel, "d[n]"]); 
1 — T T T T T 
0.8 - 
_ 06 oS SereTe (er eie/e 
= e 
(ta) © pe 
0.2 
0 A 1 fe ney 1 
0 2 4) 6 By 10 42 14) “16 
n 
(04) 


Solution to Computer activity 60 
(a) First, define the sequence. 
(fat) elisa 
(hot) 1 
(%i2) r[n]:=1.5*r[n-1] - 0.1*(r[n-1])*2; 
(402) ryet@1.5r,-1-0.1 x2, 
Next, create a list of points to plot, suppressing the output. 


(4i3) rpts:makelist({n,r[n]] ,n,1,50)$ 
Finally, plot the graph. 


(4i4) wxplot2d([discrete, rpts], [style, points], [xlabel, "n"], [ylabel, "r{n]"]); 


= 


(ht4) = 


1 
0 5 10 15 20 25 30 35 40 45 50 
n 


(h04) 


The terms of the sequence seem to be getting closer and closer to the 
value 5. 
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(b) First, define the sequence. 
(415) s[n] :=n*3-10; 
(%05) s,:=n3-10 
Next, create a list of 50 points to plot, suppressing the output. 


(%i6) spts:makelist([n,s[n]] ,n,1,50)$ 
Finally, plot the graph. 


(447) wxplot2d([discrete, spts], [style, points], [xlabel, "n"], [ylabel, "s[n]"]); 
140000 So 
120000 
100000 
80000 
60000 
40000 
20000 
i) 
-20000 


s[n] 


Cht7) 


0 5 10 15 20 25 30 35 40 45 50 


n 


(407) 


The terms of the sequence seem to be continually increasing. 
(c) First, define the sequence. 
(%i8) tn) :=1-1.14n; 
(%08) t,:=1-1.1" 
Next, create a list of 50 points to plot, suppressing the output. 
(/i9) tpts:makelist({n,t[n]],n,1,50)$ 
Finally, plot the graph. 
(4i10) wxplot2d([discrete, tpts], [style, points], [xlabel, "n"], [ylabel, "t[n]}"]); 
0 — 


t(n] 
3 


(%t10) 


0 5 10 15 20 25 30 35 40 45 50 
n 


(4010) 


The terms of the sequence seem to be continually decreasing. 
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Solution to Computer activity 61 
(a) (%i1) sum(n*2, n, 1, 100); 


(fot) 338350 
100 
So > n? = 338350. 
n=1 
20 


n 
(b) The sum can be written by (5) and so can be calculated as follows. 
n=1 


(c) 


(d 


(412) sum((1/3)%n, n, 1, 20); 
- 1743392200 
(402) 


3486784401 
P S 1)" _ 1743392200 
) =| ==> 
iat 3 3486784401 
Notice that the sum is a geometric series with first term, a = 3, 
common ratio, r = $ and with 20 terms. So the sum can be calculated 


using 
20 
a(i—r") _ 3(1 - (3) ) 
1-r 1-3 


ll 
tol cantoleore 
F cite e. 4 
_ 
| 
oe 
| + 
8 
~~ 


379-1 
(Ss) 
370-1 
a x 320 


3486784400 


2 x 3486784401 
1743392200 


~ 3486784401 
~ 0.4999999998566. 


(413) sum(1/(4%n), n, 1, inf); 
y - 1 
(h03) Sm 


n=1 
(414) simpsum:true; 
(f04) true 
(hi5) GR. > n, 1, inf); 
(%05) 3 
it 1 
So ms You could check this result using the expression for 
=1 
the sum of an infinite geometric series. 


00 
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(e) 


Solutions to Computer activities 


(%i6) sum(1/(n*2), n, 1, inf); 


Solution to Computer activity 62 


(a) 


(b) 


(c) 


(%i1) taylor(cos(x), x, %pi/6, 3); 


m2 1r\3 
vi_v_V8(xs) , (3) 
2 a 4 12 
So the cubic Taylor polynomial about 7/6 for for f(x) = cos(z) is 


ww_2-% v8(e-2)"_(e-2) 
o Pl 4 12 


(401) /T/ aR 


(4i2) p(x)='' (taylor(tan(x), x, 0, 5)); 
7 pel Oe 
(h02)/T/ p(x) Peis + Fei te 
So the quintic Taylor polynomial about 0 for f(x) = tan(r) is 
en oe 


(hi3) ply); a 
(03) /R/ = 


(%i4) float (p(0.1)); 

rat: replaced 0.10033466666667 by 75251/750000 = 0.10033466666667 
(404) 0.10033466666667 

So, p(0.1)= 0.10033 (to 5 s.f.). 

(/i5) float (tan(0.1)-p(0.1)); 

rat: replaced 0.10033466666667 by 75251/750000 = 0.10033466666667 
rat: replaced 0.10033467208545 by 9334601/93034649 = 0.10033467208545 
(405) 5.418783991614422 10° 

So the difference between the function and its quintic Taylor 

polynomial at x = 0.1 is 5.42 x 10~® (to 3 s.f.). 
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Solution to Computer activity 63 
First, define the function. 


(4il) £(x) :=exp(sin(x%2)) ; 

(Hol) f(x) :=exp(sin(x?)) 

Next, calculate the required Taylor polynomial. The command is ended 
with $ to suppress displaying the result. 


(H4i2) px) :="' (taylor (£ (x) ,x,1,5))$ 
Finally, plot the two curves over the required range. 


(%i3) wxplot2d([f(x),p(x)], [x,0,2], [legend, "£(x)", "p(x)"]); 


4 : : 7 
2 oe 
0 a = 
a 
(ht3) 4 
6 
8 
-10 4 4 —___1—— 
0 05 1 15 2 


(403) 


Solution to Computer activity 64 


(a) (4il) a:3+4*%i; 
(hot) 4%i+3 


(b) (442) b:2-2«hi; 
(%o2) 2-2 %i 


(c) (hi3) atb; 
(403) 2hi+5 
(hid) a-b; 
(ho04) 6 %i+t 
So, atb= 5 + 2i and a-b= 1 + 6i. 
(hi) ab; 
(ho5) (2-2 %i) (4 %i+3) 
(%i6) a/b; 
hit3 


i 4 
(h06) Soni 
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Solutions to Computer activities 


(d) (4i7) rectform(a*b) ; 
(ho7) 2%i+14 


TH + 


So, a*b= 14 + 2i and a/b= —} + fi. 
(e) (4i9) realpart (a); 
(ho9) 3 
(%i10) imagpart (a) ; 
(ho10) 4 
So the real and imaginary parts of a are 3 and 4 as expected. 


(f) (4i11) conjugate (a); 
(holt) 3-4%4 


So the conjugate of a is 3 — 4i as expected. 


(g) (%i12) abs(b); 
(4012) 29/2 
(%i13) carg(b); 
(ho13) -2 


So, |b] = 23/2 = (v2)° and the principal argument of b is =F 


Solution to Computer activity 65 
(a) (hil) c:3*(cos(5*/pi/13)+%i*sin(5*%pi/13)) ; 


a eS ee 
(hot) a(n sin (=) + cos (2) ) 


(4i2) d:4*(cos(7#*/pi/13)+%i*sin(7*%pi/13)) ; 
(02) 4(4i sin (2) +cos (™ 
hi hi sin| 7, cos | 73 
(b) The modulus of c*d can be found as follows. 
(/i3) abs (cd) ; 
2 2 3 2 
ca) s2y/sin (52) seas (82) yfsin (72) sean (72) 


Simplify this using the trigreduce command. Remember, % refers to 
the result of the previous calculation. 


(4i4) trigreduce(%) ; 
(404) 12 
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The principal argument of c*d can be found as follows. Again the 
initial result requires simplification. 


(4iS) carg(cxd); 


3 as = yg 
sini | —= sin | — 
13 13 

Tr Sa 

cos | — cos || — 
13 13 


(4i6) trigreduce(%) ; 


(05) atan 


" 127 
(%o06) er 
So, eed] = 12 and Ae (ora) = am 


Similar methods can be used to find the modulus and principal 
argument of ¢/d. 


(hiT) abs(c/d); 


(4i8) trigreduce(%) ; 
3 
(408) z 


(%i9) carg(c/d); 


es is 
sin {| — 
13 


(4i10) inter emae ; 
7 
(4010) =a 


(409) ~atan 


Note that although carg gives a principal argument, that is an 


Solutions to Computer activities 


argument in the interval (—7, 7], the trigreduce command simplifies 
the expression to an equivalent argument which is outside this interval. 


naa . 240 2a 
The principal argument is 3 2 = “3° 
; 2: 
So, |c/d| = 2 and Arg(c/d) = — 


168 


Solution to Computer activity 66 
(a) (il) solns:solve(z*3+z+1=0) ; 


vB%i 1 
ee ee i en ae 
23372 2 2 ) (= 
3 


237 2 

— 1 

(sy v3%i_1 is -- 
2gt/2 a peo 

(eh 


——— 7! 
V31 1 1/3 
8 23972 2 
(b) (%i2) float (rectform(solns)) ; 

(%02) [z=0.34116390191401-1.161541399997252 %i, 


z=1.161541399997252 %i+0 .34116390191401, 
2=-0.68232780382802] 


So, to two decimal places, the solutions are 


z=0.34-1.16%, z=0.34+1.16i, and z=-—0.68. 


Vat ay 1 
*\232 2) 


(c) (4i8) solve (z*4-20*z"3+171*z"2-626*z+962=0) ; 
(403) [z=3-2 %i, z=2 %i+3,z=7-5 %i,z=5 4i+7] 
So the solutions are 3 — 2i, 3+ 2i, 7 — 5i and 7 + 5i. 


Solution to Computer activity 67 
(a) (4i1) s:solve(z*4+4*z+5=0) ; 


Solutions to Computer activities 


(hot) tees | a, aon pe SALES sa 


(b) (%i2) v:makelist(rhs(s[k]), k, 1, length(s)); 
(02) t-ya- es oui, 4-4 os 


i 2 


(c) 


+i) 


(443) pts:makelist ([realpart(v[k]), imagpart(v[k])], k, 1, length(s)); 


V¥28/244 /28/2-4 V¥28/244 _\/28/2-4 


(ho3) [l-—Saa— Samm Pale maesS/Oi= 2) aaa ema 2 
28/244 V'28/2-4 25/244 \/28/2-4 
7 1g3/2 4+ ps8 7? 9872? 9372 +11] 
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Solutions to Computer activities 


e (414) wxplot2d([discrete,pts], [style,points], [xlabel,""], {ylabel,""]); 
1.5 = is T = ow. 
1b 
0.5 F e 
(ht4) OF 
0.5 - e 
Ae 
15 n 1 — = 12 


(h04) 


Solution to Computer activity 68 
(a) (4i1) s:solve(z*7=1) ; 
2 i 1 4hin 6 i + 6 hi 
(Hot) [z=%e 7 , zehle 7 , zxhe 7 , zeke 7 


4 hin 2% x 
z=he 7 +, ze 7 , z=1) 


The solutions are given in terms of exponential functions. You will 
learn about the exponential form of complex numbers in Section 4 of 
Unit 12. 


The solutions can be written in Cartesian form as follows. 


(%i2) float (rectform(s)) ; 

(402) [z=0.78183148246803 %i+0.62348980185873, 
z=0.97492791218182%i-0 . 22252093395631, 
z=0 .43388373911756%i-0 . 90096886790242, 
z=-0.43388373911756%i-0 . 90096886790242, 
z=-0.97492791218182%i-0 .22252093395631, 
z=0.62348980185873-0. 78183148246803%i, z=1.0] 


So, to two decimal places, the seven seventh roots of unity are 
z=0.62+0.78i, z= —0.22+0.97i, z= —0.90 + 0.43%, 
z= -0.90—0.43i, z=-—0.22—0.97i, z=0.62 —0.78i 
and z = 1.00. 
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Solutions to Computer activities 


To plot the solutions in the complex plane, first create a list 
containing just the solutions of the equation (without the z=). The 
command is ended with $ to suppress the output. 


(4i3) vimakelist (rhs(s[k]), k, 1, length(s))$ 


Then create the list of points to be plotted. 
(Ai4) pts:makelist([realpart(v[k]), imagpart(v[k])], k, 1, length(s))$ 


Finally, plot the points. 
(hi5) wxplot2d([discrete,pts], [style,points], [xlabel,""], [ylabel, ""]); 


1 Se 
0.8 ° 


(%t5) 0 


-0.8 e 
4 (Seees oe ee eal pe re eh 


-1 08-06 -04-02 0 02040608 1 


(%o5) 


Notice that the complex solutions occur in conjugate pairs, and hence 
lie symmetrically about the horizontal axis. The solutions also lie on 
the unit circle (which does not seem circular here due to the different 
axis scales used). 


171 


Maxima reference guide 


Mathematical operations and functions 


Maxima reference guide 


Operation or function Syntax Example 
Addition + 443; 
Subtraction m= 4-3; 
Multiplication * 4*3; 
Division i 4/3; 
Brackets Cand ) 2* (3+4) ; 
Powers, for example 2° “ or ** 243; 
Q**3; 

Square root, for example V5 sqrt (() sqrt (5); 
Exponential, for example ee re Kero: 

or, exp(§) | exp(3); 
Natural logarithm, In log () log(8) ; 
Magnitude (absolute value) of a real number, | abs ({) abs (-3) ; 

for example | — 3|, or modulus of abs (1+2*%i) ; 
a complex number 

sin sin(™® sin(1); 
cos cos () cos (3*%pi/2) ; 
tan tan(§) tan(%pi/4) ; 
cosec csc() csc(2); 
sec sec((§) sec(%pi/3) ; 
cot cot () cot (%pi/4) ; 
sin~! asin(§™) asin(sqrt(3)/2) ; 
cos! acos (jf) acos(1/sqrt(2)); 
tan7! atan(™§) atan(1/2); 
Equality (within an equation) = 5=2*x+3; 
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Maxima reference guide 


Other mathematical symbols 


Operation or function Syntax | Example 
Matrix addition + A+B; 
Matrix subtraction = A-B; 
Scalar multiplication * 2*A; 
Matrix multiplication : A.B; 
Matrix powers on AS: 
Determinant (of a square matrix) determinant (matrix) determinant (A) ; 
Matrix inverse (of a square matrix) | invert (matrix) invert (A) ; 

or tatrix’’(-1) A” (-1); 
Identity matrix ident ((Bizé)) ident (2) ; 
Real part of a complex number realpart (j§) realpart (1+2*%i) ; 
Imaginary part of a complex number | imagpart ((™)) imagpart (1+2*%i) ; 
Conjugate of a complex number conjugate ()) conjugate (1+2*%i) ; 
Argument of a complex number carg(™) carg(1+2*%i) ; 


Symbol | Syntax | Example 


e he he2; 

T dpi 2Q*pi; 

i hi 2+3*hi 5 

co inf sum(1/n*2,n,1,inf); 
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Operation Command Example 
Get help on a command ? (command ? float; 
or, describe ((€emimiand! ) describe (float) ; 
Find help information whose 7? Hf ?? flo; 
title contains the given text or, describe(§, inexact) | describe(flo, inexact); 
Force the evaluation of a "@ g(x) :='' (diff (x%4,x)); 
command 
Load a package load ( (package! mame) load(implicit_plot); 
Assign a value to a variable : a:23; 
Display the value of a variable variable a; 
Define a function := f(x) :=2*x+3; 
Evaluate a function at a value ‘function ( ) £0 
Remove an assigned variable kill (Waray) kill(a); 
or function kill (fanetion) kill(£); 
Remove all assigned kill (all) kill(all); 
variables and functions 
Reset all system variables reset () reset (); 
Convert to a decimal number float (§) float (sqrt(2)); 
Expand brackets expand (jf) expand( (x+1)2 ); 
Factorise factor (™) factor( 2*x+4*x‘2 ); 
Simplify fullratsimp(§) fullratsimp( (2*x+4#x%2)/x ); 
Simplify something involving radcan({) radcan( log(x*2) ); 
exponentials and logarithms 
Combine logarithms logcontract (() logcontract (log(a)+log(b)) ; 
Expand trigonometric functions trigexpand (jj) trigexpand(sin(A+B)) ; 
of sums, differences and 
multiples of angles 
Express powers of sines and trigreduce (jj) trigreduce(sin(x)%2); 
cosines in terms of sines and 
cosines of multiple angles 
Simplify trigonometric expressions | trigsimp (jj) trigsimp(sin(x)”*2+cos(x)%2) ; 
Simplify algebraic fractions trigrat (™) trigrat ((sin(2x))/cos(x)); 
containing trigonometric 
functions 
Substitute subst (Wale) , Wardabie) , subst(4, x, x*2+1); 
‘expression ) which substitutes 4 for x in x2+1 
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Operation Command 
Left-hand side of an equation _| Lhs ((@qiaitiien!) 
Right-hand side of an equation | rhs (@quation)) 
Solve an equation, exactly solve (equation ) 
Solve an equation for a variable | solve (@quatien , (wardabile)) 
Solve simultaneous equations solve (Mist |6# equations’ , 
‘list of unknowns ) 
Solve an equation numerically | £ind_root ((@xpressiion', 
variable , 
interval start value , 
interval end value) 
Make an assumption about a assume (() 
variable 
State two things are equal equal (jj, 
State two things are not equal | notequal (/§,) 
Forget a property forget (property) 
List all known facts facts() 
List all known facts about facts (Wariable)) 
a particular variable 
Form a list CE. 8. --- 1 
Create a list makelist (generallterm , 
dummy variable , 
‘start value, end value) 
Extract an element of a list ‘list [index ] 
Find the length of a list length (a8#)) 
Specify a matrix matrix (ow, Fow,...) 
Size of a matrix matrix_size (matrix) 
Row of a matrix matrix [row] 
Element of a matrix matrix [row , column] 
or maerise (fou (leeinit] 


Example 


lhs (4*x+1=2*x-2) ; 

rhs (4*x+1=2*x-2) ; 

solve(2*x*2-1=0) ; 

solve(2*a*b-3*b=0, a); 

solve([2*x+y=4, 3+#x-2*y=-1], 
—x,y]); 

find_root (2*x*2-1,x,0,1); 


assume (a>0) ; 


assume (equal (n,3)); 
assume (notequal(n,-1)); 
forget (a>0) ; 

facts(); 

facts(a); 


A: fa,b,c]; 
makelist (2+n,n,1,100); 


A[2]; 
length(A) ; 


A:matrix([1,2] ,(3,4]); 
matrix_size(A); 

A[2]; 

A[3,4]; 

A[3] [4]; 
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Operation 


Command 


Plot a graph 


Plot several graphs 


Plot a sequence 


Plot a curve represented 
by an equation in implicit 
form 


Plot several curves 
represented by equations 
in implicit form 


Set the default plotting options 


Differentiate 
Differentiate multiple times 


Integrate 
(find an antiderivative) 
Evaluate a definite integral 


Find an approximate value 
of a definite integral 


Example 


wxplot 2d (SSRGSSSH 
horizontal range ,...) 


(use plot2d similarly if not 
using wxMaxima) 


wxplot2d ( IiSt/6f expressions), 


Orizontalyrangs,...) 


wxplot2d([discrete, 
list of coordinate pairs], 


(style, points],...) 
wximplicit_plot ae. 
: . 


wximplicit_plot( 


List of equations , 
Wertiealrange,...) 


(use implicit_plot similarly 
if not using wxMaxima) 


set_plot_option (option) 

ditt), WARIS) 

ditt (\§, WAEKABIS , 
positive integer ) 

integrate (\\j, Variable!) 

integrate (fj, Variable), 


lower limit, upper limit ) 
quad_qags (jj, Wardablel, 
, ) 


wxplot2d(x*2, [x,0,1]); 
which plots the graph of x? 
for « between 0 and 1. 


wxplot2d(x*2, [x,0,1], [y,0,2]); 


which plots the graph of x? 
for x between 0 and 1 


and vertical axis values between 0 and 2. 
wxplot2d([x*2, 2*x], [x,0,1]); 


which plots graphs of «? and 2x 
for « between 0 and 1. 
wxplot2d([discrete, 
((1,5], [2,10]]], 
[style, points]); 
wximplicit_plot (x2+y*%2=1, 
[x,-1,1], ly,-1,1]); 


wximplicit_plot ( 
([x*2+2*y%2=4, y=x°3], 
(x,at dy, fys-tet)s 


set_plot_option( 
(gnuplot_preamble, 
"set size ratio -1"]); 

diff (sin(x*2) ,x); 

diff (log(x) ,x,3); 


integrate (x*2,x); 
integrate(sin(x),x,0,1); 


quad_gqags (exp(-x2) ,x,0,1); 
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Maxima commands (continued) 


Operation Command Example 
Define a sequence ‘sequence [n] :=expression in n a[n] :=n“2; 
using a closed form ‘ 
Define a sequence SRGRGHCE (HAFSECSERTMRBEE] -GMiENANWANe §=©— 1] :1; 
using a recurrence ‘sequence [n] := expression in sequence [n-1] b{n] :=2*b[n-1] ; 
system 
Calculate a term in a ‘sequence [term number ] b[50] ; 
sequence 
Sum a series sum( general term, index variable, sum(n*2,n,1,10); 


lower limit , upper limit ) 
Find a Taylor polynomial | taylor (#aHEtion, Wariable), Centre|, Gegrés|) | taylor(sin(x),x,0,3); 


Express a complex number | rectform((§) rectform(2/%i) ; 
in Cartesian form 


Se a 


Maxima system variables 


Variable 


Description Example 


The number of significant figures of 
a decimal number to display 
functions The list of all user defined functions 
values The list of all user assigned variables 
simpsum Whether to simplify summations 


fpprintprec fpprintprec:3; 


functions; 
values; 
simpsum:true; 
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The following options can be added as additional arguments to plotting 


commands. 
Option Argument Example 
Vertical range ly. &. [y, 0, 5] 
Curve colour [color, {color, red] 
Legend text [legend," el [legend, "Car A"] 


Turn legend off 

Horizontal axis label 

Vertical axis label 

Plot points 

Change the relative lengths/scales 
of graph axes. 

(A positive number specifies the 
ratio of axis lengths, a negative 
number specifies the ratio of axis 
scales.) 


[legend, false] 
{xlabel, "§§"] 
[ylabel, "§§"] 
[style, points] 
[gnuplot_preamble, 
"set size ratio Humber") 


(a label should be listed 
for each curve plotted) 
[legend, false] 
[xlabel, "time, t (h)"] 
[ylabel, "displacement, s (km)"] 
(style, points] 
{gnuplot_preamble, 
"set size ratio -1"] 
(which sets the axes to have 
equal scales) 


Maxima packages 


Maxima package | Functions added 


implicit_plot | Plotting curves given by equations in implicit form 


wxMaxima keyboard sequences (Microsoft Windows) 


Zoom in 

Zoom out 

Copy 

Interrupt a calculation 
Re-evaluate the worksheet 
Close wxMaxima 

Paste 


Key sequence Function 
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wxMaxima menu commands (Microsoft Windows) 


Restart Maxima Maxima 
Maxima help window | Help 
Re-evaluate worksheet | Cell 


Insert text cell Cell 
Insert title cell Cell 
Insert section cell Cell 


Insert subsection cell | Cell 
Enter a matrix 


Operation Menu name | Menu option 
Close wxMaxima File Exit 

Save your work File Save or Save As 
Print the worksheet File Print 

Open a worksheet File Open 

Configure wxMaxima | Edit Configure 
Copy Edit Copy 

Paste Edit Paste 

Interrupt a calculation | Maxima Interrupt 


Algebra 


Restart Maxima 

Maxima Help 

Evaluate All Visible Cells 
Insert Text cell 

Insert Title cell 

Insert Section cell 
Insert Subsection cell 
Enter Matrix... 


Maxima command-line interface commands 


Turn off 2D formatting of maths 

Display output left justified 

Start writing a transcript 

Stop writing a transcript 

Playback all the input and 
output in a session 

Playback a range of input and 
output 

Save the state of a session 

Reload a session 


Operation Command 
Close Maxima quit(); 
Interrupt a calculation Ctrl-C 


display2d:false; 
leftjust:true; 
writefile("$ilepath") ; 
closefile(); 
playback(); 


playback( [Start/linelntimber, jend/inejmimber) ); 


save ("(filelpath", all); 
Load ("aIeypath") ; 
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Index 


abort calculation 14 

using command-line interface 100 
abs 15 

for complex numbers 89 

using the command-line interface 101 
absolute value 15 

using the command-line interface 101 
accessibility 94 
acos 47 
addition 15 

of matrices 74 

using the command line interface 101 
approximation 

numerical 50 
argument of a command 15 

using command-line interface 101 
argument of a complex number 88 
asin 47 
assign 21 

using command-line interface 106 
assume 68 


assumption 
about a variable 68 
equals 68 
forget 68 
list all 68 
not equals 68 
atan 47 
carg 89 


Cartesian form of a complex number 88 
CAS _ see computer algebra system 
cell 12 
marker 11, 12 
changing a sequence 80 
circles 
plotting 55 
closed form of a sequence 78 
closing the command-line interface 101 
closing wxMaxima 15 
command-line interface 8, 97 
input line number 98 
closing Maxima 101 
editing cursor 98 
evaluating an expression 98 
interrupt 100 
output line number 98 
plotting graphs 115 
reusing commands 104 
saving work 110 
scientific notation 104 


complex numbers 34, 87 
Cartesian form 88 
argument (principal) 88 
conjugate 88 
imaginary part 88 


modulus 88 
polar form 89 
real part 88 


computer algebra system 5 
configure wxMaxima 13, 16 
conjugate 89 
conjugate of a complex number 88 
copy 20 
cos 47 
cosec 47 
cot 47 
create a list using makelist 83 
csc 47 
cube roots 16 
using command-line interface 102 
cursor 
editing 19 
using command-line interface 98 
horizontal 20 


decimal approximation 17 

using command-line interface 103 
define see assign 
definite integral 64 

approximate value 69 


delete 
function 43 
variable 23 


variable, using command-line interface 
derivative 60 
assign to function 63 
second, third 61 
determinant 76 
determinant of a matrix 76 
diff 60 
differentiate 60 
division 15 
using the command line interface 101 


editing commands 19 

using command-line interface 104 
editing cursor 19 

using command-line interface 98 
element 

index 35 

ofa list 35 

ofa matrix 73 


108 


Index 
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equal 68 

equals 33 
assumption 68 

equations 33 
changing the subject 37 
left-hand side 33 
manipulating 36 
plotting 38 
rearranging 37 
right-hand side 33 
simultaneous 37 
solving 33 
solving numerically 50 
trigonometric 49 

erf 66 

evaluate 
an expression 11 
force command to 64 
using command-line interface 98 

exp 45 

expand 32 

expanding brackets 32 

exponentials 44 


factor 32 
factorising 32 
facts see assumption 
facts 68 
file types 27 
findroot 50 
float 17 
using command-line interface 103 
forget 68 
forget assumption 68 
fpprintpres 24 
using command-line interface 108 
front-ends 7 
fullratsimp 32 
functions 41 
delete 43 
trigonometric 47 
trigonometric, inverse 47 
functions 43 


gamma_incomplete 67 

gnuplot 56 

graphs - 
using command-line interface 115 
colour 39 
legend 39 
plotting 38 
vertical range 39 


help 29 
command-line interface 113 
horizontal cursor 20 


ident 77 
identity matrix 77 
imaginary numbers 34 
imaginary part of a complex number 88 
imagpart 89 
implicit_plot package 55 
indefinite integral 64 
index (of a list) 35 
inf 85 
infinite sequence 83 
infinity 85 
input line number 11 
using command-line interface 98 
input prompt 11 
inputting matrices 70 
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